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IN SEDIMENTATION | 


By H. A. -Einstein,1 M. ASCE 


‘The answers to the question of which direction 1 the ‘research ina given 
- field should go will depend, to a large degree, on the answerer ’s attitude. If he 
the attitude of a pure scientist, he direct the research towards the 
basic understanding of the processes, towards the analysis of the problem in | 
_ general terms. If he is more of anengineer, he will first find out which phases 
“of the problem are likely to become important in the future for the solution 
of the most pressing practical problems. He will then work on these phases, _ 
developing theories only as far as necessary for a most efficient of 
these problems. The engineering approach is adopted in this paper. 


The ‘most important ns to have taken place in the — mae since ai 


~ about 1910 that have had aninfluence onthe sediment motion in our rivers are | te 
Bead follows: : 1. During this period, the United States has essentially changed — 7 
from a predominantly | agricultural to a predominantly industrial economy. 
2. During the same time, the total population and the water requirements per ¥ : 


person have increased to such a degree that in almost all parts of the nation © 


= naturally available local year-round water supply has become insufficient. | 

_ Both of these developments, the first one by its quest for energy and pro-_ 
‘a water and the second directly by its inc reased water demand, have led > 
to a large-scale development of all available fresh-water resources. The key 
to this development is the provision of storage space in which the available — 

o a —Discussion open until August 1, 1961, To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE, This paper is part | 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American - 
_ ciety of Civil Engineers, Vol. 87, No, HY 2, March, 1961, =e ae 

Prof. of Hydr. Engrg., Univ. of Calif., Calif, 
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water can be stored during times s of 0 over- - supply y for cc omnauanien during the 
_ dry periods. Most of the artificial storage is provided in the form of surface 


reservoirs. These reservoirs rs change the regime of the rivers drastically, 
particularly with respect to their sediment 

1. Any large storage reservoir will permanently store practically the the 
entire sediment load of the 


2. Part of the sediment storage occurs in the storage volume proper, ‘part 
in the channel and valley bottom upstream by backwater effects. jj | 


“a The sediment load leaving the reservoir is usually nil or is restricted > 
4 a very much reduced releases that have the purpose of flushing some of the 


sediment through the storage volume. 


4, ‘The hydrograph of flows “available in the stream channel downstream 


from the reservoir is entirely changed and with the annual sediment load 


in 1 and Above the Reservoirs.—Sediment deposition inthe reser- 
--yoir causes loss of storage volume; deposition above the reservoir causes | | 


_ back-water effects, swamping and flooding in the valley. These damages have © 
hitherto been accepted as unavoidable and toa large degree as uncontrollable. 


Many examples of actual reservoirs and lakes have shown that the deposition _— 
q of the sediment may follow many different patterns. Often, changes of patterns 
are clearly | caused by particular influences, such as snags or vegetation, that 
suggest possible engineering control of the deposition process. If it is ‘realized 
that both the storage volume in the reservoirs and the area of the extended 
valley bottoms above the reservoirs represent large and often irreplaceable © 


values, it appears to be extremely wasteful to accept these losses, as they 
may occur, without question. A planned program of valley sedimentation, on 
- the other hand, may not only reduce the loss of storage space in the reservoir 
but also create large valuable valley areas that are sufficiently high to pre- 
_ vent flood damage if used agriculturally or industrially. Such a program may > 
_ also prevent the development of swamps with their health hazards. pone aie 
. What type | of research is needed in connection with the _ development of 
such areas for sediment deposition? The first stepis probably the systematic _ 
7 observation and description of the various deposition processes as they occur | a7 
in existing prototype conditions, - covering all materials from the bed load to > 
_ the finest, for an extended range of river conditions. However, even such a 
_ very extended survey is only able to describe various chance combinations of __ 
the many variables that may | or may not be important in their effect on the _ 
7 final result. Only a a systematic study of the problem under the controlled condi- - 
> 


tions of a laboratory will permit the separationof the effects of the individual dual 
Such laboratory studies, in turn, become effective only after the similarity a 

conditions for deposition processes are known. Today, most model studies © i 


with movable bed deal with an approximately uniform bed material and try to | j a 
7 find the geometric patterns of scour, transport, and deposition for that ‘material ; 
inclusion of all the fine components of the load in order to study the more - : 


_*< given flow conditions. The study of reservoir deposition would need the | 

After” the development of methods operating with a large range of grain 
sizes under deposition conditions, and after scale effects for such processes 7° 


known, one can com specific questions. 
=e 
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_ Under which flow conditions are s sediment « deposits segregated : as to aay 
sizes and under which are mixtures 
What are. the relationship; 
7 What are > the > agricultural values of the v various types of deposits? +2 


flow is not able to produce? 


These are some of the most pertinent questions in in connection with reser- 


voir and valley ey deposits. questions may arise in connection with the 
engineering control of large-scale deposition processes and may be answered 
studies. It i is certain that many more such problems will arise 


The entire field ‘ot engineering control of deposition processes in and 
above reservoirs is practically unknown, today. It is unquestionably possible 
: to develop rules, methods and other tools for that purpose in the same way 
that similar know-how has been developed for the deposition of sediment in 
- near-equilibrium streams for the purpose of channel stabilization. apa 
= There are ‘mainly » two reasons for the lack | of such dev elopment. First, 
most of the reservoirs of significant size are located upstream of the area in 


which, at the time of valley development was important. 4 


tion or two to warrant the expenditure of funds and prs to prevent it. That 
means that the incentive is missing in cases in which new areas are developed. — 
Secondly, where reservoirs” are built into an already developed area, the | local © 
landowners above the reservoir do not benefit from it, but usually incur a loss 
of usuable land in various ways, and always feel that they are not mn el 
= fors such losses. The operators of the reservoir, private or public, ; 
are then in the precarious position n of having hostile neighbors in the saul 
aren in which such controlled deposition may be practiced. = 
It is ne difficult for an individual to prove legally that he nw ili 
damaged by very slowly progressing "deposition process. In ‘contrast, 
engineering measures are undertaken to control deposition in the valley bot-— 

_ tom, many local owners who may actually benefit from the operations, will go 
to ) court and claim damages. Thisisa rather costly undertaking for all parties” i 
ge It can 1 probably be avoided by o organizing all landowners into a district 
that can then as a unit negotiate with the operator of the reservoir, aftera _ 
> concrete plan has been established for the operationo of the deposition process. = 
Unfortunately, ~ sufficient general knowledge. is not available today on these 
processes on which such a plancanbe based. The economy of any such plan is — 
usually even more in question. Only specific r research i in this field can can provide 7 ; 

the necessary background for such developments. 

_ Sediment Problems Downstream from the Reservoirs.—The construction of 

* reservoir will have its greatest benefits - in the river valley downstream | 

through flood control and increased all- year | water supply. The river valley — 

downstream must be expected to develop economically at a much increased — 

rate whether it remains agricultural or changes to urban and industrial de- — 

velopment. The river that brings the fresh water from the reservoir will also” 
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be called on ate drain the valley bottom and, depending on its size, may even | 
provide inexpensive transportation. Because the river channel is — 
provided. the laymen are apt to accept its shape and course as given and 
stable. _ When the majority of the high flows are held back by the reser reservoir, 7 
most or all of the flood danger : appears to be eliminated. Te 
To the engineer who is familiar with the management of alluvial rivers, © 
the picture | looks less rosy. He recognizes that the construction and — 
of the reservoir severely reduces the sediment supply and radically changes 
_ the hydrograph of the river. These changes are so drastic that in most cases 
7 there is very little similarity between the new and the pre-reservoir river. 
Only the channel is the same. It is a sediment-lined channel that had usually 
reached at least near-equilibrium shape for the old river condition. It would 
be an extreme coincidence if the channel were still in equilibrium at the new 
flow and sediment conditions. The sediment supply i is reduced. But also the 


p 
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4 flood flows and the increase of medium flow durations. Usually the sediment . 
Supply is” reduced more than ‘the carrying capacity, causing the channel to — 
erode. Often this development may be reversed in lower reaches after some 
sediment carrying tributaries have joined the main stream with the reduced | 
sediment capacity. The immediate response of the river will be scour in the - 

= reaches. It will have acquired a fullload when it arrives at the junction 


with any sediment carrying tributary. It may beginto aggrade at and below that _ 


junction. ‘If the river channel is left free to develop at will, the upper reaches» a i 
L will, ‘after a given time, adjust themselves to the new conditions. ‘They w will, a 
stabilize at a lower sediment rate, andthe sediment load to the lower reaches 
will be reduced. It must be expected that during that process of adjustments 7 
the sediment supply to all parts of the river system will steadily « change and 
- After a channel section has finally stabilized its profile, for instance by > 
anmaie its bed sediment, it is, unfortunately, not yet completely stable. 7 
The banks, that usually consist of a finer material than the bed, also begin - 
-y fail. They also seem to be stabilized only as long as finer material is — 
deposited from the stream at about the same rate. at which it scours. ‘With 
the finer sizes missing in the load, the banks retreat steadily and the river 
eT widens. The flows begins to meander between banks and attack them = 
with increased force where they impinge on them. ~ | is not si sufficient to sta- 
 pilize the bed alone; the same must be done with the banks. _ 
_ The valley bottom below a reservoir has been described as the area in 
; which the maximum economic development must be expected to take place. 7 
_ The valuable ground will be in great demand and developments of all kind 
: will: soon begin to encroach on the river. A gradually enlarging river channel 
that is usually already too wide for the reduced flood flows is, therefore, a 


Under these changed flow conditions, the engineer will ask again a a 


questions when be tries to design a more desirable channel below a major 


methods may to stabilize the profile of various 
+ What methods may be used to stabilize the banks of various : ene — 
a What cross section | _ be stable for such a changed ‘stream’ 
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"SEDIMENTATION 
"What changes will oceur in compo ion of the 1m bed. due to the 
Is it possible to. stabilize channels for a range of a supply rates 
Is ‘it more advantageous to stabilize the river channel | before, immediately 
after, or any time after closure of the dam? 
What types of structures are most effective for the of 
channels and which materials should be used in their design? 
a Should stabilization ‘begin at the dam and progress +2 or should — > 
it begin simultaneously over a considerable length of stream? | 


Should the changed stream after stabilization follow a meander pattern 


similar to the original stream or may it be straightened? = a 
Many more such questions can and will be asked, particularly if ‘special — 


“interests such as navigation, water quality control, and others are involved. — 

Now one must translate these engineering questions into research problems. a 

‘This has already been attempted at the University of California, Berkeley, — 
Calif, sediment laboratory for a number of years. Answers have been obtained 

from the e resulting studies that, in many instances, point the Way) to the solu- | 

= tions of practical problems. First, the general conclusion was drawn from the 
- previously: stated arguments that any designed and stabilized channel must 
have artificially The width between these banks is not auto~ 


z= contrast to the alluvial - river in its natural equilibrium state, in which : 
the width reaches an equilibrium value depending on omen of the flow, 
sediment load, the vegetation and other factors. So 
_ The basic problem of the changed river as it arises in its most general q 
form can then be stated as follows: the and the sediment 


: depths and flow velocities will result as a long-range equilibrium fora given 
chosen — channel — width? Or the slope may be given and the channel width © 


= But a a first step to. its solution has been | made a at the 
University of California, 
_ An | equilibrium flow with a very heterogeneous load was disturbed by both 
an increase and then by a decrease of the sediment load. It was found that an 
- increased load will behave quite rationally as may be expected from trans-_ 
portation laws for uniform sediment, The decreased load, on the other hand, a 
resulted in a . strong « coarsening of the bed with progressively reduced scour. ¥ 
a it was found that for each flow there exists a range of sediment sizes that 
“move more slowly than the bulk and that all sizes above a limit practically 
i a = not move at all. These “slow- moving” and “non-moving” sizes appear to 
5 be the tools that nature uses to retard and even to prevent excessive down- | 
7 ward scouring of channels. This process of bottom shingling and the resulting 
very” limited bed scour has been observed various rivers below large 
reservoirs. More work is needed for its quantitative prediction, particularly 
on the effect of flow variability, A, 
In many | rivers such non-moving sediment sizes are not contained in the a 
river bed in sufficient amounts to limit down-scour to an acceptable depth, 
} How can a similar protective layer of loose stone be introduced into the river | 
bed artificially to limit down- scour? ‘How does the necessary amount of rock — 
depend on the rock size? These are all most important problems that can be : 
‘solved and for which answer does not exist. 
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In river reaches in which there is still an ample sediment supply in the 7 
bed sizes, but a lack of fines, the banks may need ‘Stabilization whereas the 


nk bed is still in fot alalinge The question arises in such rivers how deep the 


‘foundation of the banks must go. In any event must the toe of the bank be be- a 
low the deepest instantaneous scour holes of the bed that may develop near the 


bank. This calls for a systematic. and quantitative study of the shape of bed 


4 


plex problem has recently found some interest in various quarters, but no 7 
“a is developed as yet by which quantitative river predictions can be made, 
_ These are some of the most important developments that may be expected 


7 irregularities for the various flow and bed conditions. This extremely com- © 


to occur in our rivers and an attempt has been made to predict the mostim- 


portant research problems | that will arise with them. _ It is most interesting to- 
remember that many of our notorious sedment carriers were clear and deep, 
ae flowing, and tree lined streams at the time agriculture was brought to 
land. changed into what they are by reckless methods of soil 
exploitation. it may now be necessary to change many such streams back to 
“their original state because their waters have become 


to conserve them. 


reservoir is seen as the most important cause for new 
sediment problems is our rivers. Sediment transportation research should 


concentrate on the study of river conditions aren ae and downstream from 
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ALLOCATION OF WATER IN CALIFORNIA 


By Amalio Gomez, -F. AS 


Pst 
SYNOPSIS 
This describes briefly the the costs- remaining benefits 
7 of allocating the cost of multiple-purpose water projects; summarizes the © 
allocations made by the Corps of Engineers for several reservoir projects in 
7 the 2 Central Valley of California; compares the results obtained; and makes 
some observations concerning the of the method, based on actual 


NTRODUCTION 


functions that the p: project. serves, It must a remembered that no fully | satis-_ 
factory procedure has ever been devised for assessing taxes equitably in all 
“cases. Similarly, no fully satisfactory procedure has ever been devised that Jf 
will automatically yield equitable results in allocating the cost of a multiple- 3 
purpose water project, in all cases. The problem, therefore, concerns the 
computation of approximate values that reasonable men consider fair and just. a 
Sound and ‘experienced judgment is still one of the basic ingredients | of any 7 


_ Note.—Discussion open until August 1, 1961. To extend the closing date one ah, 
a written request must be filed with the Executive Secretary, ASCE, This paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the > American So-- 
ciety of Civil Engineers, Vol. 87, No, HY 2, March, “1961, 
1 Chf., Planning and Reports Branch, Engrg. Div., U. S. Army Engrg, Dist., Sacra- 
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py tself. It probably began with the @ 
“| ei Of Siaring we cust Or Making war or Of running the government. It 
- be said that cost allocation is a new name for describing a process quite : 
Similar to that of assessing taxes. It is now used primarily to describe the 
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back some 30 yr. Many methods have been Geviect and used with more — 


less results. of the most co common methods used | include: 


Proportionate use method 
Priority of use method © “— 


Alternative justifiable expenditure method 
‘Separate projects method 
Equal apportionment 


ethod in use by Federal agencies is the costs- 


m 


7 


TABLE 1,—COST ALLOCATION 


4, Sey orale costs 
‘Separable costs 
5. benefits (items 4) 


Allocated residual cost a’ 
1. Total allocation (items 4+ 6) 


project functions share in the savings to be realized by the 
- multiple- purpose development. In the application of the method, project bene- — 
_ fits for each function are first limited by the cheapest alternative cost of ob- 
taining the same benefits. ‘The method then assigns to each function all = 
“separable cost,” that is, the costs that were incurred by virtue of adding each ' 
_ function to the project. _ The remaining costs are called “joint costs,” or 
“residual costs,” and are distributed among the various project functions in 
proportion to the benefits remaining after deducting each separable cost from _ 
each project benefit. eae method is best described in a report to the Inter 


House Committee Print No. 23, 82nd Cong., 
flee, Washington 25, D. C. 
7 


2 “Report on the Allocation of Costs of Federal Water Resources Development Pro- 7 


2nd Session, U.S, Govt, Printing Of- 


Separable cost-remaining benefits method 
3 
ow Flood | Power,| Irriga- | Naviga- “Total, 
in ‘| 1,000 | in’ | in’ | $1,000 
$2,000 $1,000 | $2,000 | 
i. Benefits |= | 500 | 1,500] 350 | 100 | 2,450 
400 | 1,000] 600 | 80 | 2,080 
1,000 350° so 1,830 
20 so} 200 | | 
— ii 
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COST ALLOCATION 
Agency Committee on Water Resources.3 This report was prepared by the © 
Sub- Committee on Evaluation Standards, that developed the separable costs- 
benefits method. A fair understanding of the method can be 
from the following quotation and Table 1 taken from the repost: peal 


— 2 recommended method of cost allocation is: illustrated below 
for a multi-purpose project for which the total project costs — | 
to § $1, 765, ,000. These include investment costs and operation, main- = 
tenance, and replacement costs, ine to a common time basis, 


— 


of the Corps of Engineers is to preliminary 
allocations in connection with preauthorization reports, The preliminary cost. 


allocation is brought up-to-date when funds for initiation of construction are 
being requested, and again as soon as after construction is initiated. 


Folsom i Completed 


Success Tule River Under construction 


Calaveras may er Under construction 


‘The final cost prepared pest the is placed in opera- 
7 ~ tion. The _ practice followed during the last few years varied somewhat with 
each project, but, generally, it was to make a preliminary cost allocation 
_ before initiation of construction and fix the percentage distribution of project 


: cost among its various functions at that time. Later, these percentages were 
_ applied to the actual construction costs. 


& 


Multiple- built by the Corps of Engineers 
Central Valley of California serve two or more of the following functions: _ 
(a) ‘navigation; (b) flood control: (c) irrigation; + (a) municipal and industrial 
ae supply; (e) hydroelectric power; (f) recreation; and (g) fish and wildlife. ; 
- Cost allocations for Corps of Engineers projects are made by the Chief of 
_ Engineers, generally with the concurrence of the Bureau of Reclamation ‘and 
; other affected Federal, State, and local agencies, With some impor tant ex- 
a 8 “Proposed Practices for Economic Analysis of River Basin Projects,” U. s. Govt. :? 
Printing Office, 25, D. C., May, 1958. 


' 
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NSTATUS 
‘Black Butte 
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HY 
ceptions, costs allocated to navigation and flood control are e borne by th the 
Federal Government. Costs” allocated to irrigation are paid for by the water 
users, generally without interest, over a 40-yr period, under a contract with © 


; the ete of the Interior, made pursuant to the provisions of reclamation — 
TABLE 3.—SUMMARY OF ECONOMIC _ 


= Pine Fla t, Kings F Isabella, Kern River 


_Flood-control storage 

Average annual runoff 1,6 729° 


Average annual new water 
Ave, ann, new and rereg. "water 


Interest and amortization 39000 

Total annual cost 


benefits 
Power benefits 

Total annual benefits 


First 


Flood control «15, 


4 573 


1 


21, 093 


| 


a $590,0 000 $59, 000 acre-feet of municipal and industriel water supply. 
b After deducting $2,140,000 creditable to generating, switching, and tr ansmission 
facilities built by the U. S. Bureau of Reclamation, the cost of which is not included in _ 
the costtobe allocated, 
© Includes $8,810,000 allocated to municipal and industrial water supply. = 
2 Exclusive of cost of power facilities built by the Bureau of Reclamation. . 
Includes $26,000 allocated to municipal and industrial ‘water supply. 


law, including the 160-acre limitation and other restrictions. . Costs allocated © 
to municipal _and industrial water supply are repaid with interest, under a 
contract approved by the Secretary ¢ of the Army. - Costs allocated to hydro-— 


peng power are repaid with interest. The power is sold by the Secretary 


| g 
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“Until recently, recreation and fish and wildlife been considered as 
byproducts of the primary project functions. How- 


ever , recreation and fish and wildlife are ‘becoming more important every 


DATA AND ALLOCATED costs 


Folsom, Ameri- Success, Tule | Terminus, Ka- Hogan, Ca-| Black Butte, 

can River River weah River | laveras River Stony Creek 
(5) 


65,335 22,350 


=e 


404 19,190 11,238 11,000. 
* = 3,160 6 7,300 
= 
18, 18,300 


1958 (P.L. 85-624) provided that in certain oe projects and in all new : 
projects, fish and wildlife will be considered for possible inclusion as a pro- 4 
ject function on the same basis as any other project f function, and that the cost 
allocated thereto may be considered wholly or in part -non-reimbursable. 
the projects discussed here, no project costs have been allocated to recrea- 


| 
18,300 
5,350 
Cost, in $1,000 
4 6,765 
 42,180° 
65,3385 
1 
gt 
— | 


- multiple- purpose reservoir ‘projects in the Central Valley of California. 
des: names of these projects pao the construction status are listed in Table 2. 
"The Pine Flat cost allocation was not based on any specific method of cost 
gave weight to the results obtained _ by the use of several 


4. .—COMPARISON OF COST ALLOCATION TO IRRIGATION. 


2 
Pine  Isa- Ter- ew 
| Flat | bella | Folsom|Success| minus| Hogan | Butte — 


Allocated first cost per | | 
acre-foot of average annual 


runoff, in dollars —_ 8. 48 | 6.27 


-acre-foot of newandre- | | 


i. 
Allocated first cost per. a 


acre-foot of reservoir — 
storage usable for 
irrigation, in dollars oo 


& 


sf Annual cost to water user 

on basis of repaying 
cost in 40 equal payments 7 
without interest but includ- 
ing M. O, and R, charges: 


acre-foot of 


_ 

uo 


Per acre- foot: oft new 
reregulated water, in 


dollars 


me thods. all other have been. made by the separable costs- 


~ remaining benefits | method. Table 3 contains some descriptive data on the use - 


of storage space as well aS a summary y of first and annual costs for or each — 
~a- It also contains a summary of allocated first costs and maintenance, _ 
Operation, and replacement costs. A further summary fo for each ‘Specific pro- = 
ject with pertinent comments is given later, | 
It is extremely difficult to pass judgment onthe fairness of a given alloca- 
tion. For this reason, results obtained at one project are often compared with 
results obtained at other projects. Such comparisons can serve a useful pur- 
pose. _ However, they must be properly qualified because there are too ney 
unrelated factors that affect the results. Among these factors are: (a) physical 
characteristics of the dam and reservoir ‘site; ite; (b) characteristics of the 


& 


| — 
tion an in some cases when the 
me costs will be charged against these — fe! 
A 
q 
7 
q 34.58| 49.85 | 57.45 88,36) 128,07 
2,53 | 2.99 
_ 
= 
— 
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ALLOCATION 
stream such as runoff and extent to which already developed; . (c) service are: area 
7" in which water is used; (d) developments inflood plain; and (e) the price level — a) 
_under which the project was constructed. . With these qualifications in mind, — 
the comparisons involving the cost allocated to irrigation given in Table 4, 
may be of some interest, 
me preceding comparisons have been made on the basis of cost allocated 


to irrigation, because such costs are repayable and there is a tendency to 
_ judge the fairness of the allocation on the basis of sapayene costs. | 
— 


"PINE: DAM RESERVOIR 

The cost allocation for the Pine “Flat Project4 was made in 1947. It was 
‘based on an estimated first cost of $35,760,000, of which $33,500,000 was 
dam and reservoir, and $2,260,000 was penstocks for future 

power at the site. The allocation was as follows: 

Flood control $18, 250,000 


2,260,000 
$35,760 ,000° 
2 The War Department Civil , Appropriation Act, 1947, appropriated funds for 
the initiation of construction of this project subject to the proviso that: 


of the appropriations for the Kings and Lake 
project, California, shall be used ... until the Secretary of War. . 
_ with the concurrence of the Secretary of the Interior, shall have made 
a determination as to what the allocation shall be: provided further . 
that the agreement of concurrence shall be made not later than 9 months 


the date of the enactment of this Act.” 
Because of the increase in price level and some design changes, when the © 


dam was finished in 1954, the total cost was about $39, 385,000, of which about — 
$1,000,000 was for penstocks for future power. Because the allocation had 
been | made at the specific request of Congress, and it had been furnished to 

and “accepted by Congress, it could not be revised or done over again, on the 
_ basis « the actual construction cost, without a new directive from the Con- 
gress. This has not been r with any of the subsequent pro-— 

The project is for irrigation on the basis of interim con-_ 
tracts between the water users and the Bureau of Reclamation. A permanent 

oct is expected in the near future. The project has been eminently suc- 


made 


cessful as illustrated by the following accomplishments during the 6-yr period 
since its completion in 1954: 
payments received from water users. $4,627,000 
Cash payments received for s storage of power ‘water. - $704, 000 
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— ISABELLA DAM AND RESERVOIR: i 
The cost allocation for Isabella Dam and Reservoir was made in 1955, 
when construction of the project was essentially completed. In this allocation, 


annual cost and benefit data were capitalized on the basis of 2.5% Federal — 
- interest rate and a 50-yr amortization period, andthe allocation was made on 7 
a “present worth” basis. Table 5 contains a summary of the allocation, 
> Even though no power generation at the site is involved in this project, the 
7 regulation provided by the project does benefit existing downstream power 
plants, and a cost allocation commensurate with the power | benefits was made. 
_ The project has been operated for irrigation since April, 1954 under terms 
7 of an interim con! contract between the water users and the Corps of Engineers. 
5. —SUMMARY OF THE COST ALLOCATION. 


_ Flood control, Irrigation, a. in 
$1,000 ~ | in $1,000 $i, 000° 


«38,374 7 403 47,479 

costs | 
Separable costs 

Remaining benefits 

Allocated joint cost 


allocation 
Allocation of first cost 


4 permanent contract is. being r negotiated ed between the water users — the 


United States Bureau of Reclamation (USBR). Payments 1 received so far for : 


the interim irrigation se service total $1,134,000, and flood damages prevented 
to date amount to $2,500, 00¢ 000. 


FOLSOM DAM AND RESERVOIR 

Folsom [ Dam was built by the Corps | of “Engineers. Folsom power plant and 
‘the Nimbus dam and ‘power plant were built by the USBR. The entire project - 
‘is being operated by the — as apart of the Central Valley Project. A pre- © 
liminary cost allocation for Folsom Dam and Reservoir was made by the 
Corps of Engineers in 1953. The results of this allocation have been used by — 
the USBR in various repayment analyses of the Central Valley Project. The — 
preliminary — cost allocation is given in Table 6. Note that the power benefits | 
and the cost allocated to p power are relatively small. The reason is that the 
*s power benefits were assigned to the power plants and to the dam and reser- 
; = such a manner | as to provide for an equal: ratio of benefits to cost 


credited to the dam and reservoir are shown in Table 
‘iin An interesting ‘sidelight in connection with this project is that no sooner ~ 
the dam completed the greatest flood of record, that of 
_ December, 1955 occurred, and the Sacramento metropolitan area was spared 
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— 
- 
q 
— 
a 
r 
— 


COST ALLOCATION 
a great cialis It has been estimated that during the 1955 flood, the dam 
_ prevented about $20,000,000 in flood oe that is more than the total cost 


were 


allocated to flood control. 


_ SUCCESS DAM AND RESERVOIR > 


A for the Success Project was made in 1956. 


Mei _ The percentage distribution of the project cost between flood control and irri- 


TABLE 6. .—PRELIMINARY COST ALLOCATION, “FOLSOM ‘DAM 


Flood 

_ control, in $1, 000. 

@i@ ii 


benefits 
Alternative annual costs _ 
Annual benefits limited by 

alternative costs 

Remaining annual benefits 

"Allocated joint costs 
Total annual costs 
O, & R, costs” 

Interest and amortization 

Firstcosts 


2.475 
65, 


; gation was fixed at that time. Such percentages will be _— to the —e 
construction cost when the project is completed. ‘summary of of the cost 
allocation is given in Table 7 (a). 


mig A contract is now being negotiated between the USBR and the water us sers” 


- for the amount to be repaid. The project: is now under construction and ad 


- scheduled for completion in 1962. 


_ * preliminary « cost allocation for Terminus Dam and Reservoir was" made 


| flood control was fixed at that time. These percentages will be applied later — 
t to the actual construction cost. . A summary of of the actual allocation is given 
A is por negotiated between the water users and the USBR for 


NEW HOGAN AND ! RESERVOIR 
a A preliminary allocation of ‘cost for the New Hogan Dam and Reservoir 
was made in 1959. This allocation fixed the percentage distribution of the pro- — 


4 
AND RESERVOIR 
1,610 | ~—-2,160 4,720 q 
1,530 | 1,400 | 400 3,400 
750 | 1,370 360 | 320 | 2,800 
4 «16,765 | 33,370 | 8,810 | 6,390_ 
> 
a 
= _TERMINUS DAM AND RESERVOIR 
under construction and is scheduled for completion in 1963. | 
a Re 
| 


‘ject cost between ‘flood control and irrigation, but the amounts will ‘not 
finally determined until the construction costs are known. A summary of the 


is interesting to observe that construction wore appropriated for 
- fiscal year 1960, with the stipulation that construction not be initiated until a 
_ firm irrigation repayment contract had been executed. Such a contract between 
- the State of California and the USBR was executed in March, 1960, and actual 


(a) Success Dam and Reservoir 


| 


Irrigation 
in n 000 


_ Alternative annual costs © 
Annual benefits limited 
parable annual costs (0 
7 Remaining annual benefits 
Allocated joint annual { 
Total M.O. and R, costs 
Interest and amortiza- = 
First costs 
Percentage of total costs 


the prime objective of permitting initiation of construction at an early date. 
The contract allows the USBR seven years to negotiate directly with the water 
users for the repayment of the amount allocated | to If the 


BUTTE DAM AND D RESERVOIR 


a 


po flood control and water conservation. Water conservation includes irri- 7 
- gation and possibly municipal and industrial water supply. i The fixed percent-| 


ages will be applied later to the actual construction costs. A summary of the 


ae in the case of the New Hogan Project, funds appropriated for fiscal 
ar 1960 could not be utilized until a firm repayment contract was executed. 
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: Such a contract was executed in March, , 1960 ‘by the USBR and the State of of 


California. The contract is identical to. the described for the of 

FUTURE PROJECTS 
“a In addition to the projects described heretofore, a cost allocation has just 
been ‘completed for the Oroville Project of the State of ~The 
ow 


Total, | Flood | Water 
in in am, 000 


_ tion in the interest of flood control, sucha contribution having been authorized 
by Congress” ‘subject to the amount being fixed by the Secretary of the Army 
and approved by the President. -é preliminary | cost allocation has also been 
made for the proposed Hidden Reservoir on Fresno River. The reports con- | 
taining these allocations are now being considered by the Chief of Engineers, 
If approved, the results of the allocations will soon be made public, pee 
Preliminary cost allocations are also being made in connection. with pro- 
spective multiple- purpose storage projects on _ Stanislaus, Tuolumne, Chow- 

chilla, Merced, Mokelumne, and Yuba Rivers. 

_ These possible future projects are noted enue to illustrate that the sub- 
ject of cost allocation in the Central Valley a ecuee is important and 
bears on many current ar and f future 


A ‘superficial examination of the costs- remaining benefits method 


4 may lead to the saaeeemeemee that there is finally available aspecific and definite 


| 
| 
| 
: 
— _ and Reservoir} (c) New Hogan Dam and Keservoir | (d) Black Butte Dam and Reservoir _ | a 
Total, 
899 = 258 188 312 | 500 
105 93 | 57° 81 54 135 
711,213 6,362 | 17,575 | 11,000 7,300 18300 
— | 63.8% 86.2% 60.1% 39.9% = 
_ 
— 

— 


a 


/ method that will will produce the same ‘the same results no matter who ‘is making ‘the cost 
allocation. — Far from it. The x method assumes that a great deal of data are 
available that were gathered ‘during the project formulation stage. Survey 7 
_ report investigations, that include project formulation and economic feasibility 7 4 
- ‘studies, are carried on generally with limited funds and often short cuts are —_ 
“needed. Therefore, it is seldom that all the data required for cost allocation © 
purposes are available. Even the data that were originally available have be- 
come partly obsolete due to the time elapsed since they were collected and — 
due to the change in n economic conditions. The data most frequently deficient or 
unavailable are the cost of the cheapest alternative projects and the separable 
costs. The development of these data at reasonable cost requires the exercise > 
of mature judgment. The most difficult part involves the determination of the 
magnitude and location of the physical works required to obtain the specified _ 
_ services in each case, with reasonable assurance that the combinations 
selected are the most economical. ‘Normally, this is a time-consuming trial 


and error job. Hence, short cuts by experienced men are almost a practical 

yr _ However, in spite of these drawbacks, the author believes that, as a general 
eule, the ‘separable costs-remaining benefits method is by far superior to any 

other method with which he is familiar. It also has the advantage that it is : 


Bi nog to be better understood ar and generally accepted. The allocations made 


by the Corps of Engineers in the Central Valley of California, by using this 
method, have been generally accepted as being fair and equitable by the water 
_ users and by the organizations involved, "Ww : 
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‘VECTOR ASPECTS OF DYNAMIC SIMILARITY __ 
By R. C. Kolf,! M. ASCE, 


SYNOPSIS 
A study of ‘the lowe | of “dynamic similarity emphasizing tt the ‘refinements 
imposed by the vector nature of the "fluid flow variables is presented. As 
examples of the method proposed, the problems of turbomachinery and geo- 
ea metrically distorted models that have needed separate treatment in the past, 4 
‘because of their eee distortions, , are shown as logical extensions of the 
"same theory. — 


INTRODUCTION 


The ba bases for the prediction of the characteristics of flow in a fluid me- 
chanics design problem are usually grouped in three general categories: (1) 
That of “engineering experience,” (2) the laboratory method of studying each 
ed by means of scale seedela, and (3) the process of theoretical analysis. 
In most of the problems that face hydraulic engineers, the shape of the physical a 
‘boundaries that confine the flow are complex and they are often in motion. 
Consequently, the prediction of the flow characteristics by purely analytical — 
means is many times impossible because of. the lack of presently available 
a theoretical methods that are adaptable to such | situations, The second method | : 
of prediction from the study of hydraulic models, utilizing whenever possible 
the best features: of the. other methods, has proved very useful in such situ- 


‘Note. .—Discussion open until August 1, 1961. 2 - To extend the closing date one month, ; 
a written request must be filed with the Executive Secretary, ASCE. This paper is part _ 
_ of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- in 
ciety of Civil Engineers, Vol. 87, No. HY 2, March, 1961. 

1 Assoc. Prof. of Theoretical and Applied Mechanics, Marquette Univ.. a , Milwaukee, 


2 a. of Theoretical and Applied Mechanics, Marquette Univ., is Milwaukee, Wis. 7 
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Notation. —The letter for use in this paper 


a 
are 


When only geometric similarity is of the prediction 


equations can be written in terms of a single scale ratio, in 


‘ 2, and 3 V is the volume, L refers to the characteristic Tength, 


A second class of variables involves time measurements as well as" 
length measurements. . For kinematic similarity a second s scale ratio, namely — 
Pe the ratio of characteristic times, tr, must be established. — - 


Vp 


in which V V denotes the velocity and a is the acceleration, 
A more complete ‘listing of these prediction ratios is available in many 


f standard watemnes works,3 The significant characteristic to be noted at this 


absolute magnitudes of vectors. If this is to have any pertinence 
in terms of similarity, it is necessary that each velocity or acceleration - 
s vector in the model have the same direction as the corresponding vector in| 
s @6=—séthe prototype. Then, in addition to Eqs. 4 and 5 being fulfilled, the flow pat- 
3 terns must identical. This condition is ordinarily fulfilled, or at least 


approximated, insisting | on ‘geometric similarity as a prerequisite for 
_ ‘Kinematic similarity because, for each boundary geometry, the ideal flow 


ASCE Manua 


alphabetically, for 

— 

of dynamic similarity are the criteria used in the design and 

™ | analysis of hydraulic models. Depending on the type of model to be built and — 4 

a 

— 
— 

3 — 
a _ by another, as indicated by these vector ratios, has no known physical sig- m ane 
— 
— 
— 
| nay 4 


net is defined. of the occur because of. shear 
ae must be accounted for by the further conditions of dynamic similarity. » 4 
Dynamic ‘similarity requires that the ratios of all homologous forces in the 
_ two systems be the same. Limiting this study to cases in which only the — 
d _ pressure force, F,, the viscous force, Fy, the gravity force, , Fg, and the - 


inertia force, Fj, are significant yields 


tios of the absolute cdesibines of the vectors involved, and that the prope dell 
vectors | on model and prototype must have the same directions in order to be 


_ 4 similar, it can be seen that Eq. 6 is a limiting solution of Eq. 8. Other pa 


tions of Eq. 8 are conceivable, but to define similarity of the resultant (in- 
-ertia) force without specifying similarity of the component forces would re- 
- gult in a distortion of the geometry and the kinematics of the systemina | 
manner that has thus far been undefined. There is no present utility of such — 


in which wo is the weight, y refers to the specific weight, p ) indicates the 7 
pressure, u denotes the dynamic viscosity, V is the velocity vector and p re- 
fers to the mass density. Eqs. 9, 10, 11, and 12 in the final form shown, all 
involve a proportionality constant. For instance, in the case of the gravity 
force, F, y In order that the force ratios in Eq. 6, have the same value 
regardless of the particular force studied, these 
should be the same for the e model el as for the protot 


— 
or, because the resultant for led the inertia force 
" 
4 
| 
| 
4 
| 
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a ‘’ the two - geometrically and kinematically similar. A hierarchy of simi- 
_ larity conditions is therefore as follows: 


phy Geometric : similarity may be obtained exclusive of the other = | 


— In order to obtain complete kinematic similarity, geometric similarity 
‘must also be present, in the ide ideal flow sense, dynamic not not 

ime _ For complete dynamic similarity, oth geometric and ki kinematic ; simi- 


larity are ncess 


dynamic similarity design criteria are fashioned rear- 


Other combinations are, _ of course, possible, but they would not be inde- © 
pendent of those shown in Eqs. 13, 14 and 15 and are, therefore, unnecessary. : 
If a geometrically similar model is fashioned according to o one of these equa-_ 
_ tions in such a way that it does not violate the other two, it will generally be 
_ dynamically similar to the prototype both in the sense of the absolute magni- 
tudes of the variables ‘and in the sense of the direction of the vector quanti- 

_ ties. This is facilitated in many problems by the absence of one of these types 
of forces as a motivating factor, that obviates the necessity of conforming to 
all of the three equations (Eqs. 13, 14, and 15). Also, because the inertia q 

- foree has been included in this study as the composite representation of a 
the other forces, one of the forces (and, therefore, one of the three equations) _ +] 
_may be considered as being dependent onthe others. Thus, a particular model ' 


design will generally neds observation of only one or possibly two of 


The vector conditions implied in the development of these laws of simi- 
az larity limit the sole use of the Froude, Reynolds, and Euler numbers to prob 
lems in which the direction of the velocity and acceleration vectors is > 
"termined | by the ‘boundary geometry. In many problems this is not the case. 
e without 
the added complication of geometric distortion, have often been handled in nal 
b past as separate problems, and not as logical extensions of the same theory. r 
ey limits of the possible extension of the general laws have not been estab-— 


me 
he 
4 
Se 
= 
— 
3 
q 
4 


=. 


has not been ‘fully 


_ An example of such a problem that has been successfully handled 2d by ex- 
- tending the Jaws of dynamic > similarity is the vortex chamber. Essentially, 
. device ‘consists of a circular or spiral shaped chamber with tangential 7 
> and axial connections . It may be an enclosed (pressure force) system# as in 
_ the Thoma counterflow brake or an open system (gravity force) as in the | 
- Portland combined sewer diversion. 5 ‘If only the Euler number and onc al 


_humber conditions must be satisfied, a particular value of the one parameter | 


be accompanied by the same value of the other ir in each system. The 
following equations then apply: 


Ons the 


in which Q is the volume rate of flow, C denotes the coefficient of discharge 7 


and H refers to the head. In the vortex chamber problem, this analysis was 
insufficient because, for a particular boundary geometry, there exist an in- 


finite number of possibilities for the direction of the velocity vectors | 
pending on mn the strength | of the vortex developed. This then is a case in which 


geometric similarity does not “necessarily fix flow net similarity. In order 
to fix this similarity of vector directions - 6 was amended to 


pression for —— force (that As a component 
‘becomes 


ins which m is the mass, u denotes the tangential component of the velocity — 
vector and K is the Kolf vortex number, This expression relates the normal 
component of the inertia force vector (and, therefore, also the acceleration 


vector) to the tangential component. Along geometric similarity 
. « ‘sf ©... 


4 “An Inv of the Thoma Counterflow Brake,” Transactions, ASME, Munich 
Hydr. Inst., 1953 


— 5 *Vortex Flow Through Horizontal Orifices,” by J. Cc. Stevens and R.C. Kolf, Trans—— 
actions, ASCE, Vol. 124, (1959, 
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conditions and the other dynamic similarity ertteria F,R, and E 


similarity 1s is achieved. In In this and 18 are 
=f(R 


Another study that falls into this category is that of hydraulic machinery. 
va similarity laws have been developed for specific use in this field. How- 
ever, these criteria (such as specific speed, unit speed, unit discharge) have 
not been represented as amplications of the previously described methods, * 
but rather as substitutions for them, thus restricting the view of the dynamic - 
similarity concept. The following presentation, on the other hand, follows ~ 
directly from the basic theory presentation. 
: te The pressure force and the viscous force should be considered in the. case 

of a centrifugal pump, with the inertia force being again included in the bond + 
as the resultant of the two. If complete dynamic similarity is assumed, it > 


should be possible to show the Euler number as a function of the Reynolds - 7 


number, but this simplification of the problem would not prove useful. Be- 
; _ cause of the variability of the angular velocity of the impeller, kinematic 
in the vector cannot be fixed by | geometric alone. 


The “speed ratio,” , that. is: the ratio of a characteristic velocity, te to 


tangential component of a characteristic velocity vector within the impeller. 

: has been used to fix this kinematic similarity. This ratio is recognized to 
be the reciprocal of K from the vortex chamber consideration. It is, there- -) 
fore, proportional to the Square r root of the ratio of < Maheaed force to to centrifugal — 

force. Using the relationship uaoND and V Q + 


in which N N ‘is ati scsienatiaes velocity and D denotes the diameter. ‘Most of the 
problems of engineering interest involve high Reynolds numbers, The pres- 
sure force is then the predominant force, and in the same way that Eq. ae 
modified to the form of Eq. 22 for the vortex the Euler number 
‘must be modified E4. 24 or an 


— 6 “The Vortex Chamber as an Automatic Flow Control Device,” "by R 


’ 
B. Zielinski, Proceedings, ASCE, December, (1959, 
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DYNAMIC SIMI LARITY 


ig. 1 shows a plot of | 2 


Eq. 26 for a vertical end-suction pump operated as a 
laboratory model, with an impeller diameter of 20-5/16 in., and a geometri- 


cally similar ' prototype with an impeller diameter of 50 in. “The speed of the 
_ model was 1267. 5 rpm and the prototype speed was 514 rpm. The character- 


istic om. ¥ V, in Eq. 26 was arbitrarily chosen as the average velocity in 
« exit pipe 


or 


‘ 


a= To the scale shown it was not possible. to differentiate the data for the two 
7 ‘pumps, the WA spread” being less than the thickness of the line. Because of the 
= 7 dynamic similarity defined by the curve, each point represents a particular 
v efficiency for this design. Curves of this nature have also been used with the 
success represent the data for a ‘single pump at various: 


speeds, 
From Eq. 20 
From Eq. 


a 
FIG.1 
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will be recogr that the Eqs. 29 through the “laws of 


homologous Dus operation” ’ of turbo machines. 


GEOMETRIC DISTORTION CHANNEL MODELS 
_ The problems involved in applying the laws of similarity to the vortex 
"gamer and to turbomachinery might be classified as resulting from a 
“kinematic distortion.” That is, true vector similarity is not assured by geo- 
metric similarity. The situation encountered when geometric distortion is 
"necessary, _ although falling in a separate category, has many of the ne same 


Ory Geometric distortion is often a necessary evil in river hydraulics. The 
‘Reynolds and Froude criteria must be observed inthis case. Although viscous 
shear is not a predominant force except at low values of Reynolds number, — 
this restriction places a lower limit on the size of the model in order that the 
flow will be well into the turbulent range, and the modifications of the velocity 
profiles due to shear will be like those of the prototype. This lower limit on 
s the size is, in many instances, unreasonable for economic reasons. When this — 
“is the the case, , geometrically distorted models are used in the hope that, although 
complete Similarity is impossible, , similarity with respect to an isolated 
f phenomenon may be achieved. If the phenomenon to be studied involves a a 
“movable bed,” the conditions of of sediment eed also ‘be: con- 


design of such models. severe have been noted, 
however, using these presently accepted techniques. This may be illustrated 
by considering the Froude number in which the characteristic length is con-_ 
— to be the depth. The reason for this is that accelerations of the fluid | 
are accompanied by changes in the elevation of the water surface. The el 
fication is evidently an assumed two- - dimensionality | of the flow with the velo- 


city variations in vertical, longitudinal planes more important to the when 


7 “Similarity of Distorted River “Models,” by H. A. ‘Einstein, Chien, Transactions, 


ASCE, Vol. ot, 1956. 
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DYNAMIC SIMILARITY 
_ than the transverse variations, This bas basic assumption can be shown to be a 7 


oti All ratios of homologous accelerations should have the same value if if kein kine- 


ap = Vis 


Let us suppose that the grain size of moveable bed for the ‘model is 


Froude conditions based on depth, so that approximately similar sediment 3 
transport conditions are produced. The most common type of geometric dis- 
tortion is to a depth ratio, that is greater than the length ratio, | 
_ Good engineering practice requires that the model and scales chosen be veri- _ 
flea by reproducing known prototype developments with similar 


hon have been observed in such model studies include® 
(1) excessive scour at bends and the beginning and end of obstructions such | 
as deflecting dikes, (2) improper location of scour, (3) excessive bar building, 
and (4) formation of eddies in which none occurred in the prototype. All of 
these difficulties occur in regions of significantly curvilinear flow. The trans- i 
verse velocity distribution in such regions (scour) is related to the normal | 
force experienced by t the ‘fluid. . The transverse currents (bar building) can_ 
also be linked to this phenomenon. If the ‘model has been designed to give — 
similarity of inertial conditions in the straight portions of the channel, the 
distorted normal acceleration will result in dissimilar conditions in the curved 
sections. Because ‘the acceleration ratio | is fixed, Eq. 35, another criteria — 
that should be observed wee 


— 


factor of 6, that is not uncommon, and that the Froude condition has been ob-. 


"Supposing tt that the depth ratio, ‘i is greater than the length ratio, ie by a 
then 


: It should be ‘noted “that, because Rr = = : Le, both of these equations c cannot = 
satisfied simultaneously. Eq. 36 will have been violated by the same factor 
6, _ resulting in| severe distortion of the transverse conditions and greatly in- = 
creasing the effects of the “normal acc eleration. Several possibilities of 
alleviating this condition exist, but there is no published record of their at- 
_ tempted u: use. _ One would be to model the radius of curvature of eevee 


= eliminating the distortion in the zones that normally cause difficulty. 

7 This type of distortion would introduce a psychological problem to the secret 
because it would not “look” like the prototype condition. It would not eliminate 


"Transactions, ASCE, » Vol. 104, 1939. 
“Transactions, ASC 
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m in Hydraulic Models,” by K. D. Nichols, 


In rivers that ‘normally flow a at depths (small Froude num- 
bers) it seems that the Froude number as computed from the depth loses its — 
significance. In such channels, small accelerations of the fluid do not greatly 
affect the depth. Here then, it would seem better to evaluate the ‘parameter 7 
using the e length ratio, thus satisfying Eq. 36 in preference to Eq. 
conceusions 
ait ~The limitations imposed by the vector nature of the fluid flow variables" 
necessitates reformulation of the statement of the laws of dynamic similarity — % 
in order to extend their utility into new realms of application. In some cases 
this extension may be accomplished by using the “Kolf vortex number” in 


addition to the classical Similarity criteria. 


2. A large class of problems in which variable amounts of “vorticity” a 


- Amposea by “moving portions of the flow boundaries (such as pumps, turbines, 
and so on) | may be included in the general study of similarity by adding this _ 
parameter (ratio of centrifugal to inertia force) that assumes like directions | 
3. Many of the discrepancies noted in geometrically distorted models 
arise from a violation of the normal acceleration requirements in zones of 


- curvilinear flow. Recognition of this of this fact indicates possible model adjustments — 
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APPENDIX.—NOTATION _ 


; The following symbols adopted for use in this paper and for the guidance . 
of _discussers, with “American Standard Letter Symbols 


| ii 
2 
4 
7 
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7 
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a a 


= acceleration 7. 
= cross- -sectional area ' 


= = coefficient of disclarge 7 


eo 
ive 
=diameter 
= Euler eel 


= force 


inertia force 


= pr “essure force 
y 


= viscous force | 


= 


- — 


volume rate of flow 

= radius of curvature 7 
= Reynolds number 

< _ = distance along streamline 


= tangential component of velocity vector | 


= velocity vector 


weight 


a 


= specific weight 
dynamic v viscosity | 
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THEORY OF WAVE AGITATION A HARBOR 
By Bernard Le Méhautél 


iw A first order theory for the whee of the \ wave agitation ina simple basin | 
_ caused by an incident periodical gravity wave is presented. This ~— is a 
natural result of experimental results and previous theory, = 
AO harbor may be considered as a combination of discontinuities separated. 
— by a definite distance. The ultimate one, often in the form of a beach, is a 
de total obstruction that may dissipate a part of the incident energy. The value 
_ of the agitation in such a basin is computed when the motion is predominantly | 
two dimensional, as a function of the amplitude and period of the incident wave. 
A practical example applies the results of the theory directly. 
_ The influence of the shape of the entrance, beach, wave traps, and so on gq 
on the value of the the agitation are vino cag The obtained formula for the magni-_ 
th of the basin to 


ws. 


_ Note,—Discussion open until August 1, 1961, To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE, This paper is part | 
of the copyrighted Journal of the Hydraulics Div “ery Proceedings of the American 1 So- : 
ciety of Civil Engineers, Vol. 87, No. HY 2, March, 1961 
I Ingenieur-Docteur, Special Lecturer in G1 Graduate Studies and Research Associite : 


Tes 


at Queen’ s Univ Kingston, Ontario, Canada, 


March, 1961 HY 2 
: 
1% 
@ 
vat 
| 
= 


the theory of in a harbor could be considered as a natural 
come of the experimental results presented? | in 1954 and3 1956 and of the - 


author’ s theory of periodical gravity waves on a — 


A harbor is considered herein as a combination Lof separated 
by distances that define the lengths of the basins. The terminal discontinuity 
is, of course, a total obstruction. This total obstruction, frequently formed by 
beach, dissipates a noticeable part of the incident wave energy. 
ie _ With this assumption, the value of the agitation in a harbor may be com- 
= from a general theory, that considers the behavior of waves acting on a 


combination of a number of discontinuities. This broad theory may be applied 


=> 
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wave length is important because of multiple reflections. Resonance 
___ Results obtained by theory are compared with those obtained experimentally. § [7 
| 
q 
7 = = 
&§ 
— 
f 
— 
“Mouvements de Resonance a Deux Dimensions dans une Enceinte sous l’action 
d’Ondes Incidentes,” by F. Biesel and B. Le Mehaute, La Houille Blanche, July, August, — 
Le Méhauté, Journal, ASCE, Vol. 86, — 


to a 1 number of such phenomena as waves passing over a shelf, agitation be- = 


an However, _ because of its importance in engineering practice, and because — 


wnat as shown by Fig. 1, is considered herein. It is, of course, a ema «il 

_ case in the more general theory of waves on a combination of discontinuities, 7 

in the” _author’s previous” work 4 great u use is made of the complex. 
; ig It will be noticed that this method of complex | number ete provides a =i 
simple and ‘powerful tool for the analysis of gravity wave motion. Solutions 


will be obtained by this method with a minimum of computation for cases 
“% which would usually be solved by finding a complicated potential function, go. 
The method of the potential function, 9, is always tedious if there are disconti- 
= inthe boundary conditions, 
An example is presented which is consistent with the theory 

‘The theory - is limited in that only termsof the first order of approximation 
are considered. That is, it deals only with the linear solutions in which, for 
example, the convective inertia effect, ‘that is proportional to the square ia 
the velocity, neglected. This ‘particular assumption seems valid for ‘the 
majority of cases. Also, the theory deals only with the main vertical and 7 
} horizontal components of ‘sinusoidal motion. This stipulation tends to make — 
the theory invalid as the ‘discontinuity is: approached, but it seems sensibly 
valid at distances of the order of two or three times the depth from the dis-_ 


| continuity. Three dimensional components, if existing, must have a limited 


effect. "They are theoretically negligible when the width of the basin is less” 
_ than half the length of the wave inthe same domain. However, even for greater — 
width, experiments indicate that this theory provides a fairly good approxi- 
: mation for the value of agitation, as long as no transverse resonance occurs. J 


Transverse resonance is likely to. occur when the width of the basin is very 
close to a multiple of half wave lengths.4 ee en 
ee A linear periodical gravity wave may be definedby a rotating vector whose 
modulus is equal to half the wave height and whose argument is equal to the 
; 4 .- of the free surface at agiventime and at a given point. In the same z= 


4 


a clapotis is defined by two vectors rotating in Opposite directions, meee 
tation.—The “letter symbols adopted for use in this paper defined 
_@ where they first appear, _ in the illustrations or in the text, and are arranged | oe | 


alphabetically, for convenience of reference, in the Appendix, — orl 


_ COMPLEX VALUE OF THE AGITATION IN A BASIN 
4 : oss Consider a | basin n aligned in in the ‘direction of the incident wave as shown by  &g 

: = The incident wave, coming from the open sea arrives | at the mouth of this | 


basin with a half amplitude Aj and a phase VE at a given time. Ba Pulte oe 


This wave is considered to — by the complex number Ay modu- 


Jus — Ay. . Similarly, the wave reflected by the over. 
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; all and basi) is is detined by a complex m number By of 


wave Ai is partially transmitted and partially 
transmitted wave is ay Ay of half amplitude a phase(4, 
4 
ra] 


+ Ay). ‘The terms and hi to oft 


reflection, respectively. They are defined by complex numbers of moduli 


and and arguments aj and a, respectively. 
These moduli (a a and py): are equal to the ratios of the transmitted and re- 
“fected wave heights to the the incident wave height, it, respectively. — Similarly, | 


_ FIG. 2,—NOTATION OF PRIMARY WAVES 
and BY are equal to the c change e of phases, caused by the 
_ discontinuity, of the transmitted and reflected waves, respectively 2 
Now, consider a wave defined by a complex number Dj arriving from the 
ov at the discontinuity. _A part 89 Dy is reflected within the basin, and © 
another part a2 Dj is transmitted to the open sea. The terms @» and Bo are 
two complex numbers defining the transmission coefficient towards t! the open 
sea and the reflection coefficient towards the basin. 
The four complex mumbers or eight parameters moduli 
a9, B 


and Ci. Dy inside the basin at he entrance to the harbor. - 
waves have the same period pe 
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AG AGITATION 


wave moving inside the basin a away the entrance, arrives at 


end of the basin, with a different phase r. This phase is a function of the rela- 
of the basin with to inthe basin 


inwhich 1 is the length of the basin an aad Lat the length of the wave in the 


7 tween the 1 wave » at the entrance. and the s same wave arriving ng at the back. ‘The 
_ wave height remains constant if friction is neglected. However, if friction is 
: _ included, as produced f for example by a wave trap or by lateral beaches, then 
the amplitude of the wave is reduced in the ratio r(z <1). Asa first on 


in which K is a friction coefficient. algae 6: 
mt Hence, the incident wave at the back of the basin intel 


A 
in which r is a complex number of modulus r, and argument Fr a. 


_ In the case of a single basin, the obstruction i is complete, and the wave Co 


eas 


me 


which p is a coefficient of reflection ieee 4) that depends on the steepness of 7 


rs beach. The ‘coefficient . Pp may be computed by. use of theory 5 or be obtained 
_ experimentally as in the chart given by Messrs. Greslov and Mahe.® Y As a 


a) “Etude du coefficient de reflexion d’une houle sur un obstacle constitue par un — 
incline,” and Mahé , Proceedings, ¢ Coastal 5th Conference, 
ty) “Le pouvoir “nie des ouvrages maritimes exposes a v action de la houle,” 


M. Miche, Annales des Ponts et t Chaussées, Mai, Juin, 1951, 


4 
| 
“4 
— 
| 
4 
7 _ _ However, the loss of energy, such as, that due to the breaking | &§ 
ar = on a terminal beach, must often be considered. Hen¢ j 
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first approximation, p may be constdered asa a function m of the ratio of th ~?_ 
length of the beach to the width of the basin, assuming that all the wave — 
arriving at a beach is dissipated in the breaking zone. — ena pape ate 
_ However for a long wave, such as a seiche type wave for which the iia a 
is of the order of 100 sec the reflection on a =" is ermal total even when 
it has a very gentle slope : so t that we have Eq. 


1, 2 2, 5, Qa, and 10 of tive ‘that 


cacao of the five unknowns Cy, © Do Dy and By as functions of the 


incident wa waves and the coefficients: 7 ij/is each a complex 


inciden 
number of modulus 1Tl = Ls and eee 0, that is the real number $. 


7 (It should be noted that i is, in this case, hcaiaaie from the ‘puny num- 


- following formulas are obtained from the previous equations: 


equation —_ 


at the bac back bythe the equation 


Do) = os 


value in the sea, at the mouth, is 


— 
f the harbor becomes Dj, such 
7 
G 
— 
3 


d difference ¢ of phase Ty + By) + + Dy), ‘that ‘is 


———, gives the difference of phase of the agitation between the hacia 


of the entrance and enables an estimate to be made of the current at the « en 
_ Verification.—A simple verification of these formulas consists in comput= 
‘ing the absolute value of By when the loss of energy in the overall layout is_ 
zero, that is when p= 1 and r = 1. In this case, By should be equal to Ay. j 


ao 
Bg 


Moreover ‘the ‘value of the in at the back of the harbor 
, when it is assumed that all the incident ee is a bya terminal 
beach such ait is found to be 
2 +D, | =2 = 2 r 
At the entrance hab. = 


from the open. sea into the basin. 


ABSOLUTE | VALUE OF THE “AGITATION 
, The value of the agitation at the rear of ‘the harbor is given by th the absolute 
value 2 | Co + Do |; at the entrance in the basin by the absolute value 2 | Cy 
+D! ; and at the entrance ! inthe open sea by the absolute value 2! A Ai + By # 
The amplification factor is given by dividing these values by the amplitude i 


= 


the incident wave, Ay 


the back o f the harbor, ‘the factor is found from Eq. 17 to be 


[1 +( + (p 2 -2p cos | (2 4, 


q 
ag 
ag Introducing Eqs. 19 and 20 and utilizing t iously pre- 
: 
a represent tne transmitted Wave 
4 
| 


‘sinh sinh 2 my hy m, 4 
‘Hence, it is possible to express the value of the agitation as a function of Z, 
"The term a, may be expressed as a function of the shape of the disconti- a 

by the empirical formula he; ‘Al 


in which b is the width of the entrance of the basin. 


"PRACTICAL FORMULAS 


ravelling 
from the + pavetranes to the | back of the basin; p is the coefficient of reflexion of 7 
ix: the back of the harbor (p is small when this back is a beach); ay is the coef- 7 
.. ficient of transmission for the wave arriving from the open sea into the basin, _ 
| ~ and is computed from Eq. 29; fg is the coefficient of reflexion of the entrance _ 
; | for the wave travelling inside the basin in the direction of the open sea. - Be is ae 
_ given as a function of a1 by Eq. 24; y is the difference in phase in the a 
Between a wave at the entrance and a wave at the back of the harbor. oll 


terms lg and 1, define the widt — 
mm the basin, respectively and Lg and Ly are the wave lengths in hese two do- aa ie 
term L is a function of the period Tandthe depthh: 
A 
q 
computation of the val- 
n .-’ ue of the agitation at the rear of the basin for a given incident wave of ampli- — f : 
Ly 
The value of the 1 is given by the theoreti- 


—— - 


WAVE AGITATI ON 


given by Eq. 3; 3; My is th the difference in phase — wave reflected by 
the entrance inside the basin inthe direction of the basin and the wave travel- 
ling in the direction of the open sea also inside the basin; > is zero in the 
case of a total obstruction and equal to 7 in case of a total opening in a wide 
and deep open sea; and Br is given by Eq. 25. ie ae 
In these formulas, use is made of coefficients: — the shoaling coefficient 
= 2 as defined by Eq. 28 and the change of domain coefficient Z 
ponent by any change from the open sea to the basin, and defined by Eq. 26. 

_ Example.--Consider a basin as shown by Fig. 2 in which the length of the 
basin 1 = 1450 ft; the width of the basin 1g = 500 ft; the length of the mouth b — 
= 187 ft; and the depth of the basin hg = = 32 ft. A terminal beach exists witha — 
very gentle slope and a length of 340 ft. _ The open sea (or another basin) is ; 
66 ft deep (hy) and 625 ft wide (ly). The incident wave has a period of 10sec i 
2 a height of 2 Ay = 15 ft. 
: _ Preliminary Computation. —A 10 sec period wave in the basin in a depth — 
* = 32 ft, has a wave-length Lg = 300 ft. A 10 sec period wave in bad open 

1 = 66 ft has a length: Ly = 400 ft. | Hence: 


elg=. = 375 ft- >= and and ly = 750 ft > the the theory is “not strictly ap 


"plicable should a fairly good to the of the agi- 
tation. A mieten resonance is unlikely because lp and ly are not close to 


any multiple « of 


computation of the shoaling fi by means of Eq. as fol- 


7) (32), 


sinh hg: = 1.7 


2 my =2 ~ hy 
400 


Substituting 


‘into. Eq. 26 yield th the ¢ ame of domain ¢ coefficient 


= 0.75 x 0.80 = 


— 
— 
4 
= 
— 


= (1. = 1. 06 

pa 

©. 


= 1, 1.06 x0. 0.61 x1. 1.08 = 0.700 

From Eq. 24 the coefficient is found 


2 x 0. 


the phase e of the entrance is found to al 


_ From Eq. 3 the phase 1 lag o of the basin is found to be _ 


€08(2'F r+ B2) = = cos (120° + 45°) = =- - 0. 


a4 


a. The terminal beach reflection coefficient p is then ‘computed. | When the 
beach length A is small relative to the width of the basin lo, the coefficient of 
reflection of the back of this basin can be computed by assuming that that the w wave 
energy | dissipated is proportional to the beach length sothat; 
a 


sen A is of the same order as lg, it is more accurate to compute the reflec 


It is assumed that ‘the friction inside | the is by putting r 


tion coefficient as a function of the beach length, its slope and the steepness 
In this example, it is 


of the incident wave — this example, it is assumed that ee 0.1. - 


the value of the amplification fa be the us use e of 


the Eq. 23; © 


A 10 x 0,00 x 0.82)” - 2(0,10 x 0,00 x 0,82) (0.97) 

of the amplitude of the agitation at the back of the harbor is 


= 0.622 Ay = 0.6215 = 9 ft 


and 
and substituting into Eq. 29 yields the transmission coefficient ‘ 
q 
36 
a 
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WAVE AGITATION 
4 tt is. likely that this value is ‘slightly over- estimated. Because lg > 5 


and yy >= it is is that some secondary three dimensional motions 


ona superimposed on the primary two dimensional ‘motion. 
‘The Importance of the | Ratio of the Length of the Basin to the Wave Length.—_ 
It be that a in the term r) the 


2 p pr2| Ro cos cos (2% r+ + By) us usually h has a a large effect. on the > value of the agitation. 7 
ae. If 1 is of the same order as Lg, a slight variation n of T, does not change the 
value of the « cosine given by the previous formula nor, consequently, the agi- 
tation very much. In practice, this is the case for long waves such as tides” 
‘in a bay or seiches ina basin. But as far as ordinary wind generated waves 
and swell are concerned, lis. generally much greater than Lo, and, a 
—_ change of the period of the incident wave T, causes a ee ¢ in the voles ? 


of x that introduces a great variation in the v value of cos (27 r+ Bp) an and ‘conse- 


Because waves in the field are always ‘irregular it is necessary to consider — 
“degree of agitation” that could be found by the mean value’ of 
- for f» ), or or ‘more simply of 
the wave length, , that is the wave period, varies. 
A it is interesting to note that the “degree of agitation” does not depend on ae 
whereas the value of the agitation for a given period is very much 


RESONANCE 
“4 Maximum agitation occurs when: 


nay 
gives 
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the: 


are all 


The corresponding a agitation is to 
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verified by a 


perimental work have: been presented previously. For case 


‘ot the relative length of the basin for various shapes of the entrance and vari- 


loss of energy in the is 


No damping of the wave: re = 
These e theoretical curves are compared with the experimental results in _ 


4, 6, and 

if ‘The maximum value of the agitation | that occurs, of course, at siege 

nance has been analyzed and presented in Fig. 5 for the case of a basin closed > 
“by a simple obstruction. _ The computed values presented a: assume 


and b) a constant value of the friction terms p and r. 


‘These two theoretical curves are compared with the experiment. 


Successive have been made for the following entrances: 


1. Obstruction i in ; a flume. e. (The domains 1 and 2 have the same charac- 


teristic. ) (Figs. 3, 4 and 5) 
2, Change of width, constant depth, and full opening. (Fig. 6) 
Change of depth, constant width, and full opening. (Fig. 7) 
_ has also been verified that this theory could a applied to a number of 


_ other combinations with the same degree of accuracy. 


i 
— 
The minimum amplitude of agitation is equal to 
4 The theory presented herein h number of experiments 7 —— 
= 
le 
a is they are given as functions of the opening of the basin and the depths and a oe : 
a : yes are computed assuming — 
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6. AGITATION IN A BASIN WITH A CHANGE OF WIDTH 
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T= 2sec. 


Period : 
Incident Wave Height: 2a=Iicm) 


03 


Amplitude of Agitation:Ay| = 
Incident¥pWave Height :A,| Period: T= 2sec, 
Incident Wave Height: 20=!cm 


Curve 


mod. 
| Length of the Basin (f) 
| —~Wave Length (Lo) 
eng 
05 


— 
> 
4 
FIG, 7.—AGITATION IN A BASIN WITH A CHANGE OF DEPTH 


WAVE AGITATION 
The ‘experimental and theoretical data presented in Figs. 3, 4 , 5, and 6 are 
considered to exhibit fairly good agreement. ad 
— In the” case of a basin closed by an obstruction, the friction has, as would 
seem logical, an important effect particularly on the maximum value of the 
- agitation that occurs at resonance. Here it is necessary to take into account — 


‘determine accurately. theory, particularly concerning. the values of 


‘causing the resonance, is well verified, 
For the case of a simple change of width, the friction effect may be neg- ’ 
lected, (p = 1, r = 1) except when a terminal beach exists: p<1. The com- | 

y = values of the maximum (resonance) agitation is particularly well veri- 4 


fied by the experiments, but there isa slight discrepancy between the theo. 


retical values of ‘Ba causing the resonance and those obtained by experiment. 


The means that the value of cos fg by introducing ay Ag) 


is incorrect as far as the v value of Bo is concerned. AS previously seen, the 

“degree of f agitation caused by an irregular \ wave does not depend on fg. Hence, 7 


aff computation of the e degree of agitation caused by a an irregular wave | is not 


affected by this discrepancy, 
In the case ofa simple change of depth, the o occurrence of resonance asa 7 


function of 1/Lo is well verified. _ However, the magnitude of this maximum © 
seems —* than that computed from the theory. 


verified by experiments and | can be successfully applied to ‘compute the period 


and the amplitude of resonance, particularly if friction effects are small. 
Above all, the d degree of selectivity of resonance given by the theory i is 3 always — 
accordance with the experimental results. 


ickness of the peaks of the agitation curves. 
‘The is found both by theory anc and experiment: 


| 


4. 1. That the o occurrence * resonance increases with a a dene of depth and 
a change of width, and decreases with an obstruction such that it tends to zero — 


That the value of the agitation at resonance increases as Lj 
or r hy /hg and l9/b increase. It will be noticed that the value of agitation tends _ 
: to infinity with lo, /o, that is in the case of total obstruction, when the friction — 
is neglected, but at the same time the occurrence of resonance tends to zero. 
- However, when friction effects are taken into account, the value a the aeme- 
_ tion tends to zero when lp/b tends to infinity. 
3. ‘That resonance occurs when 
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in the case of a complete at 


3) 
in n the case of a change > of depth: without ot obstruction; 
4. That the minimum m agitation occurs when: when: 


ai - ‘That the agitation is proportional to the amplitude of the incident wave, _ 
_ even inthe case of viscous friction as long as the quadratic friction is consid- 


It is relatively easy to tiie approximately | the reflection coefficient of a a 
beach in nature, that is, it is easy to find ead value of p. It is certainly more = — 

. difficult to determine the damping of a wave. . Damping is particularly impor- : 
tant in laboratory ‘scale tests in which the motion in the boundary layer is 

- likely to be laminar. In nature, as a first approximation, this effect could - 
_ probably be neglected, that is r = 1, except when a wave trap or the presence — 


of lateral beaches could produce important frictional effects. _ iJ 


CONCLUSION 
‘The motion of a wave on 1 a discontinuity has been analyzed p previously.4 . 
In this: paper, _ considering the entrance and the back of a harbor as disconti- q 
nuities, theory, has bees: formulated that permits the computation of the 
value of the agitation in a basin of relatively simple shape. Such an idealized _ 
basin is rarely met in practice, but a number of cases may be considered as © 7 
; very” similar. This theory is valid to explain not only some phenomena of 
agitation in harbors but seiche motion and even tidal motion in bays and on 
‘the continental shelf, 
main interest of the theory is to show quantitatively in an in idealized 
_ case and qualitatively in some more complex cases the importance of a num- — 
ber of | parameters such as the shape ofthe entrance, the slope of the oor ence 
: beach, and also the importance of multiple reflections and the number of \ wave 
- lengths contained in a basin. It has been shown that for the case of resonance, 


1 the agitation in a basin could become greater inside the harbor than in the | 


Only the main topics have been emphasized in this paper, but it 1 is possible aa 
by the set of formulas presented herein to study the previously unknown in- _ 


of a very great number of other parameters. 


This” theory is limited to the computation of the agitation in a harbor sub- ; 
jected to a wave of constant amplitude and constant period. In the field, waves ‘i 


a 
mm 
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in which the probability of ‘obtaining a degree of of agitation is is to be computed 


asa function of the incident wave spectrum. ies gh | 
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ENDIX. 


reference and for use by discussers: 
Aq = (semi-amplitude of and argument, 
= (phase of the incident wave at the entrance to the harbor). 


= at ‘the back of the basin. basin, 
= agitation with ‘resonance. 
minimum agitation, opposite resonance. 
min 
= (semi- -amplitude of the reflected wave seawards) and argument, 
= of the reflected wave at the entrance to the 


b of t the entrance. 


= number of modulus. 


= shoaling (see Eq. 28). 


complex number of modulus: Co and argument: Co: half the the one: 
height and phase of the incident wave at the back of the basin, re an 


complex number of modulus. a 


= (semi- -amplitude of the wave arriving from the basin tothe entrance) 
and argument-phase. 


= at the entrance. 


= complex number ae modulus: : Do and argument Do: half the wave 


hetight and phase a the he reflected wave at the back of the basin, re- a 


Al 
— 

here for 
— 

= 

f 
g towards the basin at the en- 

e 1 
Cc 

v4 | 

> ac - 

— 
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1 = depth of the open sea. 
depth of the basi 
complex number of modulus: i=41, argument: : 


coefficient of friction, 
L; = wave length at the open sea and: ™-* 


Le 


‘= wave length in the bast operas 
wave len n the basin an m9 


width of the open sea. 
1g width of the basin. 


=: = length of the basin, 


_ coefficient of reflection of the back of the t basin. _ «2 1 


= of the wave- -height in the basin r= =e 


= complex number of modulus: and argument: Te 


= difference of phase betwee between the the waves at the t two extremities a the 


see . 26), 


ay number of modulus: and : @ 
= coefficient of transmission of the wave Ay: ‘ratio of the transmitted 4 


wave height to the incident wave height: Ay. 


=... complex number of modulus: ao 


op = ‘coefficient of transmission of the wave Dy: _ ratio of the transmitted 
wave height to the incident wave-height D4. 
= = change of phase of the transmitted wave at ~ entrance - seawards. 


complex number of modulus and argument 4. 


os coefficient of reflection of the wave: Ay : ratio of th the reflected wave 


os change of | phase of the reflected wave height at the « entrance sea- 


be complex nur number of modulus: Bo 5 and ‘argument: 


= = coefficient of reflection of the wave: Dy: ratio of the reflected bod wave 
towards the basin to the wave height: Dy. 


= @eeee @ of “of phase of the e reflected wave towards the basin at the en- 


trance. 


back of the rr 


seawards to the incident wave Ay 


f: 
rey. 
a 
“| 
4 
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Proceedings of the American Society of Civil Engineers 


wae The se sediment deposits most characteristic of estuaries 2 are muds which con-_ 


: tain a large percentage of clays. These clays givethe mud its particular qual- 
ities. The transport of these muds is largely determined by flocculation char- 


alyzed. The process of deposition and scour of the bed of these sediments is 
also described. All of these processes are finally combined in a description 
of the process of shoaling in Mare Island an deepened side 


acteristic of suspended clays. The flocculation process is described and an- J 


Sediment a or bay system to an ‘inter- 


ested observer to have several curious aspects. Sediment materials, trans- 


ported by land drainage, enter an estuary in varying amounts with the inflowing — 
waters. Except during exceptionally high river inflows, little sediment may 
leave the estuary tothe ocean. ‘It is apparent that river-borne sediment is re- | 


- tained by the estuary. During periods of low sediment inflow to_ an pene 


however, shoaling inchannels and harborareas may occurat much larger rates © 


bs Note.—Discussion open until August 1, 1961. To extend the closing date one month, 

a written request must be filed with the Executive Secretary, ASCE, This paper is part _ 

of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- — 
ciety of Civil Engineers, Vol. 87, No. HY 2, March, 1961, | Vi 
1 prof., Hydr. . Engrg., Dept. of Engrg., Div. of Hydr, and San San, Engrg., ‘Univ. of Calif. 

Asst, Research Engr., of Div. of and San, Engrg., ‘Univ, 


a 
— 

ESTUARIAL SEDIMENT TRANSPORT PATTERNS 
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_ than are at s of hi, high th inflow : rate. The amount t of re- 
moved for maintenance of navigable channels often exceeds the total sediment - 
inflow to the estuary. That the estuarial sediment movement consists of more 
than the deposition of the currently river-borne sediments is shown also by 7 
varying degrees | and locations of turbidity, often visible at the water ‘surface, 7 ‘ 
that occur at times unrelated to the varying river inflow. 

~The estuarial | environment obscures even the simplest information on wend 

“ment movement. — . Tidal fluctuations at t at the mouth of an estuary or bay are pro- 

_ Pagated throughout the estuary with 1 velocities and amplitudes that depend on 

a local bottom and shoreline configurations. These tidal changes, often modified 
resonance together with river inflow, result in roughly periodic, , continuously 


* changing circulation patterns withinthe estuary. Water current velocities ee 


‘from zero at slack water to high velocities, particularly at narrows, during the 
Salinity changes are also characteristic of an estuarial environment. During 
_ “high runnoff periods, fresh river waters may extend to the ocean. During nor- — 
mal or low flows, however, estuarial waters are commonly saline, with a 
| . ing 1 rivers. The > salinity often varies with tidal currents, particularly 1 near the a 
- river mouth. Regular wind patterns are often characteristic of the coastal a- _ 
reas (of which estuaries ar area part. 1 In addition t to’ modifications of) water — 


peste alter the character of the sediment transporting medium in sallow 
water. Extensive shallows formed by accumulated sediment and cut by Sparse — 


dendritic drainage channels are characteristic of bays that might by their na- 
ture be expected to have an important role in the estuarial sediment regime. | : 


_ The San Francisco Bay system in California, for example, has all of the 
characteristics mentioned. Drainage from the entire Central Valley area 
discharged at the northeastern end of the bay, and numerous small streams 

: _ draining the neighboring area enter the bay at widely distributed = 
bay system covers an area of 485 sq miles, and 70% of this area is less than 
15 ft deep at mean lower low water. Large areas less than 6fft deep occur. 
‘The bay is about 55 miles long, extending both north and south of a narrow 
mouth. Resonant tidal variations, excited predominantly by diurnal tide, con-_ 
_ tribute to water interchange between the north and south bays, and river dis- 


charges modify the northern current patterns. 


a this bay should also be , applicable to them. Recent studies of the sediment be- 
havior in the portions of the bay near the entrance of the predominant land 
drainage have revealed the sediment transport processes that comprise the 
4 sediment regime. . This paper presents a brief description of the factors that , 
largely determine this regime and the oe of nn factors in over- = 
‘THE 
of sediment from the bottom of potate distributed around the bay 
system indicate that, except at the mouth of the pe See channels, — 


"particle : size ranging up to about 100 1 mv. . The of the fines" 
ogee — watershed includes many clay minerals of ae most are — 
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sified a as illite and smaller amounts: of the kaolinite group 
added. The bay sediment shows the same composition. = | | 
These minerals respond rapidly to changes | in their hydraulic and chemical 


their crystal structures. The minerals are commonly thin plates, 
i: ‘Sheets, or tubes, and range | in effective diameters from a few microns to small ; 
"fractions of a micron. The arrangements of atoms in these minerals are such > 
- 2 that the large surfaces of the minerals have a negative charge, and when _ 
7 rounded by water containing salts, cations accumulate near the surfaces ina © . 


cloud. In waters containing small amounts of salt, these clouds are large and 
interact electrically, preventing the particles from coming close together, Such 
“a suspensions are called “dispersed” and are typical of suspensions in many — 
streams. Since small particles settle only slowly, | dispersed particles tend to 
remain permanently in suspension in flowing water, and represent there part | 
When these particles are suspended in water containing a gram of sea a salt 
per liter, or more (sea water contains approximately 34 g per l),t ¥ 


the cation cloud is repressed by the abundance of anions in solution, and par- 
_ ticles do not repel one another at the larger distances. They can approach one 4 
another freely, and an attractive force, that operates. only at short range, 
_ _ predominate and mutually attach the particles. _ They then stick together, and 
if successive particle collisions occur, a multi-particle aggregation, or floc, | 
a form. The shapes of clay particles cause the structure of a floc to be open 
and largely to contain water as does a sponge. The reduced density and the in- 


creased size change | the e hydraulic character of suspended sediment. nh 

The two requirements for the formation of flocs from | dispersed suspensions 
gy are the reduction of the repulsive forces and the movement of one particle to 


another. first requirement, the salinity at least 1 g per occurs 

most of an estuarial system to the extent that nearly every collision results in 7 

adherence. Collision of particles can occur in several ways: (1) by Brownian — 

motion, or random thermal movements of small suspended particles; (2) b . 

_ internal shear motion of the water; and (3) by differential particle settling ven 

locities. Since flocculation determines the hydraulic transport character of 

. 7 suspended clayey sediment in an estuarial environment, it is useful to consider i 
the relative importance of these different ways of bringing particles together. » | 
_ Atheory of rapid flocculation was first given by Smoluckowski in 1917-18. | 

‘This theory, and extensions by later workers, is given by J. Th. G. Overbeek. 4 


The probability of collision from Brownian a Le can be approximated 7 


83 “Kinetics of Flocculation,” by J. Th. G. Overbeek, Colloid Beience L Scien 


lishing Co., , 1952, PP. 278-3-1 -3-1, 
af 
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wh which k Biltzman’s constant, 1.38 x 10 ergs: T denotes es ab- 
andy temperature, in K; vo is the number of primary particles per cu cm; 


and 7 represents the viscosity, in poises. Suspended sediment concentrations _ 
* _ in channels of San Francisco Bay commonly range between a few and 100 parts 
per million. Taking 0.3 mv as a representative spherical diameter and 15 C | 
a typical temperature, , 1 = 0.038 per sec. As primary particles join ¢ other 
primary particles and flocs, the total number of multiplets is reduced, increas- ~~ 
_ ing the distance between them and reducing the frequency of collisions. _ row 
_ Each time a collision occurs to all particles, the number of particles is sre- 
, duced to one half of the | previous 1s number, but the particle weight is doubled. 
The rate per unit time at which the particles grow in mass during any part of 
the growth period is proportional to their size, but inversely i 
- the time between collisions, whichis itself inversely proportional tothe number | 


— 


of particles or proportional tothe weight of the particles. This means that the ~ 
rate of particle growth is constant, or for the prior given 0.3-mv, original size 
particles, the growth for any size floc in 26 sec is equal tothe mass of an orig- -— 


e = _ The time to build particles of 100-mv average size andof a specific ray 
of 1.10 from particles of 0.3-mv diameter and 2. 65 specific weight is is - 


(265-1) onc = = 670 days” 
while the growth from an ; average s size of 0.3-mv to a floc size of 10 mv takes 
only 16 hr. This process, which does not require any outside driving force, — 
must be expected tocreatethe smaller floc in a short period, whereinthelarge _ 
flocs will probably appear only if one of the other effects enters the picture. _ 
Only a much faster build-up of flocs could actually maintain an average size 
much larger than 10° mv if the natural breakage of these rather fragile struc- _ 
bares is taken into consideration. 


flowing water cé can be calculated by finding the probability o of one particle pass- 
ing through a 1 sphere of radius R within which a lasting bond between particles : 


is formed. The simplest spherical radius would be the sum of the radii of the 
colliding particles. Fora water shearing , du/dz, and a concentration of 

the probability of collision Jis 


motion yields 


A comparison with the probability of of collision from 


the 


The relative importance of internal shear strongly grows with the — or par- 
ticle size. For a shearing rate of 1 = sec anda particle. or floc diameter “4 : 


age 
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except n near the boundaries and in n strong turbulence. % The bulk of sediment al 


4 pended in the flow must form flocs of moderate size, therefore, without sig- 
nificant increase in flocculation rate by shear. Under wave conditions, unstable 
& flow around pilings or or - other « objects, , the flocculation rate, », enhanced by 
_ local turbulence, can “snowball” and produce flocs of sizes that rapidly settle. 
a When particles having a wide range of sizes are suspended together, larger 
ones settle faster than smaller ones, and increased frequency of collisions oc- 
cur. H. ‘Miller4 studied this mechanism and showed that a significant effect 

is found when the radius, r, of the particle 


It be 


or flow-disturbing works. | 


_ While the individual clay particles a are so small that Brownian / motion alone 
is almost able to keep them in permanent suspension, the flocs of such particles’ 
are sufficiently great and heavy to make them settle out in still water. Large : 
flocs | have a settling velocity measured in feet per hour. Even the slack water » 
‘periods in tidal areas are sufficiently long to permit considerable amounts of | 
-flocculated clay toaccumulate onthe bed. 
© _ Laboratory e experiments show that at a density of about 10 g per 1, the a. aol 
flocs form a continuous structure. density is obtained by settling of in- 
dividual flocs. Settling near this density occurs at a much reduced velocity 
compared to that at low concentrations, because of the increased resistance of 
the waters escaping from between the many flocs. One speaks here about “ hind-— 
ered settling.” ” At a density of 10 g per 1, the flocs make bodily contact and 1 
form a continuous cohesive mass of the same density and strength as the in- 
dividual large flocs. This mass can support a measurable shear stress with- 
out permanent deformation, at least for restricted duration. Upon prolonged — 
stress, it will very slowly creep. Such creep may be either a shear motion or — 
a compaction, or both, depending on the type of stress applied. Compaction 
‘particularly important because it causes the strength of the deposit to increase 


ast more a and ‘more bonds are mate. an 


Muller, Kolloidchem, Beihefte, Vol. No, 223,1928, 


which g 1s acceleration Of gravity and p 1s the density of the particles. e 
— term in brackets corresponds to a 5.0-mv diameter solid particle ora 6.3-mv_ ff 
— 
 && which corresponds to a 2.1-mv solid particle, or toa 2.6-mv floc. 7 7 

that this mechanism is most important for collection of small flocs 
a larger ones. This mechanism intensifies the formation of large flocs after ac- i 
and further enhances deposition near pilings 

— 
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| 
arn occurs not only in still water, but also in flowing water. Up toa a 
flow velocity of 0.5 fps, deposition is the same as in still water, that is, the FF a 
_ entire flocs are deposited at and bonded to the bed. _ As the » velocity increases, 
some of the loosely attached particles are occasionally torn off the bottom while = 


= tougher layer of clay develops at the bed surface, but at a much slower 


smaller flocs or parts of larger flocs can still deposit on the bed. A slower 
3 Fate. Sse ey bonds between clay minerals ca can be visualized to develop a 


qe structure. A wide range ze of ‘such structures. occurs in n nature, each being 

characteristic for a particular velocity of flow during deposition. This picture 

is further by ‘secondary compaction of | of the deposits: due to gravity. 


The scour of nein deposits has s not been s studied with as much detail 
as has their deposition. It is unknown whether scour is caused predominantly — 7 


by the shear stress at the bed or by pressure fluctuations. It has been ob- 


served, however, that scour does not occur by chipping large complexes from ; 
= bed, but rather by the continuous wearing off of small particles. The speed a 


of scour is reduced with increasing concentration of the flowing water, indi- a 
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BOTTOM OF MARE ‘ISLAND ST STRAIT 


AT BOTTOM OF MARE ISLAND STRAIT 


osition. attempt been made, tar, to define the condititine of such an i 


In returning: to the problems of the previosuly quoted San Francisco Bay, 
some particular observation m may be described that resulted from a | prototype — 
_ study which had the purpose of providing a basis for the verification of the © 
_ Corps of Engineers’ hydraulic model of the Bay at Sausalito, Calif. Oné of the 
most important problems is the very high rate of mud deposit in Mare Island = ‘ 
7 ‘Strait, the lower endof the Napa River channel, just above its junction with the 
main Bay system. This junction is at Carquinez Straits, a narrows of the hn 
in which the fast tidal currents prevent deposition of mud. The deposited sed- 
iment was identified as identical with that of the bay | proper, . not with that of © 
_ ‘The most perplexing observation was made in connection with the time “ 
which the sediment was deposited. It was found that deposition occurred pre- 
dominantly during the summer months when the bay water contains the lowest | 
_ sediment concentration and when the rivers carry practically clear water. “Tt 
thus, necessary to find immediate source of somewhere 
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"Island Strait,a number of specific observations were made. "First, a canna 
charge of bay sediment was released at the beginning of flood tide at the _ 
_of Mare Island Strait, after it was identified by radioactive gold label. This” 
trace sediment was found to have been almost fully deposited inthe straits a 
ing the first tidal cycle. It was found, furthermore, not to have moved signifi-_ 
cantly during the next few days during which the activity was still sufficient to 
; be identified. Within the Strait, distribution of the radioactivity was the same _ 
as that of the normally deposited sediment (Fig. 1). Characteristically, the a, 
highest rates of deposition occur on the two sides behind piles andother struc- 
_ tures which provide local turbulence that causes fast growth of large flocs. —” 
This is shown in Fig. 1 by the 500- count - min contour along the west shore 


2 shoaling material enters in suspension reel the mouth of Mare Island Strait, 
which enters the area from Carquinez Ll J 
- Next, it must be asked, from where may suspended s oulianens originate dur- 
ing the summer months. Immediately west (this is, downstream) from the 
- junction of Mare Island Strait, the main channelof Carquinez Strait widens into 
San Pablo Bay. This large body of water has wide, shallow, tidal flats | in which | 
the regular summer winds generate strong waves. These break near the en- 
7 : trance to Carquinez Strait and suspend large amounts of sediment. The tur- - 
_ bidity at the surface, resulting from this action, can be “observed 
we _ Finally, another test with radioactive tracer was made to study the removal — 
of sediment from the tidal flats during periods of wave action. Bay ———— 
. labeled with radioactive gold, was placed in suspension above the bottom during | 
g a quiet period. _ The tracer was placed over two shallow areas of ‘San Pablo 
‘Bay, and part of the. tracer was deposited with the natural sediment. These 
Later, a moderate 
dition ‘occurred and it was found to have diaturbed the shallow deposit some- 
_ What, but to have left the deeper deposit undisturbed. The disturbance appeared , 
the form of _elongations of the pattern of 


in the but does not move ita significant distance. Transport 


of from mud flats is predominantly by in moving water. 


>, 


 CONCL LUSIONS — 


_— studies of the motion of iiiniaian is sediment in estuaries, , both in 
_ the laboratory a: and in the San n Francisco Bay system, it was found that the clay 
sediments flocculate slowly as they advance from the rivers into water with 
y ‘Salinity of 1 g per 1 and more. A faster development of the largest flocs which 
settle rapidly is caused by a general turbulence which exists in the the exten 
tidal flats. This sediment supply occurs during winter floods. 
_ During almost daily summer winds, some of these winter deposits a1 are e stirred . 
_ upby wave action, and high concentrations of flocculated particles occur. . These 
moved by local currents” into the quiet areas of the artifically | deepened 
channels, particularly where pilings and other harbor structures cause rapid > 


| See of settleable flocs. Here the sediment deposits andnecessitates ex- 
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PROBLEMS IN HYDRAULIC STRUCTURES. 
By Frank B. 1 


SYNOPSIS 
~~. ‘The paper ‘is principally concerned with the vibration of control gates and 
regulating valves commonly in hydraulic structures, as dams. 
The principal findings of 12 yr of vibration investigations by the Army Engi- — 
civil works activities are reported. | 


‘The classical theory of vibration is reviewed to select those portions ap- 
plicable to the subject. 


the results of field tests are analyzed to determine the constant friction force. 
Both laboratory and field tests are — ee 


: _ Emphasis is placed on exicting forces expected to be found in hydraulic 
"structures. Von Karman vortex trail, self-excitation involving reflecting 
_ pressure waves and other hydraulic pulsating phenomena are examined. Cer- -— 
tain conclusions are drawn. 


‘The study a vibration theory has had an interesting history. ~The physical 
phenomena involved with musical instruments was of interest to the scientist 
of the eighteenth century. Lagrange? solved the physical problem. 


Note, —Discussion open unti August . To extend the closing ate one month, © 
a written request must be filed with the Executive Secretary, ASCE, This paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 87, No. HY 2, March, 1961, 

_ 1 Chf., Hydr, Analysis Branch, Hydraulics Div., U.S. Army Engr. Waterways Experi-_ 


- “Recherches sur la Nature et la Propagation du Son,” by JL L. Lagrange, Couvres — 
de Lagrange, Vol. 1, Paris, 1867, 39 (Miscellanea Taurinesia, ‘1, 1759), 


a 
| 
q 
— 
— 
ag 
— 
— 
a 


“of the a string in 1759. A — advance to the physics © of vibration was ' 
- contributed in the nineteenth century by Lord Rayleigh. 3 3 - Many ‘refinements to -_ 
‘Ray leigh’s work have since been accomplished. 
a It is the purpose of this paper to emphasize to civil engineers, the impor-_ 
tance of certain types of mechanical vibrations inhydraulic structures. These 
particular vibration problems are not adequately | covered, evenin mechanical 
—Special is placed on the problems of exciting forces” 
which arise from certain hydraulic phenomena, 
2x The vibration | problems arising in connection with the operation of gates 
ae valves commonly employed to control and regulate the flow of water 
oT through or over dams will be stressed. The word operation is used because | 
much research needs to be done before a hydraulic structure can actually be 
; designed with confidence in the knowledge that the mechanical elements in- 
volved will not vibrate under certain conditions of operation. ; 
The | writer is associated with prototype investigations to determine causes 
=a of vibrating gates or gate elements. Therefore, this | paper. will constitute, in. 
a sense, a summary report of the U. S. Army Engineers experiences of un-— 
desirable | vibration, which have been brought to the attention of the Waterways — 
_ Experiment Station, The problems of vibration in hydraulic machinery —?_ 
io be | treated, nor will the problems of ‘water hammer and surge tank design. | 
APPLICABILITY OF CLASSICAL THEORY — 


fe is s important to examine the classical theory of vibration to determine 
which portions of the science are useful to the solution of the specific prob- 
te lems which» h might be encountered in hydraulic structures. The theory may 
_ ibe found in various publications st such as those of C. R. Freberg, E. N. Kemler,4 
P. Den Hartog, 5s, Timoshenko, D. H. Young,6 M. ASCE, S. Jacobsen, 
A definition of “symbols is shown in Fig. 1. The case of undamped free 
7 ‘vibrations is used for the Ms of defining basic physical quantities. The _ 
device shown is a mass, defined as W/g, with an elastic suspension having 
a spring ‘constant K, expressed in 16 per in. In this case, assume a cube of 
_ steel supported by a tension spring. There will be an initial Static elongation : 


If the cube is forceably pulled downward and released, a vertical vibration 


ensue. The basic differential equation iS 


The symbol x ‘signifies the e second d derivative of ‘displacement with ‘respect to 
time. This derivative is recognized as acceleration. The first term is, _ 


= 3 “The Theory o of Sound,” by J. .W.S, Rayleigh, First Amer, Edition, New York, Do- 
4 “Elements of Mechanical Vibrations, ‘Freberg and E, Kemler, ‘Second 
John Wiley and Sons, Inc., New York, 1949. 
Mechanical Vibrations, by "J.P. Den Hartog, Edition, McGraw-Hill Book 
Coy inc., New York, 1956. and DH Yous Thi 
6 “Vibrations Problem |! in by S. Timoshenko D Young, Third 
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: therefore, the product of mass and acceleration, which represents the inertial 
force of the vibrating body. The second term has the product of the spring > 
constant and the displacement, which is referred to as the siete force. 
useful solution to the basic equation is 


denote the natural frequ ncy of the petted sapien, The symbol w is ; commonly 


ealled the angular frequency although various writers use different symbols. 


Viscous Damping 


1.—DE FINITIONS- UNDAMPED FREE IG, 2, —DAMPED FREE 


The natural period 


“since ‘the frequency in per’ second is the the reciprocal of the period. 


= 


It may be readily seen that the natural frequency is inversely proportional 


tee the square root of the ‘initial static 6s. ‘For the elastic 
described, this value is expressed 


The satin seth is thus hangar to be an important characteristic of the | 
vibratory system. Some sources give: tabulations the constant for 
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; ‘The at angle circle is included to demonstrate the initial displacement 


position. _ Certain writers, notably _ Jacobsen and Ayre, 7 have developed 
va graphical phase plane solutions which appear to be useful in certain aypen of | 


damped vibrations, 


Damped Free Vibrations. —When a restraining force is imposed on a vi- 
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brating sy ystem, the type of damping must be known before analysis. The basic» 


differential equations for the three common types of damping are: | 
Constant friction (no ‘single complete solution) 
§ 
28 K+ KX + Ce (sin 
g 


y= 
We: 


A cos wn t). 


"important to this ‘eager. It should be usted that each of the differential equa- 


tions contains inertial force involving mass and acceleration plus vanseerienll 


“opposes the motion and, thus, changes sign at each half cycle. Therefore, no — 
_ single complete solution can be written (Eq. 8). Each part of the solution must _ 
be considered between largest. positive and negative displacements of 


_ The second term of the viscous damping equation (Eq. 9) contains the 

: damping» coefficient © which is multiplied by the first power of the velocity. 

ot The general equation has a relatively simple and complete solution. The first 

pt factor on the right contains an exponential function which involves the damp- | 
ing constant Cy and mass. The second factor contains a cosine function of i 
Hf the natural frequency. A large part of the classical theory of vibration is based — 
on viscous damping. Unfortunately, the viscous damping solution does not 

ge to be generally to problems head hydrau- 


tl Turbulent damping (Eq. 11) is the third case to be c considered. The turbulent 
_ damping force is a function of the square 0 of the bneenenag A ‘This is the common 


gree differential equation and no direct simple solution. Milne8.9 has 
obtained an elaborate exact solution. He has also prepared tables for the 


evaluation of the solution which have been found useful in surge tank problems: 


“Damped Vibrations W. E. Milne, Univ. of Oregon Publication, Vol. No. 2 


1923, 
“Table of Damped Vibrations, ” by W.E. Milne, Mathematical Series, Vol. 1, No, 
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volved is beyond the 
“scope of this paper. As 
- The pr oblem of attenuation of ‘displacement in damped free vibrations is con- 
—Sidered next. Constant friction characterized by an arithmetic decrement. 
as indicated in Fig. 2(a). It is important to note that only the constant friction 
force and the spring constant are involved in the determination of the arithme-_ 
_ tic decrement. This is an important relationship as will be shown later. eo 
In the case of viscous damping, the attenuation for. each cycle is called the 
logarithmic decrement and can be evaluated when the viscous damping co- 
efficient, the mass, and the frequency are known. _ The attenuation is affected 
by an. exponential law so that the displacement i is ‘damped rather rapidly in the — 
= , as demonstrated on Fig. 2(b)._ After the damping factor Load 


negligible, consider that a steady state exists. Reference will be made to Steady 


- 3. FORCE VIBRATIONS- ‘FIG. 4. —FORCED VIBRATIONS- 
‘VISCOUS — CONSTANT FRICTION 


th exciting force which can increase the energy and hence, the amplitude of 


= 


the vibrations. The basic equation for forced vibrations with the classical — 
iscous damping is 


i 
: ‘The solution shown in Fig. 3 is is fort the case of meine state vibration. In other a 
words, the damping effect of the exponential term has become negligible. The 


F is considered to be applied with a cosine variation, The 
solution for the steady state vibration is 
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The first factor r of the e solution is witeetsiaiaes F/K, ‘which represents the @ static 
"elongation that would be produced by the maximum exciting force with the. 
| existing spring ‘constant. If the cosine function in the numerator is placed 
= to 1, the factor one would be divided by the radical, which is ordinarily _ 
_ known as the Magnification Factor. Note that the Magnification Factor involves , 
of to natural” for the various 


The. critical damping coefficient involves only the 1 natural frequency of the 
vibrating system and its mass. 

In considering the subject of ‘damping friction, it would be appropriate 
Jacobsen and rel0 as. follows: “Friction is therefore, a a 


A graph c can be. constructed dof the Magnification Factor as a function of | the 2 
frequency ratios. The various damping r ratios represented by a family 
of curves. For a damping ratio of zero, ‘the curve represents an undamped > 
_ system, Near the condition of resonance where the ratio of frequencies ” 
equal to one, the Magnification Factor and the amplitude approaches infinity. — 
Dy _ The system is safe against wildly fluctuating vibrations in the resonant range 
the damping ratios are substantially greater than0O.3. 
_ Constant Friction.—As previously mentioned, the differential equation for 
the case of constant friction has no simple continuous solution. . Nevertheless, 
Hartogl1 developed a remarkable solution for maximum 


2 w 
a The solution again contains a ratio of 7 force to spring constant. The > 
- variables A and B under the radical are two different functions of the ratio of 
forcing frequency to natural frequency of the system. The second term under > : 


the radical contains a ‘square ¢ of | the ratio of the constant friction force to the to oll 


“Engineering Vibrations,” by L.S. Jacobsen and R.S, Ayre, First Edition, 


Book Co., Inc., New York, +1958, p. 196, 


ficient which can be evaluated according 
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exciting een The magnification can be expressed as as the ratio of the maxi-— 


i ‘When ‘the second term under the radical (Eq. 15) becomes larger thanthe Jf 
first, the evaluation of maximum displacement involves the square root of 
negative number, Eq. 15 is only of value in the determination of magnification 
_ above the dashed line shown in Fig. 4. The family of curves represent i 
ratios” of the: constant friction to maximum exciting force. For 
ratio forces greater ‘than 


ing high force ratios. yo is ‘sufficient in 1 the region 1 of high force 
ratios to cause a stoppage in the vibratory motion. It is interesting to ie 
_ that the two lines representing force ratios of 0.85 and 0.95 from the 
Hartog investigation have maxima at frequency ratios substantially ae - . 
than the resonant frequency. It might be expected that, in the case of constant _ 
- friction, if the force ratios are considerably greater than 0.8, there would be 
little likelihood of developing large vibratory amplitudes. Inthe case of viscous - 
damping, we surmised that the situation would be sted safe if the damping | 


The engineer is then faced with the problens of determining the » constant 
- friction force and the exciting force. Some test data are available on friction 
factors, which should be applicable to hydraulic gates. On the other hand, 


: experimental data on the eaeaee « of exciting f forces in the various hy draulic 


structures will be presented 


_ Hydraulic Gate Friction.—A substantial amount of work has been conducted 
on sliding x friction and rolling friction research. However, , much of the recent _ 
experimentation has been concerned with the effect of various types of lubrica-_ 
; tion on journal friction. Itis believed that ordinary lubricants are not effective 
cover a long period of time when used on submerged high head hydraulic gates. 
_ The U. S. Bureau of Reclamation!2.13 has conducted tests on both sliding 
and rolling friction, in connection with their studies of gates for hydraulic 
structures. These laboratory tests were conducted in the dry rather than w with — 
water surrounding the element as would be expected in an outlet works, oe 


example. It appears that the coefficient for Ee | is in the —— 


high head gates, provision must be made for adequate capacity to break the 
seizure of metals after the gate has been closed under a high head for a 
substantial period of time. Therefore, with application of a judiciously a 


factor of safety 


the coefficients mentioned would probably be entirely ade- 


q 001 depending on the type of metals used. In the field of hoist design for. 


\ Hartog, Transactions, ASME, Paper APM 53-9, presented at the June 1931 National Con- 
' — of Applied Mechanics at Purdue Univ. in Lafayette, Ind. in tee es 
12 «Tests on Rollers,” by N, G. Noonan and W, H, Strange, U. Ss. 

Tech, Memorandum No, 399, Separate No, 26, September 26, 1936, 
_ “Report of Tests on Coefficients of Friction,” by N. G, Noonan and W. H, Str Ww. 
‘S. Bur, of R Reclam., Tech. Memorandum No. January 30, (1935, 
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quate for purposes of hoist capacity design. On the 1e other hand, : it is believed 
that further investigation be made of friction forces in the prototype to 


bs been found to be practical and useful in ina ‘a study c of exciting forces. The 
Wate rways Experiment Station]4 made studies for the Omaha District beginning _ 

j in 1947, _ These studies simulated the early Garrison control gate design as 

well as the Fort Randall control gates. 

: ba: mass of the gate in the prototype was simulated in the model. An ex- 


a 
Mode! Studies .— —Hydraulic model studies of amass on an elastic suspension 


tension spring was selected which produced a a static elongation similar to that ’ : 
in the care was taken to reduce friction in 


_» 2 VORTEX TRAIL FROM DOWNSTREAM EDGE. 


FIG. 5.—FLAT BOTTOM 6. —DEVELOPMENT OF STANDARD 45° 


>= 


One of ‘the inten of the tests was to determine whether a vortex trail 
pemen a partly open flat bottom gate could produce vibrations. The basic theory 
of this investigation is demonstrated in Fig. 5. The initial static elongation of 
the model gate on its elastic suspension b= 0.0126 Eq. 5, 
natural frequency of this system fn = 8 cycles per oy ees lia, 


vortex trail is from each of the plate. “The frequency of the 
of these vortices is defined by the Strouhal number (St) which involves the fre- 
“quency, _ the length of the plate, and the velocity of the fluid, 


=z “Spillway and Outlet Works, Ft. Randall Dam, Missouri Rive River,” ? Tech, Rep mayer No. 
2-528, Ss. Army Engrs., Vicksburg, Miss., Oc tober, 1959, 


668 March, 1961 HY 2 
il | (©) STANDARD 45° 
4 
ig 
J 


VIBRATION 


Fig. 5 shows the flat plate to have half of its length protruding ¢ into a gate, 
chamber. This represents the partial opening of a control gate. The image © 
of the flat plate above the roof of the conduit completes | the normal concept of 
the vortex trail. During the early phases of the tests it was not known whether, 

_ the vortex trail can be shed from half of a plate as shown in this situation. 
How wever, ‘it was soon learned that this phenomena can exist for the eo 


of a gate leaf into a conduit. ag 
_ Various experimenters Dh have determined the Strouhal number for a flat 
_ plate normal to the flow to be approximately 1/7. If the vortex trail can pro: a 
duce an exciting force on the bottom of the gate which wouldcome into reson-. 
_ ance with the natural frequency of the elastic system, it would be expected that 
the gate would produce vibrations of large amplitudes. For a 17-ft prototype 
gate opening, L,, = 0.333 ft, and a model velocity 18.5 fps, the theory in- 


dicated that there. would be essential resonance. For a normal plate St = 1/7 


It was actually determined that when the model gate was set to simulate a 17-ft 
_ prototype gate opening, violent vertical vibrations were produced. Therefore, 

: it was 0 obvious that resonance between the frequency of the vortex trails from 
a gate projection acting as an exciting force on an elastic sy stem could ieee 
be produced in the laboratory. 

Split Leaf Gates.—In connection with the subject of the vortex trail as an X= 
citing force, mention should be made of split leaf spillway gates. With this sys- _ 

tem, one leaf rests on top of the other. In times of flood-flow both leaves are 


After the top leaf has been removed and the bottom leaf is in the process ss of | 


‘raising, water can flow over under the partially leaf. The full vor-. 


bottom of the gate — 
Old River Control Structure gates. 

Improv ement of Gate Lips. i er from the studies of Ft. Randall 
gate vibration that a flat bottom gate had high downpull. The tests at the Water- | 

_ ways Experiment Station further showed that a flat bottom gate also was more 

susceptible to vibtations. lip | extensions from the downstream edge 


nk from ‘the are shown in Pig. and 
A fairly regular motion of 8 cycles per | sec may be seen in Fig. 6(d), whereas 
no such regularity is shown for the case of a lip extension (Fig. 6(e)). 
Py The vortex trail was probably springing essentially from the lip extension 
rather than from the. upstream edge of the flat bottom gate. . The engineers of. 
_ the Waverways Experiment Station tried the next logical step which was a gate 


with a 45° bottom surface (Fig. 6(c)). It was found that the vortex trail would | 


Clarendon Press, Vol, 2, 1938, p. 571, 


= 0.143 and 1 ff = = 7.92, thus the ratio of forcing frequency to frequency 
« 


a 


q 


—a “Modern Developments in Fluid -Dynamws,” by S, Goldstein, First Editio Edition, , Oxford, 4 


~ 16 «Downpull Forces on Vertical Lift Gates,” Tech. Report No, 2-477, Report No, : 
U. S. Waterways Experiment Sta., Vicksburg, Miss. , December 1908. 
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spring from the downstream edge and minimize pulses: on the bottom of 


for vibration was found for the 450 gate lip. Furthermore, the 


is indicated and there is towards vibration of t the gate. 
Prototype Tests.- — Before proceeding with a report of the prototype tests, it 
a: appropriate to briefly examine the important problem of the spring constant 
for cable supported control gates. Tests were made by Waterways Experiment 7 
Station engineers at Knightville Dam in New England in the fall of 1955. In the 
_ analysis of the data, it was desirable to compute the spring constant of the cable ; 
supports. the total metallic cross-sectional area of the cable and the modulus 
of elasticity are known, , the spring constant can be determined. A fairly accu- 


rate value of metallic cross- sectional area can furnished by the manu-_ 


BP only information on ‘on modulus of elasticity of wire dept at that time yrs 
by 
one of the wire rose manufacturers. ‘The use e of this value did not es ; 
theoretical frequency which agreed with the measured frequency at Knightville 
= Dam on tests at three different gate openings ranging from 10 1/2 ft to11 ft. 
a _ The natural frequency measured at Knightville was about 6 1/2 cycles per 
} sec. . Knowing the other physical data pertaining to the vibrating system, the 
+ modulus of elasticity was 3 estimated | to be approximately 21,000, 000 psi. , Soon 
_ thereafter, information became available from tests on prestressed wire rope. 
_ These results indicated that the modulus of elasticity should be in the neighbor-— 
a hood of 20,000,000 psi. This value is now considered to be a fairly reliable a 
7 modulus of elasticity for wire _— which has been stressed repeatedly i in oper- 


Le It was clearly evident ‘that the mass of the water contained within the gate 


-must be added to the mass of the gate metal for a proper vibration. analysis. 7 

_ Fort Randall Tests. The Waterways Experiment Station!” conducted ase- 
"ries of three tests, for the Omaha District, 7 the es gates for Ft. Randall — 

in (1954 and 1 1955, while the reservoir filling. Fig. ' 7 pertains to 

test sade with a head of 110 ft. The s ‘spring constant for the cable suspension _ 
1.7x 10° psi and the weighed 47 tons. 


large ‘motion | of the gate ry upstream direction. Half of the ¢ guides were 
‘spring mounted and the other half were 
‘The oscillograph records show that for large ¢ gate openings, there were oc- - 
casional vertical vibrations of approximately 4 cycles per sec. These vibra-_ i 
tions continued for several cycles and were then damped out. After an exam-_ 
ination of the record, it was apparent that the frictional force involved in the 
; rollers was rather high, but that occasionally a surge against the downstream 
ai of the gate moved the 1 upstream so that momentarily the rollers were 


Pas 


came in | contact, , the vibration was 


River,” Tech, Report 1 No. 
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a - less for the sloping gate bottom than for the flat bottom. erefore, the 4 z ae 
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that ae damping v was of a character. ne of the damping records is 
shown in Fig. 7(b). It may be > noted that straight lines connect the peaks and 
troughs ofthe curve, 
ey This was a significant finding, he: although experimental evidence was at 
hand to show that the damping was predominantly caused by a constant friction 
force, the _ theoretical problem was known to be complex. It is pertinent to 
quote again from. Jacobsen and Ayre18 as follows: 
“It is therefore clear that the presence of constant friction greatly com- _ : 
plicates the analysis of steady state response to an alternating force, ie 
‘especially if if the friction is large enough to enforce > Stops in the motion.” 


arithme ic decrement was then computed for damping 1g periods. 
_ Knowing | the spring” constant and th the arithmetic decrement, | the comping force 


The ¢ section of in Fig. 7(b) is fora a 20 ft gate opening the 
damping force was computed to be 2,740 Ib. _ An average of four age 
(gate openings of 19 ft to 21 ft) indicated a damping force of about 2, 800 Ib. 7. 


should be mentioned that in two cases studied for the same gate opening, one 
_ damping force appeared to be double the damping force for another period. It 
is possible that the lower force could have been ci caused by only one roller track 


‘coming in 1 contact - with the guide. Int the case e of the higher force, it appears 
that both roller tracks came in contact with the guides simultaneously. a 
tt It was realized that although the theory of Coulomb damping is complex, , the 
simple determination of the arithmetical decrement affords a means of n meas- 
uring friction. It is believedthat a study of damping records can furnish valu- 
able information on the total damping force for any particular system under bas 
observation. There has not been an opportunity to attempt to compare the total | 
damping force measured in the prototype with a computed force based on fric- 
tion factors as measured in the laboratory. An estimation of the frictional 
_ forces of the individual components of a rollertrain may not be a simple mat- 
Sluice Oates. — The hydraulically operated sluice gate is a common device © 
used for controlling the flow of water through concrete dams. _ Prototype t tests” be 
were made on the sluice | gates of Pine Flat Dam, Calif. for the Sacramento 
District in 1952. Measurement of the vibrations of the gate leaf was a part of : 
_ these tests. The gate size was 5 ft k by 9 ft. Fig. 8(a) st shows the schematic sec- 
Bs The 6 ton sluice gate was supported by a 6 1/2 in. diameter steel stem. A 
typical oscillogram is shown in Fig. 8( (b). The maximum amplitude of vibra- a 
5 was of the order of 6/1,000,000 in. This” emphasizes the fact that we are 
dealing with a relatively light weight gate on a stiff suspension. SON Pen ee inl 


The in Fig. 8(c) represents a significant observation on the 


“Engineering Vibrations L.S, Jacobsen and R.S, Ayre, First Edition, McGraw- 


Hill Book Co., Inc., New York, 1958, 226, 
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. It seems ‘reasonable t to believe that this larger vibration am- 
4 plitude was ‘caused by alternate seizure and relief of the metals on the slide 
bearing. This is sometimes called chatter, | 
a Although the hydraulic head for the test it just described was 110 ft, this gate 
was subsequently operated under aheadin excess of 300ft. No detectable dam-— 
age to the gate was cursory study was made of the fatigue as- 


‘taate: any further consideration of fatigue. No study has been made of the pos- - 
: sible fatigue om aa the supporting elements at either end of the stem. 4 


Flexural Vibration.— The horizontal beams ofa gate are normally 
stiff. The Fort Randall gate beam is shown on Fig. 9. 
sa a For a concept of the stiffness, it may be noted that the initial static deflec- - 
tion based on the mass of the gate and contained water is only 6 x 107 4 in. A 
_ comparison is shown between the theoretical frequency of the beam with an with- 
out contained water, 
- The frequency for the dry condition may be seen to be 205 cycles per s sec 

and 140 cycles per sec for the condition of inclosed water, 


the importance of accounting for the mass of the contained water for gates of 


Howe ll- -Bunger a years ago it wa was learned that weld seams on : 
4- -vane Howell-Bunger valves at two different projects had failed. Fig. 10 isa 
schematic representation of the vibration characteristic of both a 4-vane and a ; 


-vane ne valve Timoshenko treats the subject of vibration | of thin walled 


_ Fig. 10(a) shows schematically a 4-vane valve shell vibrating with the fun-_ 
damental mode and Fig. 10(b) a 6-vane valve shell vibrating with a secondary — 


a - mode. | Eq. 21 gives the natural frequency for the fundamental mode and any | 


It n may be seen that the ‘eaee oniinitiie radical i is concerned with the m mass , 
of the factor r (Mj) is a mode factor. 


‘the ‘mode ( Fig. 10(a)), i= 2, and the resulting mode factor 
_ = 2.68. For the secondary (Fig. 10(b), i= 3, and the value of the mode 
factor Mg = 7.59. Since the naturalfrequency of the shells are proportional | to 


ratio of the fundamental mode is 2.82. 
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‘This ignores the existence of vanes within 1 the e shell. " However, it appears | 
that by the use of 6 vanes and producing a vibration of the 1e secondary r mode e type, 
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the frequency could be ‘expected to be much higher than for the 4 vanes s and 7 
- fundamental mode. When the frequency is substantially higher, it would nor-_ 
mally be expected that the amplitude of vibrations would be lowe er and that for 
each cycle the total amount of energy, is less for the higher modes. For this 
- reason, it is believed that a 6-vane valve is less liable to know failure of the 


seams which connect the vanes sto the ‘shell. 
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To the wri the writer’s » knowledge, the problem of the Howell- Bunger valves with 

_ the mass and elastic characteristics of both shell and vanes has not been 
7 _ treated. It would be interesting indeedto see a thoroughanalysis of the 4-vane : 
valve compared with the analysis of the 6-vane valve. Furthermore, it is con- : 
ceivable that a 5-vane valve would cause the vibration to go toa still higher | 
mode and hence have less liability for failure of weld seams. aa 7 ee oe 

ae: Rayleigh-Ritz Method.—It is appropriate to mention the Rayleigh- Ritz method 
by which problems of this nature can be solved. _ Lord Rayleigh3 treated only ‘ 

the ¢: case e of the fundamental mode. . He devised the method of writing the basic 
equations for both potential energy and kinetic energy. _ By equating these two 


quencies of the various types of elastic systems. Some years later ‘Walther 
 Ritz19 extended Lord Rayleigh’s method to treat the problems of the vibration 
characteristics with modes higher than the fundamental. _ Since that time, nu- 
2 -merous contributions have been made to the theory of vibration of elastic ele- 24 
ments using the Ray leigh- Ritz method. = 
- Gate Skin Plates.—The vibration of the skin plate on a submergible gate at a 
Cheatham Dam on the Cumberland River offers an example of this type of vi- 
fe motion. The operators reported a loud noise issuing from the Cheat- 
ham gate when the overflow head was about 3 ft. The gate is shown in Fig. 11. 7 
_ The Waterways Experiment Station investigated the cause of the noise in 
cooperation with the Nashville District. Accelerometers were placed to meas- — 
= the principal freedoms of motion of the gate itself as well as those of the | 
upstream and downstream skin as The specific source of the noise was not 
_ It was found that the e frequency of the downstream. skin plate approximated 
the frequency of the sound. The frequencies of other freedoms of motion were q 
substantially different. It was therefore concluded that the downstream skin 
_- The resonant air chamber inside the gate was analyzed b’ but it found to havea | 
frequency in comparison. Similarly, the vortex trail from the strut arms 
- possesseda low fr equency. As the exciting force of the vibratory motion could © 
not be readily determined, itwas recommendedthat the skin plate be restrained : 
bya support in the middle of the panel. it was reasoned that this would greatly 


increase the natural frequency of the skin plate and eliminate the vibration which 


“Self Excited Oscillations. The action of the bow ona violin string is the 
classic example often cited. good example of, self-excitation from the hy- 
draulic engineering viewpoint is the interaction of a reflecting pressure wave 


conduit and vertical a . This is: shown schematically in in 


small vertical of the gate wont cause a pressure wave to. 
travel upstream in the tunnel. When the wave reaches the reservoir, it is re- 
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i If the speed of the wave is is a and the distance from gate to reservoir is L,the 
period of the pulsation against the gate bottom would be L/a sec. The pres- 


_ the sake of simplicity. When the natural period of v vertical vibration of the gate 
‘is close to the natural period of the pressure wave ‘condition of self- 


sh 


Free Water Surface Phenomena .—Two free. water sur face phenomena | are 
own in Fig. 13. Either could constitute an exciting force on an elastic str uc- 

_ The fluctuating nappe for a ‘small head ona sloping weir has been noted by 
various my pon This phenomena was extensively observed and analyzed by 
- Bruno Leo. 20 The phenomena was also observed by the Bureau of Reclama-— 


tion21 on the drum gates at Black Canyon Dam. In this case, the vibration was 


by under the nappe. 
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PRESSURE WAVE 


Channel Waves. —The binodal wave which has. a clapotis action f from 
side wall was “observed in the fish ladders at McNary Dam. OA trinodal 
wave was also observed in the Bonneville fish ladders. These waves can con- 
stitute an exciting force on adjoining structures. = 
Other Exciting Forces .—Various other phenomena can possibly constitute 
_ exciting forces in vibratory 1 motion. The successive formation and collapse of 
vapor cavities in the phenomenon of cavitation seem to have a periodicity. ; 
Little experimental information is available for the multitude of possible 
metrical boundary situations. . Some experimenters have suggested a Strouhal 
. number which defines the frequency of shedding of cavities. 


20 “Self Excited Vibrations at Overflow Weirs 


1 «Report on Vibration Studies Made at Black Canyon Dam,” by R. E. Glover, C. W. 

F, Hammett, U. Bur. of Reclam, Hydr. Lab. Lab. , Report No, 58, Denver, 
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Tn eee cavity which forms on the trailing side of an obstruction | in high > 
velocity flow has an intermittent pressure pulse. However, little research has 

been conducted on this problem. 

- ‘The toe of the hydraulic jump is known to pulse in an up and downstream — 


direction. The possible effect of this is unknown. 


by measurement of vertical vibrations. 
_ 2, The type of friction whichaffects the damping of hydraulic gates has been 
shown to be constant friction or Coulomb damping. = — 
3 The hance of the ' constant friction force can be evaluated by an anal- 

: 4, The Von Karman vortex trail has bee been observed to be the exciting force 
which can cause vibration of gates. 
5. A gate leaf with a sloping bottom sheds the trail from down- 
edge and therefore minimizes vibratory motion, 4 
Extensive research is neededon the magnitude of the exciting forces be- 
gates can be designed against vibration. 

=. Considerable research is needed on the character rol vibration of both fixed 7 


a “a traveling ca\ cavities and other hydraulic pulsating phenomena, = 


a by Leiland M. _ Duke a George C. , Downing ¢ of the Instrumentation Branch at © 
the Waterways Experiment Station. It was only - with the cooperation of many © 
engineers in the Division and Districts of the Army Engineers, that the field 
tests were possible. _ Engineers of the Waterways Experiment Station, too nu- 
merous to mention, also participated inthis work. E. P. Fortson,Jr., F. ASCE, 
is Chief, Hydraulics Dibision and J. B. Tiffany, Jr. , F. ASCE, is Technical Di- ms 
rector, Waterways Experiment Station. — The field work was performed over a 
period when C. H. Dunn, F. ASCE, P. Rollins, Jr., F. ASCE, and H. 
F. ASCE, were successively directors at the Station. _ 
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“HYDRAULICS DIVISION 
aaa: of the American Society of Civil Engineers 


STREAM- NETWORK IN THE UNITED | ED STATES?, 


The national stream- gaging program is with regard to size, 
comparison with other countries, and between states, Its deficiencies are pin- 


of areal stream gaging is described that would overcome some of the existing 
deficiencies. Progress made toward ird adapting such a aig! is 3 reported and 


INTRODUCTION 


At least four thorough studies2,3,4,5 have been made since 1949 to evaluate 
“the adequacy of hydrologic data in this country. The strengths and weaknesses 
of existing data programs were analyzed at great length and a number of rec- 
were made for needed improvements. of these 


s Note.—Discussion open until August 1, 1961, To extend the closing date one month 
_a written request must be filed with the Executive Secretary, ASCE, This paper is part 
“of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 87, No. HY 2, March,1961, 
publication authorized by the Director, United States Geological Survey, Dept. of the 
b Presented at the March 1960, ASCE Convention in New Orleans, se t ¢ 


a Dist, Engr., U. S. Geological Survey, Trenton, N. J. 


_ “ Adequacy of Hydrologic L Data to Meet — al Needs,” ” Federal Inter- Agency River = 

“A Water Policy for the American People, ” ‘President’ s Water Resources Policy 
Comm., U. S, Govt, Printing Office, Washington, D.C.,1950,. 


yeaa 4 “Water Policy,” Presidential Advisory Committee on Water Resources Policy, U. 


= Printing Office, Washington, D.C.,1955, 
’ “An Engineering Appraisal of Hydrologic Data,” Task Group on Hydrologic Data ata of | 
the Committee on Hydrology of the Hydr, Div., Proceedings, ASCE, Vol, 85, No. BY 7 
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also a to estimate the required size of program or number * stations a ; 
needed to provide adequate hydrologic data of the several principal types. The 
estimated number of stream- ~gaging stations needed, as determined in these 

three surveys, is shown on Fig. 1. Also plotted for is actual 


number of gaging stations showing the growth rate since 1903. 


7 *) _ It can be seen that all three of the surveys concluded that the number of © 


gaging stations should be expanded rapidly. The surveys differed only slightly 

as to the total number of stations recommended by given target dates. ‘Never- 

theless, the growth rate of total - stream-gaging stations slowed down about a 
1952 in the face of these carefully considered recommendations. To under- 

stand this paradox one must look at the existing program from several different _ 


7 =e and then follow the evolution of new concepts aimed at getting the 
‘maximum amount of hydrologic information for each dollar eee in stream © 


PROGRAM FROM SEVERAL 


Size and Cost of Program. —Operation costs for 2 a stream-gaging 
in the United States range from $600 per year to as much as $5,000 a year 
- for a few remote stations, with an average of about $950 per station per year - 
; in 1960. Installation costs range from $1,000 to $15,000, and average about 
$3, 000. The average life of a structure is on the order of 30 yr, making de- 7 
-preciation, without interest, about $100 per year. 
= cost to the nation averages $950 for operation plus $100 for amorti- 


¥e on a 50-50 matching basis and nearly one-fourth is transferred to the USGS 
from other Federal agencies. Perhaps another million is spent annually by | 
other Federal agencies within their own organization for for stations closely re- 
as The Mississippi River Commission | annually publishes about 60 records of 
- river stage and discharge and another 100 records of river stage only. 6 The 
International Boundary and Water Commission, United States and ‘Mexico, © 
‘publishes approximately fifty streamflow records collected in the Rio Grande > 

basin.‘ The USGS also publishes some records collected by other agencies. _ 


The total national investment, including state, county, and city expenditures 
q 


Interior (USGS) alone. Of this, about one-third is contributed a the states = 


-made— in n cooperation with | the ‘USGS , is about $9 million each year for basic ; 
‘streamflow data. This amount represents less than 0.5% of the amount spent — 
in the country annually for dams, levees, irrigation works, water supply 
sources, and other works dependent on streamflow data for safe and economi- 
eal design. n.8 Expressed another ¥ way, it costs about five cents per person | 
per year for operation of streamflow stations at all levels of government, This 
xpenditure modest considering extent to which man is dependent on 


er for luxuries, necessities, and indeed life itself. Yet continuing infla- 


? _ 6 “Annual Report of Stages and Discharges of Mississippi River and its Tributaries,” — - 
River Comm., Vicksburg, Miss, 
Bra — T «annual Report of International Boundary and Water Commission, United States and 
8 «Water Facts for the Nation’s Future, Ww. Hoyt, Ronald 
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tionary trends and mounting Federal it to economize 
even more if that is feasible. 

Comparison with Other Countries .—How does the ‘United 

‘ with other countries? Fig. 2 shows the number of gaging stations per thou- | 
is - sand square miles in relation to total area for various countries, as listed by 
7 ‘Walter B. Langbein. 9 It is evident that area alone is not a reliable guide to 

the ‘number or density of gaging stations needed. The range is infinite. Libya 

has” no gaging stations in an area of 633,000 sq miles, while Israel has 85 


Ze in an area of 8,000 sq miles. It wales be inferred from this plot that 


station density in “relation to population density is illustrated in 
Fig. 3. Including: population as a parameter along with area makes the re- 7 
lationship take a vague form that can begin to be interpreted. Note that un-— 

‘ developed countries plot low and usually toward the left. The correlation is 
still too scattered to be very useful. However, on the basis of population, sol 7 
. United States does not seem unexpectedly high or low in gaging-station Goneity. 


_ plosion” will presumably continue, the position of the United States might be 
expected to continue to rise and drift further to the right. The change in posi- 

_ tion since 1930 is shown by the dashed line. Of course, most other countries 
are _ tending to: -move in the same general direction at a rate dependent on the 
speed of development oftheir water resources, 2 
_ Comparison of States Within the United States.—Perhaps better evidence 
that the number of gaging stations is related to population density can be 

seen in Fig. 4, as reported by Langbein.* 9 ’ This graph shows the relative posi- 

: tion of each state in continental United States as of 1955, after empirical ad- _ 

_ justment for irrigation and waterpower. The adjustments were determined by _ 
trial, their effectiveness being measured by reduction in scatter of the plotted © 
‘points. Langbein’ S analysis indicated that, insofar as stream-gaging onl 
are concerned, each acre of irrigated farmland was equivalent to 4 people. a 
7 Likewise, , each hydro- power plant of 1,000 kw or more capacity was found to 
be equivalent to 10,000 people. With these adjustments the points tend to de- 
fine a rather consistent relationship between number of gaging stations and 
7 population density. That is to say, that in the United States the stream- -gaging 
program in all states is affected by the same national policies and that the _ 
standard ¢ of living or stage of economic development does not vary appreciably © 
between - states. . New Jersey with 10.5 gaging stations per thousand square > 
‘miles and Nevada with 0.5 station per thousand square miles represent nearly 
_ the two extremes in both population density (2nd and 48th) and industrialization 
in conterminous United States All other states line up reasonably well on 
Effect of Administrative Policies.—It is interesting to note that Alaska, 
. which has replaced Nevada as our least developed and populous state since 
~ Langbein’ s analysis (1957) would plot off the sheet to the lower left but about 


in line with the other states . Alaskais the only state that has never cooperated 


_ 3 “Numbers of Stream-Gaging Stations in Various Countries with Analysis of their 
- Distribution in the United States,” by W. B. Langbein, Bulletin de 1’Association Interna- 
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ie - more intensively than the larger countries, and that gaging-station density | a 
= 6€=—CStsérefilects this in a general way. Note that the United States is considerably _ eg 
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a 50- 50 financial basis, and all stations in Alaska prior to 1959 
were financed ‘completely by the Federal government. Hawaii would plot off 
7 - the top of the sheet, having twice the gaging-station density of New Jersey, 
é the previous leader in this respect, but only one-tenth the population rae _— 
- of New Jersey. While Hawaii is hardly typical, being made up of eight major 
and many smaller islands, it seems likely that the administrative policy of 
the Territorial government was a major factor in reaching such a high level 
in gaging- station density. The huge pineapple plantations probably exerted 
major influence in this because of their irrigation and drainage needs. Hawaii, 
‘ Bi in contrast to Alaska, has for several decades furnished half of the stream- “a 
‘matching of state or 
territorial - funds for stream gaging. This principle of 50-50 matching has 
been followed ever since and the Federal-state cooperation has become the | 


: backbone and ‘major “source of funds for the work. Referring again to Fig. 1, 
a 


_the sharp rise in the curve since 1929 shows very plainly the effect of this | 
partnership policy on growth of stream-gaging 
Israel is another example of a country which , by administrative | policy, 
decided it would have intensive stream gaging and ground water exploration. - 
mh Taking advantage of all feasible means to “lift itself by the bootstraps” from 
4 backward > nation to a position of leadership in the Middle East, this tiny 
country has perhaps the most intensive and best planned hydrologic investi- aot 


gations’ program of any nation in the world (asof 1960), 
—The nny pressures f for gaging stations 


“ment. rather only surface runoff not being measured today in 
“conterminous United States is in small coastal streams and in tidal 
- reaches of the larger rivers. In fact, some water originating in the upper 
Missouri or Ohio Rivers” } may be measured a dozen or more times before it 
a _ However, in most instances this flowing water would not be measured until 
it had passed through a runand abrook, or perhaps a creek or two, , and finally 
entered a river draining several hundred square miles. Less than one-third _— 
of the gaging stations in the United States are on watersheds of 100 sq miles ~ 
or smaller. Less than 5% of all gaging stations are on drainage areas of 10 
sq miles or less despite a three-fold increase in small-area gaging-station 7 
coverage since 1945. Yet there are many, many more miles in aggregate 
length of small streams than there are total miles of the large streams. As 
the years pass, the smaller streams are becoming of ever greater economic, — 
_ Sanitary, and recreational importance. To be able to say with confidence what _ 
runoff is available in these smaller s streams, or the extent of their flooding 
‘potential and other « characteristics, more of them n must be gaged. At the same 
[ time, new technical knowledge must be developed that will let the engineer 
or scientist interpolate the characteristics ofthe still ungaged small streams, 
F: which probably will always outnumber the gaged streams a hundred-fold. This 
on > Ss large order because the flow characteristics become more diverse and 


y variable as the becomes smaller. 
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3 Assume the choice is one station for every 20, 50, 100, 200, 400, or 700 miles 
(of of stream. For small streams » choice of one every 50 miles would be the | 
om as saying there should be one on every tenth stream averaging 5 miles — 
in length. Most engineers quizzed by the writer satuativety chose from 50 - 
100 miles as a reasonable distance between ‘stations, yet even the more con- 


-servative of these two choices would require quadrupling the size of the present a 


wor 


ae From spot sampling of stream density (map studies of perennial stream 
‘length divided by area) in various parts the country, E. W. Coffay10 
found ven average to be slightly more than one mile length per square mile 

of area. Although this study has not been carried to completion (as of 1961), 
and the final answer / may be slightly different, a rounded figure of 3 million» 
miles of perennial streams in conterminous United States can be used safely _ 
estimating puspases. in Alaska and Hawaii but including 


r a It is interesting to 10 ponder | the question of the spacing ing of gaging stations, 


“state agencies (6680), the ratio | comes out ae one station n now in 1 existence for 


“miles a on the average, so it ‘must be concluded that stations a are “e still i 
ae farther apart on the small streams at present. In fact, stations now may aver- a 
age less than one for every 600 miles of aggregate length on the small streams, 
Admittedly, a mileage basis is no way to design a stream-gaging network; ‘its > 
= significance is to lend perspective as to the magnitude of the task of — 
getting adequate data on all small streams teal 
Segments of the Hydrologic Cycle. —An excellent and up-to-date 


i eter of the six principal types of hydrologic data was given in 19599 by * 7 
“the Task Group of Hydrologic Data of the Committee on Hydrology, Hydraulics 
_ Division, American» Society of Civil Engineers (ASCE). It will serve = ' 
a 4 simply to quote excerpts of the Task Group’s conclusions which are pertinent — 
5 et ‘The existing federal networks of gages measuring hydrologic data 
are producing satisfactory records. Periodic reappraisals are desirable 
_ for bringing about improvements to keep pace with changing needs and ae 
= A reluctance to change prevailing practicescanlead to over- 7 
collection and failure to begin records needed for the future.” = an 


The noted several new types of dataof current and future importance. 

_ These included natural or background radioactivity of waters and changes pe 
- £ radioactive poilutants, data to evaluate effects of watershed treatment 


In an earlier part of the Task Group report ; devoted specifically to runoff | 


it t dthat: 
“The obvious is in data from small drainage areas 

f 50 square miles or less. More data are needed on streamflow in the 
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a of tidal influence, and such data already collected should be ; as- 
sembled, analyzed and published.” 


ES AND TENTATIVE 
Looking national stream-gaging | program m several different 


rspectives leads to the followin tentative conclusions: ¥ eo 


United ‘States is” not unexpectedly high or low in gaging-station 
» d density in relation to other countries. Considering the high standard of living 
“ and extent of industrialization in the United States, it may be lower than ex- 
_— On the other hand, the United States has considerably more stations" 
any other country of comparable size. 
The need for streamflow data apparently is related to population 
- extent of industrialization, both of which are ‘projected to grow very 


rapidly in the United States in future decades. 


- ‘streams. Past methods of building permanent gaging stations wherever data 


are needed will not suffice even if funds and manpower could be doubled or 


bs ee 5. To keep pace honest the widespread needs, the program 1 must be directed 


bd on smaller streams and tidal streams. Also some new types of data should 


be observed and | published to meet recently emerging requirements. oo Loe 
The phrase “areal stream gaging” is used to describe such a program, in 

which all available streamflow, flood, and drought records are evaluated along 
- with topographic, geologic, and meteorologic information to determine regions — 
_ or areas of hydrologic similarity. Within such an area the shape of the hydro-— 
; - graph is reasonably similar for streams of equal drainage area; peak runcff 
per square mile, | average runoff per square mile, and base flow character-— 


istics tend to be comparable. Within such an area only one or two complete 
‘gaging stations need be maintained continuously, with variations from place 


to place within the area , being defined by) reconnaissance methods. an 


= REVIEW OF THE NETWORK | 

D velopment of Concepts. —The thinking and concepts previously expressed 

were actually the development of many years of careful consideration and 
c study by many hydraulic engineers. The first person to document these ideas" 
formally, as far as the writer is aware, was Walter Langbein in 1951. The 
ao pros and cons of these concepts were debated extensively within the USGS in 

the months that followed, and in 1952 a committee of six engineers was ap- 
pointed to es this matter formally. This committee, of which Langbein was 
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G 
“tributed within | the USGS for and field trial of the ideas and 


[— proposed. The concepts of the committee were stated as follows: a 


adie. “The printishes that are proposed as s guides in the development of the 
gaging- station program in order to achieve maximum information rela- — 
tive to costs are the following: (1) Apart from operational needs, stream 
_ gaging is a sampling process, and (2) within limits it is ‘possible ps 
relate the flow of one stream with that of another. A gaging station pro- 
vides information on the rates and amounts of flow during a particular 
_ time in a particular channel and therefore samples a given time and ; 
given place. The variations that occur with respect to time at a gaging <> 


- station are fairly well understood, but those that exist from place to ae 


place less so. A well- balanced network , therefore, givest the two factors — 
of time and place their proper ‘consideration. The two factors can r: 


brought into balanced relationship through use of correlation techniques. — 

“These basic principles lead to the conclusion that a record can be e- 
_ viewed not only as a measure of the flow at one site but also as an in- 


gg aan to the flow of streams in the contigrous area—in effect magnifying a 


paw) 
In November, andC. presented a which 
the theory, procedures, , and statistics it involved in n the correlation of 
a a short record with a nearby long- term gaging- -station record. . These tech- 
~ niques are very wsetul in applying the program concept of a hydrologic net-- 
work of stations, The Presidential Advisory Committee on Water Resources — 
“Policy, in December, (1955, recommendations* which added further 
(4 impetus to this concept of a hydrologic network, made up of long-term and — 
_ short-term stations. The USGS ran a pilot projectin the Colorado River basin 
to test the feasibility of the concepts and the problems involved in attempting — 7 
to adapt an existing stream-gaging program to these network principles. The 
problems were many but none insurmountable. results of the pilot project 
deemed to warrant adopting the network as a guide for future 
During the last half of 1956 and all of 1957, a formal and detailed review 


‘2 the entire stream- ‘gaging program in the United States and its territories 
: _ was made by the USGS in accordance with these network principles. Since that 
time the ‘USGS has been following these principles in its program design. 
- Naturally, _ these concepts and the modification in program goals were dis- 
cussed with representatives of other Federal agencies inthe Subcommittee on 
7 Hydrology of the Inter-Agency Committee on Water Resources. They were 
—- also discussed with local cooperating officials in state and federal agencies" 
at different times: during the consideration, testing, and adaptation of these 


sss At no time during this period of more than 8 yr that these concepts have | ‘ 
been developing have engineers of the USGS felt that ‘they have known all the 
answers. Accordingly, there has been no great rush to get into print any final © 
/_version of this type of network Coen. However, it is felt that the concepts 


have been sufficiently tested | and that enough answers are known to warrant rant 


4 il “Extending Streamflow Data,” by W. B. Langbein and C, H, Hardison, Proceedings, a 


ASCE, Vol, 81, Proceedings Separate No. 826, November, 1955, sss 
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"publication, Furthermore, i ‘it has become apparent that lack of wide dissemi- 
2 nation of the ideas and concepts has somewhat deterred their acceptance by 


the engineering profession and officials concerned with streamflow data pro- 


Eze information has been presented by Langbein and Ww. G. Hoyt. 


‘This book, sponsored by the Conservation on Foundation, co contains a compre- 
“hensive discussion data networks and programs. Therefore, 
, procedures, — 


Streamflow Station Stage Station Partial | Record Station 


ydrologi ic Network | Hydrologic Network drolog ic Network — 


Primary (long term) Primary (long term) 7 
¥ “Areal (including ne ne- Areal (Ponds & lakes) 


A ‘12 Mainstream Mainstream (River 
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Classification of Stations. _The first step of the 

a program was to review and classify all existing gaging stations as to their as 
best use in a nationwide network. Table 1 _ shows the classification | system a 
adopted. Each main category was further divided into hydrologic network or 
water management stations. A water management station is defined as one © 


needed to pene data for one or more specific operational, legal, or admini-_ 
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7 strative purposes. Its es establishment, location, and length of | operation depend -~— 


and general purpose needs. The programming problem was thereby narrowed — 
+ down immediately to the design of a general purpose hydrologic network. ¥ a 
_ However, in all cases for which a ‘gaging station was found in the original 
: review to serve hydrologic network purposes as well as certain water man-_ 
agement purposes, was as classified as as a hydrologic network station 
_ The hydrologic network was broken down into primary stations and sec- ~" 
ag ondary stations . The primary stations are the long-term stations to be oper- 
ated for an indefinite period of time. ». Secondary stations are intended to be - 
of short-term operation, say from 5 yr to 10 yr. At the end of that time the 
‘Station, and as much of the equipment as can be moved, would be relocated to | 
other point needing stream gaging. Both the primary and stations, 
classifications were further broken down into areal and mainstream stations. 
i Primary Stations.—Areal primary stations are index stations on streams © 
that must essentially be free from past regulation and, hopefully, free : from 
extensive future regulation, diversion, or other development. In parts of the © - 
United States it is becoming increasingly difficult to find stations completely y 
free from regulation and diversion. This made the job more difficult but no- 
where prevented the use of the network principles. Thus an areal primary 
station is selected for representativeness and length of record; insofar as 
possible it is free from past and future regulation or diversion, and will be 
Operated for an indefinite period of time to obtain a long range time mae 
; of the hydrology of the section in which it is located. This station record would 
"4 be used as the independent variable in making correlative estimates of long- 
; term streamflow characteristics at other sites in the same hydrologic pro- 
A mainstream primary station, on the e other hand, may be considerably — 
aun by regulation or r diversion, It serves only as as a record of actual © a 
‘flow at that point and as an index of flow at other points upstream or ‘down- Es 
stream on the same large river. An example of this would be the Mississippi 
River at Vicksburg, at which point the river flow would correlate very closely - 
with simultaneous flow in the river at other points upstream or downstream ss 
for many miles, but would not bear any particular relationship to the flow in » 
small tributaries entering the Mississippi River inthis vicinity. 
_ The areal primary stations selected were those which seemed to be most = : 
7 representative of the hydrologic characteristics of the area in which they 
= located. ‘Trial correlations of monthly mean flows for a 5-yr period, _ 
_ generally 1951-55, were used for this purpose. It was found that stations with . 
_ either very large or very ‘small drainage areas did not correlate well with — oF 
surrounding stations. Therefore, most areal primary stations, asa result, = 7 
are Stations with medium size drainage a area, . Attention was given to this 


: and some on larger drainage areas in all parts of the country in order to 
_ avoid complete bias in the final 
Secondary Stations .— —An areal secondary station would be operated where 
en streamflow information is wanted or will likely be needed in the | 
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A mainstream: ‘secondary station, likewise, would be operated for 
only aS many years as needed to define a correlation with a nearby main-— 

stream primary station. Then either type of secondary station might be dis- a 

continued or converted to a partial-record station. Because of their frequent 

relocation, the design of new secondary stations will feature prefabricated, i 

demountable, and s semi- portable type shelters that can be — knocked 

_ Partial-Record Stations. —The principal sitet on partial- records ison 
flood peak flows and base flows, the two extremes of river regimen. Flood 
¥ peak partial- records } are frequently obtained by use of the “crest-stage gage. * 
-_ A small vertical pipe on a tree or bridge ; abutment has inlet holes near the 

bottom and is charged with a spoonful of granulated cork. After a flood the 

- cork leaves a high-water mark k on a graduated st stick left in the pipe. The engi- 
4 neer visits the station periodically, preferably after each large rise, and notes | 
the peak stage since the previous visit, removes the mark, and re-charges 
with cork. _ The cost of operating such a station is from 5% to 30% of a com- 


_ plete gaging station cost, depending on need for and difficulty of calibration 

of the site todeveloparatingcurve, 

“es Base flow measurements are made several days to a week ormore after 
_ the most recent rain on an area. The flow includes no direct storm runoff and 


is essentially all from flow is 


maximum ebovation for 1 the year, , and then measured again in the late summer 
or fall, when ground- water levels and discharges to the streams are usually | 
at or near the minimum for the year. Thus, two widely separated points on 
the correlation graph are obtained each year for several years to define the 
base ‘flow at any particular site in comparison to a nearby secondary or pri- 
- mary gaging station. The cost is about 5% to 10% of the cost of a a 
_ Streamflow Correlations .—The correlation technique used in the nationwide — 
review of the stream-gaging network was the comparison of —— 
_ monthly mean flows. A graphical solution was usually made, as shown in Fig. | 
; 5. At” least 60 months (5 yr) were plotted and a smooth mean curve drawn 
through these points. The standard error was estimated by drawing two curves 
parallel to the mean curve at suchdistances from it as to exclude the top one- 
im sixth and bottom one-sixth of the plotted points or months, thus including two. 
thirds of the points, or the part included within plus and minus one standard > 
' deviation. In this example, two-thirds of the wane were between +22% and 
In general, the correlations ‘of streamflow characteristics were found 
be good in the humid east and southeast. Correlations were good to fair = 
the Midwest and northwestern parts of the country, and fair to poor in the 


arid southwest and states immediately east of the Rocky Mountains. as 
_ Other possibilities for making correlations between primary and secondary 
; stations or -partial-record stations are available. Concurrent daily flow, 
= base flow conditions, is especially useful for low flow partial- — 


stations. Discharge for equal percent duration or equal recurrence ee 
instead pad monthly mean flows to extent records at 
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it "especially where great of the relation 
would be required, to make a study of flood frequency or drought esnenr 
by other methods. There is the additional prospect that rainfall data can be 
used to improve | the correlation between the flows of two streams. Multiple — 
_ correlation using records for several nearby gaging stations may sometimes, 
but not always, improve the results. 
Scope of the Review. w.—Insofar as it was practicable, the USGS’ review of 
the “network of gaging - stations included stations operated by states or other 
agencies such as the Mississippi River Commission, International 
Water Commission, United States and Mexico, the Corps of 
Engineers, Bureau of Reclamation, agencies of the Department of 
Agriculture _ The criteria for inclusion of stations operated by other agencies” 
states required that the station records be published regularly and that. 
“the station be suitable for inclusion in the hydrologic network. If a station 
record collected by another agency was considered accurate and opt 
available to potential users, it was felt that there would be no need for dupli- = 
cation with an additional hydrologic network station operated by the USGS. 
Stations, operated by other agencies, that “would be of use for water 1 manage- 
ment purposes only, however, were not classified. These policies ‘yielded 
. _ coverage that is somewhat less than perfect for the nation as a whole. A con- 
_ siderable number of gaging stations are operated by other agencies in con- 
nection with research or operational projects, but the basic data have not 


been published though the results may have been included in 


“and useful from the standpoint of national 
; Results of the Review. —Certain aspects of the review were quite ie 

and yielded definite results. For example, it was possible to assign a classi- 

fication of best use in the nationwide program to every existing gaging station. 
= might have been expected, even this was not timeless and there have been 
a | number of reclassifications since the initial review was completed. ‘However, 

the ‘number of reclassifications has been very small compared to the number 

_ of stations, less than 1% in 3 yr, so the classification of existing — 

may be considered fairly firm. 

_ Another firm finding was that some 700 specific basins in this country did 7 

not have a single gaging station within the area that met areal primary 

requirements. Most of these basins are in the vantage Plains and Rocky Mountain 


characteristics pate than in the humid areas. ‘The correction of this this de- 
_ The review showed that more than 1,200 existing secondary s' stations could 
; be moved to other locations after a period of operation of 5 yr to 10 yr, and 
that a a good share of these were ready for immediate relocation. It was further 
shown by brief correlation tests that about 1,150 gaging stations previously 
_ discontinued had accumulated sufficient length of good records to permit ade- 
_ quate correlation with primary stations. It was also found that some 1,250 
existing stations classified in the water management categories would qualify 
for the secondary hydrologic network (except for relocation after 10 yr). Thus, 
4 total of about 3,600 points have been or are being gaged in the secondary oS 


12 “Graphical Correlation of Gaging-Station Records,” J. K, Searcy, U.S. Geol. 
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hydrologic network even ‘though | only 1, 200° stations are classified 
as secondary stations it 

result, much less firm than the ones, was tentative con- 

clusion that about 2,500 additional short-term secondary stations would be 
fae the next ten years. This conclusion was reached by adding up the ‘| 
Sa: in the judgment of the 44 district engineers of the USGS | 
across the country, would be needed in their districts. Thus, itis strictlya J 
m matter of judgment and cannot be proved right or w lite although it is the best — 


_ The last, and perhaps most hypothetical, result of the review was the indi- — 
_ cation that a | _ total of about 10,000 active gaging stations in the United States © 
is the maximum m number that would ever be needed and justifiable at any one 
time. This figure was arrived at by adding the 700 additional primary and 
& 500 additional ‘secondary stations needed, as previously described, to the * 
existing network assuming that for every station established thereafter 
another one would be dropped. The figure of 10,000 assumes a slow but steady 
increase in water management stations and a very large increase in — 


"record | stations for reconnaissance type data. No specific date was set as a 
_ target for obtaining 10,000 stations. A straight line projection on Fig. 1 indi- - 
cates reaching this level about 1975-85, much later than the goals established a 
A side-effect, if not a result, of | the review was the strengthening of the e- 
conviction of those participating in the study that a broadening of the types of © 
data coll collected would be more valuable than simply adding more years of daily 
- discharge record at all existing gaging stations. Data on flood plain inundation, 
water use, time of travel, flood profiles, bankful stage, water temperature, _ 


_ ‘Trends Since the Review.—Fig. 6 shows the number of streamflow stations 
in each classification as of January, 1957, andJuly, 1959, Of the 700 additional ¢ 
areal primary stations the review found to be needed, about 80 had been es- 
tablished by July, 1959. Only a net of 36 additional areal secondary stations 4 
were in operation but nearly 350 new secondary stations had been established - 

— during the period. Similarly, for all classifications of streamflow stations, 
= were being collected at 660 new locations, but because of selective dis- 

‘continuation the net increase was only 115 stations,  —— 
There were very few changes in number of stage stations during this: 
a "period, 34 being established and 17 discontinued for a net increase of 17. 7 
However, the expansion in the number of partial-record stations during the 
"period was particularly encouraging. The exact number in existence in 1957 
- is not known but is believed to have been about 500 low-flow and 800 crest- _ 
_ stage stations. The count in September, 1959, shows growth to more than _ 
A 2,100 low-flow stations and about 1,950 crest-stage ati 
a At the risk of misleading the casual reader, the increases in the several 
-_ categories and classifications can be totalled, though strictly speaking each | 
is different from all others. It can be generalized that surface-water 4 


Aw 


- are being collected at nearly 3,500 new points since January, 1957; the data 


: are less detailed and require analysis for full use, but the cost in funds and 
Manpower ‘Tisen only about the equivalent of of 600 new 
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to be working out very well, not only in theory but in actual practice, there 
ewin afew remaining problems. The most immediate, if not the largest, 

of these problems is that of gaining acceptance of this type of program design 
by engineers and hydrologists throughout the country. So far there has been 
no ‘outright disagreement with the concepts or opposition to the policies from 
anyone close enough to the problem to study it through. However, the eon 
and the procedures have not been well advertised and thoroughly discussed 
the literature yet, to inform fully all who r need to about this 


has arisen ina isolated instances. These can probably be re- 
- garded as the exceptions that prove the rule. Most of these instances arise 
g when it is proposed to discontinue a secondary station which is no longer 
needed for hydrologic network purposes. Better understanding of concepts 
and program goals will gradually overcome this problem. = | ~~ 
Another fairly urgent problem appears be the acquisition of funds to 
establish “bench mark” areal primary stations found to be needed. A “bench 
_mark” station is one in a nationalora state park, a wilderness area, or other 
area ea where hydrologic change due to man’s activity does not appear likely. 

. ‘There are few existing ‘stations where the hydrology of the basin will not 
change perceptibly in the next one hundred or more years. Wherever a — 
can be found in which long-term, man-made change in the hydrologic regimen 
"seems unlikely, an areal primary gaging station should be established to de- 
tect and document long-term trends as between wet and dry cycles. If an 
average of two or three such stations could be found in each state, they would - 


regional and 1 more reliable determination of the probable frequency 
recurrence of given floods or droughts, = 
Another problem arising as a result of thistype of programming is getting 7 
all basic data assembled and published for wide use. Data from partial-record | 
: stations and other reconnaissance-type information must be either fully eval- 
data themselves ‘must be published. It is to publish data 
from crest-stage gages and base-flow stations along with continuous station - 
records in water-supply papers of the USGS. 
_ The final problem which needs “mention is that of getting unc underway ‘more | 
special studies and analyses in which all available hydrologic data are as- - 
7 -sembled into interpretive reports. These reports should contain maps, graphs, 7 
= _ and tables to make the great mass of accumulated hydrologic data more readily — 
useful to practicing engineers and other data users. In addition, shorter and > 
- less technical reports should be prepared directed toward schoo! children and 
_ general public readers. _ These should contain the more significant facts to q 
promote wise public use, control, and conservation of the nation’s water re-— : 
sources. The investment in data collection will be if — 


_ The stream- am-gaging program was reviewed from several different : angles 
and found to be generally adequate. To overcome come existing defic seinen on 
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changes were found to be desirable. Accordingly, a flexible system has been 
adopted with nearly permanent primary network of index stations and 
continuously changing secondary network of short- term stations. Both net- 
works are supplemented by a large number of partial-record stations de- Ps 

_ signed to give specific data on flood peaks or low flow. © vil a. 
_ The whole system is mutually supporting and correlations between pr primary, | a 
secondary, and partial-record stations will make data collected at each much — 
logical guide for efficient expansion of the program and, if need be, for re- 
_ductions with minimum damage. Above all, the program is designed to speed 
‘up greatly reconnaissance-type coverage on small streams and to give maxi-_ 
mum information for the manpower and funds e: expended. The growth rate, in 
terms of total stream- ~Gaging stations being operated, has slowed down as a 


of network However the growth rate, in n terms of total sig- 
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more useful and valuable than if used alone. The network ‘system provides 4 
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HYDRAULICS DIVISION | 
_ ae of the American Society of Civil Engineers 


BASIN DAMAGE AT CHIEF JOSEPH DAMa 


SYNOPSIS 
Chief Joseph Dam on the River, in the state of Washington, is 
the first dam downstream Grand Coulee Dam. The dam, approximately 
200 ft high, is a concrete gravity structure with a tainter- gate controlled } 
ogee spillway, that discharges into a hydraulic jump-type, energy ag cpa 
having a single row of baffles and an end sill. Damage of the concrete in the — * 


basin was discovered ‘in 1957, two roe after the project became operative. 


The occurr 

common event. 

er the events “causing this damage are _ attributable to construction or to 

conditions and to define the specific situations r responsible. The 

repair of a spillway energy dissipating structure » for a major dam after a 

project is placed in service is often costly. Therefore, a full understanding _ 


- —Note.—Discussion open until August 1, 1961. To extend the closing date one month, 

a written request must be filed with the Executive Secretary, ASCE, This paper is part 
of the copyrighted Journal of the Hydraulics Division, —— of the American So- 

ciety of Civil Engineers, Vol. 87, No. HY 2, March, 1961, Rae a ee 

Presented at at ‘Seattle, Washington Meeting, Hydraulic Division Conference, August" 

1 chf., Planning § Sect., ae of Engrs, » Dept. of | the Army, Seattle District, Seattle, 


* 
&g 
4 : together with an evaluation of conditions that are considered to have contri- — 7 
— 
— 
7 
— 
= 
— 


me 
vd 


~ 


— 
— 


: 


M 140d390 


Same basin, that is a conventional structure constructed in accordance with 
2 conventional procedures, should be helpful. It is interesting to note that dam-— 


major structures | on the Columbia River, 


PROJECT LOCATION SPILLWAT FEATURES 


“Chief Dam is located 51 mi miles downstream from Grand Coulee 


‘sists of a concrete gate- spanning the river chan- 
tho , and | an intake structure and powerhouse located just downstream, as 
shown on Figs. 1 and 2. Use of the spillway is required annually, generally — 

over a 5- to 7-month period during which time peak flows of about 350, 000 

to 500,000 cfs can be expected. The historical flood on the 
‘River 2 at Chief Joseph site is about 730,000 cfs. The flood of record perso | 
(1948 is 640, cfs, and spillway design discharge is 1,250,000 cfs, 


HY2 CHIEF JOSEPH DAM 
circumstances leading to the occurrences of damage in Chief 
— 
il 
A 
| 
FIG, 2,—CHIEF JOSEPH DAM, FLOW OVER COMPLETED SPILLWAY, 
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Normal head at the project during ordinary flood periods is in the order of 
magnitude of 150° ft from pool to tailwater | and 200 ft from pool to ee : 
_, Flow over the spillway is controlled by nineteen tainter gates, 44.5 ft high 
and 40° ft wide. During flood periods, the pool will ordinarily be maintained | 
a: Ele, 946 to 948 with a maximum rise to Ele. 957 for the spillway design : 
a. , A paved concrete stilling basin that inducesa hydraulic jump for energy 
dissipation is provided at the toe of the spillway. ‘The stilling basin is con- 
nected to the ogee by a a 100-ft-radius bucket section and is 167 ft long, from — 
point of tangency of the bucket to the end sill, as shown on Fig. 3. A single ~ 
a baffles and. a stepped end sill are located in the downstream portion of the 
stilling basin. - The stilling basin details and general plan are shown on Figs. 
y and 5. The stilling basin floor is comprised of slabs designated in Fig. 5 5 


from Row | D upstream to Row Al downstream. 


‘TABLE 1 _-HYDRAULIC JUMP DE ‘PTHS 


_ Flow desi-| FI Q per feet 


ix wh (2) 


3-Year 


a 12, 0 114, 


Frequency 


(Dy = Theoretical depth at at toe of jump. 
= Theoretical velocity at of jump, 


‘Theoretical | required t tailwater depth to form jump. 


= = Actual tailwater depth navailable. 


‘The slabs have a of 5 ft and are tied to the 
foundation with an extensive system of grouted anchor bars. _ Only the down- | 
stream slab on which the baffles and end sills are located is reinforced; this 


slab is heavily reinforced to take the impact of high velocity flows on baffles © 


streambed, An extensively jointed granitic rock underlies the basin and is 
found at or near the “surface of the streambed downstream of the stilling 
"basin, The rock jointing varies: greatly in strength and in many areas the 
is subject to erosion ender high velocity flow. 
The stilling basin causes energy dissipation by hydraulic jump action, A 
-.., row of baffle piers, : streamlined on the sides to reduce cavitation po- “a 
tential and a stepped end sill cause the jump to occur on the basin under max- _ 
imum flow conditions and also aid materially in damping downstream wave 
"action, . The jump conditions in | the basin for spillway design, the historical | 
maximum flood and a 400,000-cfs flow (3-yr frequency), are given in Table 


ro - and end sill. “The floor of the basin is level at Ele. 743, about 12 ft ie 


“4 
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— 
; a 
4 
91.0 | 090 
91.0 
400,000} 440. | 39 | 113, | 53.7 | 095 
4 
i 
q 
s of 1.0 for flows 


— as of 1960. The: design features of the stilling basin were thoroug 
ly tested in model studies, using a 1:33 scale, 3-gate sectional model, and —_ 
vf 1:80 scale general model that included the entire spillway. ete my 
a One of the outstanding features of the design is the coordinated effect of 
_ the baffle piers and end sill in causing formation of the jump. The baffle piers 
a were located close enough to the end sill so that the end sill aids increatinga 
in of the baffles that the occurrence of 
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4, BASIN EROSION co 


pressures. — basin wa: was carefully tested in a 1: in a 1:36 sectional r model 
to obtain an adequate length | and depth and to avoid cavitation pressures on nO 
baffles and end sills. The model tests indicated an adequate margin of ale 
for satisfactory stilling basin performance for all flows» and showed no in- 
‘dication of cavitation pressures for flows” of less than 900,000 sec-ft. 
significant feature disclosed by the model tests was that a slight variation in 

the curve of the sides of the baffle piers from the theoretical shape ae B.! 
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CHIEF JOSEPH DAM 


resulted in n high at design flow not when the 


first diverted a of stilling basin in 
1982. until the was placed in service in 1985, a 
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‘FIG, 4. CONTINUED 


of t the stilling basin was utilized continuously each year ‘under construction 
conditions. Continuous use of the spillway lasted until early in 1957 when . 
sufficient hydraulic capacity was installed in the ‘powerhouse to permit com-— 
plete diversion of low aie a . Thus, initial inspection of the basin was 
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_ ; primarily of four underwater traverses by a diver from one side of the basin § 

“a 


4 ‘sill ogee This "inspection 1 revealed several 
gE in which erosion had occurred to —_ of zero to 12 in.; several ex- - 
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areas to expose the 6-i aggregate used in 
= paving. Although no evidence of rock from construction operations was found 
4 the basin. There were numerous loose pieces of aggregate in the basin - 


that had been rounded and reduced in size to 2-in. or 3-in.-diameter spheres | 7 
. _by washing action. Several loose pieces of broken 1-1/2-in. square reinforc- 


ing bai bars w were — found in the basin. — or ten baffles were reported as 


mper O1 
and marked damage to a nu ‘orl 7 
;_ = ere was exposed reiniorcing stee ill. The surface of the slab was eroded oo 
d several sections of end sill. 

4 
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RESULTS OF SURVEY OF NOVEMBER 1957 

29 A detailed underwater survey of the basin ‘was undertaken in November, 
1957. . A team team of two divers: fabricated measured 


Left training wal! 
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elevations over the floor r of the basin on a grid spaced at approximately 10-ft 
i centers and also made detailed measurements of baffle piers and end sill 
The results. of this survey are presented in detail in Fig. 4 and condensed 1 
to. show overall major damage areas on Fig. 5. Fig. 4 indicates the erosion . ¥g 


contours | with the amount ¢ of erosion indicated in inches. | The symbols indi- 


CHIEF JOSEPH DAM 
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cated in Fig. 4 are. defined as follows: SF indicated a sand finish, with no 


erosion of concrete; is pea gravel in relief; AR defines aggre- 


9 ® 


gate not eroded: T.S, denotes temporary sluice, is (1960) filled 


o~ concrete, Scarcely a single slab within the basin failed to show some 
degree of erosion. Because the aggregate used in paving of the slab was * 


- PUAN OF STILLING BASIN SHOWING MAJOR 
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CHIEF JOSEPH DAM 
"graded up to a a6: in. maximum size, any erosion | of the slab concrete resuite 
= development of an an extremely rough surface oF Sass in n Fig. 12, and $ 
“4 as reported by divers. Major damage are 
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LEGEND FORSTILLING BASIN 


= Fiat slab or end sill erosion from 12" 
and erosion exceeding 2" in bucket section 
of ogee and on D series sieb 


Slab or end sill erosion exceeding (2 


—<— Baffle piers with erosion of 6-12" 


— 
Baffle piers with erosion exceeding 


EROSION, 1957 SURVEY 
‘erosion exceeding 12 in. and more lightly eroded areas with erosion of less 
than 12 ‘in. are ind indicated on 1 Fig. 5 5 by \ various types of hatching. 
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~ sila 5 ‘shows that only in the “A” row of slabs did erosion exceed 12 in. *, 
The most severe damage appears concentrated inlocalizedareas immediately 
in front of baffle piers 37 to 26, and 3 to 9, between these baffles and 1 between ve 
” baffle. piers and end sill. The areas of greatest slab damage appear to 


LEGEND 


> Total erosion, March 1960 Increase in erosion, December 1957 to March 1960 © 


TAv Average erosion in slab. ,in feet and hundredths. Tay Average increase in erosion in slab, in inches. 
TMi Minimum erosion in slab , in feet and hundredths. -IMi Minimum increase in erosion in slab, in inches. s 
Maximum erosion in slab, in feet and hundredths. —TMa- Maximum increase in erosion in slab, ‘in inches. | 
a No erosion measurement taken tied 1957 survey; 
visual inspection only was 
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‘severe damage, Two small scour ‘holes saorntecsalit 5 ft deep were found | 
_ downstream of baffle piers 22 and 23, and 30 and 31; other areas in which 
erosion occurred were generally less” than 2 ft in depth. In the 
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of slabs, that the row with reinforcement, the surface layer of of 

2 steel is extensively exposed in locations at which exceeds 12 in. 5 aera 
; _ Detailed measurements were made and plotted of damage to each baffle. 
Of we forty-t three baffles, ‘it was found that twenty erosion of 6 in. 


‘sie 


Note: The overoge toilwoter elevation during the stilling bosin survey 
was 770 Yeet, U.S.C. BGS, 
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1960 


on and could be considered in satisfactory condition; ‘seven showed erosion 
of 6 in. to 12 in., indicating near exposure of reinforcing steel; and sixteen : - 
showed erosion exceeding 12 in. that exposed reinforcing steel. The last two 4 
- a of baffles are indicated by shading on Fig. 5. Damages to baffles 
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——~ TABLE 2,—PAVING EROSION DATA CHIEF JOSEPH DAM STILLING BASIN 

| ‘Total erosion March 19602 
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Averages 


Second row of slabs upstream from baffles © 


4444446 


LOR 


b Visual inspection only, 1957 (except several isolated measurements). 7 oe 
© Visual inspection made in 1957; no measurements were taken, ite 


602 
RA Ab feet inches inches | inches 
ih p4) | 1-3/4 | oO 
(b) First row of slabs upstream from baffles 
0 1-1/2 | oO | 3-1/2 
(1-1/2 | «41/2 | 8 
= | in 
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JOSEPH DA 111 
on the front face and sides with generally accom= 
Pa _ panying damage to adjacent slab areas. However, of the solid shaded baffles, 
_ numbers 4, 26, 27, and 28 suffered damage principally ; to the tai tail of the e baffle, — 
number a large piece knocked out of an corner. 


ae To establish the rate of progression of erosion damage and to determine 7 
erosion might be attributable to operating conditions, a secon under- 
water survey was undertaken in March, 1960. The results of this survey 
_ tabulated in terms of total average, minimum and maximum erosion in feet of | ib 

= each slab in 1960 and the corresponding increase in erosion in inches since © ; 

: 1957 are shown in Fig. 6 and in Table 2. In the two flood seasons between 

: surveys, peack flows of 377,000 in 1958, and 395,000 in 1959, were experienced, * 
_ both of which are less than 4- -yr frequency occurrences, It was not possible 
_to make measurements in the “D” row of slabs in 1960 because of the hazard 
of ice falling off the tainter gate assemblies. During the survey, the divers %) 
were able to distinguish many light colored fresh pockets throughout the — 
basin: in which aggregate had been recently eroded from the surface of the _ 
concrete. Shere was extensive occurrence of fresh impact marks on 
baffles and end sill, probably caused by small rounded balls of aggregate 
_, The average increase in erosion between the 1957 and 1960 survey in each — 
a? of rows “A,” “B” and “C” was found to be approximately 1-1/2 in., although © 
_ the | total erosion in both “A” and “B” rows of slabs appreciably exceeded that a 


| in upstream row “C,” ” Correspondingly, the 2-yr increase in damage relative | 
. 4 to the initial damage | measured in 1957 is much greater in the upstream row 
of slabs. Additional surveys in the future ‘willbe needed to determine whether w 


this trend will continue. 


In the 1960 survey, detailed measurements were made of only fourteen 


baffle piers: including ones representative | of light, moderate and heavy dam- 
in 1957. In general, baffle piers” with surfaces that were undamaged 
in sound condition in 1957 were also found to be in sound condition in 1960, 
"Increases in baffle pier erosion was limited to surfaces that were found to — 
have been markedly eroded in 1957, 
Changes in end sill erosion followed the same - condition n pattern as as the = 
baffle piers with added erosion being generally limited to some areas that | 
previously been damaged. The overall average in 
damaged end = areas was only 3/8 in. 


CAUSE OF DAMAGE 

The | oo observed in the last 2 yr of spillway service appears minor 
as compared to the pattern: of major damage found in the 1957 survey. Tt 
|e ‘that most of the damage | to the stilling basin has originated from 
conditions during project construction. The pattern of localized, deep and 

' _ extensive erosion between the end sill and certain baffles and the damage to 
“front faces and sides of these baffles indicates that much ofthe damage must > 

- have resulted from scouring of foreign materials in the basin, particularly, — 
there has ne not n any ‘development of new areas Gamage 
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CHIEF JOSEPH DAM 
1957, An of conditions during project construction has verified 
= on several occasions, rock was introduced into the stilling basin ‘ion 
river diversion, and this situation combined with diversion 


Sa construction of these cells it became necessary to dump a sub- 
stantial grate of 1-ton to 3- -ton rock immediately upstream of cells 3 through 


"struction of downstream cells. This rock groin, shown on Fig. 8, was only 
partially removed. Five low monoliths approximately at streambed level were 
‘constructed in the first stage cofferdam together with four pairs of 8-ft- wide 
= 16 ft-high temporary sluices with invert elevations 14 ft above riverbed 
for use in passage of second stage diversion flows. ee) ) eee 
‘First stage construction was completed in January, 1952; the cells were 
removed and second stage diversion, consisting of cells numbered 26 through 
_ 46, were installed and operative by April, 1952 so that the 1952 high-water _ 
was entirely diverted through the low monoliths and temporary sluices of the ; 
first stage construction on the : right side of the channel, Fig. 9 shows flow 
conditions in March, 1952, with the second stage cofferdam nearly completed 
just prior to high water. Note the white water over the remains of the rock : 
_ groin upstream of the first stage construction and also the outside rock groin | 
being placed at the _ upper end of the second stage diversion cofferdam. In 
addition to rock wash from these sources, a near failure of cell 29 in the 
second stage cofferdam during the 1952 high- water required emergency 
- dumping of rock around the outside of the cell. The high velocity flow around 
the outside of the cofferdam undoubtedly washed a considerable amount of 
; rock into the first stage stilling» basin. There was also rock wash into the 
basin placement of armor ‘toe rock in the wraparound abutments, 
particularly | on the right bank side. With this set of circumstances, the 
"pattern of concentrated damage downstream m from the etal block monoliths 
8, 10, 12, and 14 becomes understandable. 
"Oh the third stage of diversion in 1953, Fig. 10, the cellular cofferdams 
_had been removed and flow was diverted through the twenty-four temporary 
the low monoliths that were still: at streambed level. During 
“flood periods a very heavy concentration of flows resulted through the low 
-monolith diversion notch that, in turn, generated a strong eddy currentin _ 
_ the basin. This same pattern of flow was repeated in 1954, except that with © 
the raising of the diversion notch was now 100 ft above the 


ersion 1 notch and the location of eddy currents, 


illustrated by Fig. 13, shows exposed bars and erosion 
on the face of sill the end sill concrete eroded 
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construction for river diversion. In the first stage of diversion a cofferdam 
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FIG. 8,—FIRST STAGE COFFERDAM CONSTRUCTION SHOWING ROCK GROIN 


“FIG. 9,- 


~SECOND STAGE COPFERD AM PRIOR 1 20 OCCURRENCE OF HIGH WATE - 
AND WHITE WATER ROCK GROIN AN D ERN 
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_ SHOWING EDDY CURRENTS = 
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il. —FOURTH STAGE CONSTRUCTION DURING NG HIGH p 
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1961 
initial underwater inspection made in 1957, two years after the 
_— spillway was placed in full service, showed no rock in the basin other than 
small amount of dislodged aggregate, ‘This confirmed model tests that 


showed ‘the completed basin in normal operation | to be self-cleaning for 
discharges in excess of 300,000 cfs. However, flow conditions throughout the : 


that rock washed into the basin was ‘nie asa a tumbling agent to: scour 
and abrade concrete surfaces. 

pe Another damage factor during diversion was that the very heavy concen- 


7 tration of flows in the Sourth win ad diversion notch approached maximum 


4 
— 
S diversion period_were so differ from normal operation that it is probable — 
— 
- 
y ciency than that for which the basin was designed, as indicated in Table 3. : Lt Ras 


Bz of adequate tailwater depth | created a a ‘temporary, condition for 


— possibility of "stilling basin damage because of the quality of the e con- 
crete has been considered. Concrete for 


the stilling basin was designed 


FIG, 13 13, _—UNDERWATER PHOTOGRAPH SHOWING REINFORCING 
_ BARS IN END SILL BENT UPSTREAM, 


compressive strength of 3,000 psi. However tabulation of concrete 

from construction records (Table 4) indicates, this strength 

quirement was substantially exceeded in many areas. Slab and unformed — 

surfaces of the _top of baffles, end sill and ogee bucket section received a — 

steel trowel finish. Smooth- -surfaced tongue and groove lumber was used 

for formed sections of the baffles and end sill. Concrete was water- cured» 

for 14 days. A la prnpection of concrete in the north half of t the stilling © 
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was “specified. Three possible “an of damage | to concrete other than rock © 


rasion are believed possible. 
__ Erosion of Concrete by Water-Jet ‘Action. — The exposed aggregate condi- 


3.— 
‘DIVERSION ENERGY DISSIPATION CHARACTERISTICS 


city 
feet per second | 


000 
000 
— 44 
days) 
300,000 | 


66 


Compressive ‘Strength? a 


TABLE 4, —PROPERTIES or ¢ USED 


Maximum} Cement ater Entrained| 

90 days” 


Size Ag-| Content 7 

in inches |per cu yd| weight | 
. 5, 600 


cylinder tests 
= front face of many baffles could wns be attributed toe erosion caused by 
_ high velocity water-jet action. High velocity flows through the stilling basin © 
i during diversion have ranged from 30 fps to over 100 fps. Many of the changes” 
E in the direction of flows were induced by impact, either from flow over mono- 


a liths, through temporary Sluices striking the ogee or stilling basin floor, or 
from diffusion of flows by 


basin was placed from freezing pro- 
H STAGE 
| 
tical depth of tailwater for jump formation. — 
| depth of tailwater for jump formation, 
cation 
4 (1) 
4 
4 | 


rose 
first erosion due to. “water- jet action was reported during the initial 
diversion of high water in 1952 when water flowed up to 15 ft deep over mono- .. _ 
liths 7, 9, 11, and 13. Inspection immediately after the flood-flow period re- 
vealed light scour and an exposed aggregate condition of concrete on these a 
monoliths below Ele. 775. Extensive tests conducted by the United States 
Bureau of Reclamation Dept. of Interior indicates that concrete will stand 
the flow of clear water at high velocities without appreciable damage, provided e ; 
that the flow is smooth and the direction and velocity of the water is not — 
abruptly changed. However, water-jet tests at velocities from 100 fps to a. 
175 fps, showed pitting and surface roughening when the jet angle « ‘exceeded 
5°. Damage from scour at moderate to high velocities is apparently a random 
occurrence depending on the density and quality of the surface attained in i 
construction. It is expected that resistance to high velocity scour varied 
: throughout the basin as no special quality of finish was specified. a a 
7 _ Exposure of Under-Strength Concrete to High Velocity Flows.—Some con- 
° _ crete in the north half of the stilling basin was exposed to high velocity river | 
ranging to 40 fps within 60 days after placement. Concrete strength 
at this time may have been below normal because it was placed and cured in — 
freezing and subzero weather during November 1951 and December i951. ae: 
_ Some concrete was exposed 45 F and colder, high velocity 
flows less than one month after placement. 


Concrete Surfaces. — Paving of the slab using concrete with 6- in maximum 


the “slab surface. The “exposure of this aggregate created a an n immediate 
- cavitation hazard that with cavitation damage could spread progressively to 
undamaged | adjacent areas. The model tests showed incipient cavitation dam- an 


age for minor ‘deviations from the | theoretical curvature on the sides of the 7 7 


- baffles even when using a smooth curve. In the prototype, side curves were 4 
formed with narrow tongue and groove lumber, resulting it ina marked devi- 
ation from ‘the type of surface tested in the model. It is ‘apparent that any 
irregularities in the surface resulting from forming, finishing, or patching | 
would greatly increase the possibility of cavitation. In construction, a 
baffles were formed only on the sides and front with the top and back et : 
placed and finished without forms. The use of this open-type forming limited 
the amount of vibration — and compaction of the concrete because of sloughing 4 
of the back slope, , thus creating | conditions that caused rock pockets on formed — 
areas. Although all such areas” were patched, these patches are areas of r 
weakness andcould havecreated irregularities or roughness to start cavitation. PY 


REMEDIAL MEASURES TO PREVENT OR MINIMIZE 
DAMAGE IN FUTURE CONSTRUCTION 
primary cause of stilling basin damage at Chief Joseph ts 
. be scour caused by large size rock particles introduced into the basin by 
_ diversion flow conditions during construction. It is also possible that damage 
could have been reduced or minimized by some higher quality control of ll 


W AU contributing factor to the occurrence of this damage was undoubtedly 


and field levels, Full consideration of the problem during the model tontien % | 
stage led tc to in design. For example, the tendency of the 
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-stepped-end sill to retain rock — sloping t the 

. of the sill. It might have been possible to raise the elevation of the low —_ 

' — liths 10 ft to 15 ft above streambed to prevent rock from entering the i 
during diversion, but this would have necessitated higher cofferdams and the — _ 

single cell- -type of cofferdam used was already at maximum height for the 


strength of steel sheet pile that was available. Thus, a double cell- type. coffer- 


dam probably would have been required. Le 


Pee full consideration of the problem might have indicated a different plan 
of 


diversion during construction; however, the costs might well have exceeded — 
costs of repairing damage to the basin. Greater restrictions placed on the 

_ Operations of the contractor to minimize operations that introduced rock into — 
the stilling basin would also have increased the costs. A higher quality con- 

- ‘erete surface would have been desirable and steps were taken several years — 
ago by the Corps of Engineers to specify superior quality concrete of higher 7 
_ strength and better finish and employing small size aggregate for use in 
‘The steps taken to control damage during construction m must necessarily — 
be tempered by balancing ultimate costs of repair against the added costs» 


CONCLUSIONS 


iy Damage t to the Chief . Joseph Dam s stilling basin re reported herein was caused 
primarily by rock entering the basin during the diversion phase of construc- 
waren age p rs 7 


“ tion. A greater awareness of this type of damage potential by both designer 


field forces is necessary to avoid ft ie recurrences. 


data used herein are taken from records of the Seattle District, 


Corps of Engineers. Work on - stilling basin studies has been under the di- 4 
rection of RP P. Young, Colonel, U. S. Army and District Engineer; and Noble . 
A, Bosley, | Chief, Engineering Division. All of the detailed sur- 


7 
under the technical supervision of W. W. Whipple of this office. | 
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FORECASTING RIVER RUNOFF BY COASTAL FLOW INDEX 


An index method for forecasting April through September runoff for the 
Columbia River at The Dalles, Oreg., is presented, based on the ee ; 
_ between winter runoff in coastal drainage areas and spring snowmelt runoff 


of the Columbia Reliability: of is comparable with those 


| g The Columbia River is one of the principal snowmelt runoff streams in the 
- United States. Inasmuch as the flow of the Columbia River is primarily from 
iz _- snowmelt, an evaluation of existing snowpack water equivalent makes it possi- 
ble to forecast the total runoff potential several months in as 
forecasts enhance the use of multiple- -purpose reservoir storage for obtaining 
optimum benefit of water resources for all water uses. The evaluation of the 
total: runoff potential of the Columbia River is necessary for coordinated 
q 7 system operation of reservoir storage for flood control and other water use 
requirements, as well as for flood potential evaluation prior to completion of 
a an adequate system of reservoirs for flood control. The Corps of Engineers, | 
Leads United States Army, is responsible fc for flood flow regulation in the Columbia 
‘River Basin. The prime ‘basin forecast of the total runoff is for the gaging 
station oda at The Dalles, ‘Oreg. (D.A. = 237,000 sq miles), at which wil 


-Note.- Discussion ‘open 1 until “August 1, 1961, To extend the closing ¢ date on one e month, 
a written request must be filed with the Executive Secretary, ASCE. This paper is part | 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 87, No. HY 2, March, 1961, eT te 
Hydr, Engr., U. S, Army Engr. Div., North Pacific Div., Portland, Oreg. 

= Meteorologist, U.S. — Engr. D Dist., Portland Div. , Portland, Oreg, 
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the runoff the entire iver east oft the 


, -Methods of forecasting seasonal runoff volume from snowmelt have been 

based primarily on use of precipitation or snow accumulation measurements, 

The most commonly used procedures employ the concept of indexes, whereby 7 q 
point ‘measurements of hydrologic variables are correlated with runoff 
various statistical techniques. Coefficients so derived from historical data _ 
: _ provide the basis for forecasting future runoff, by applying them to observed _ 
a index values. Forecasting methods of this type depend upon an adequate period 


of record for statistical treatment, and these index procedures may involve 


States Dept. of Commerce (USWB), Cooperative ‘Snow Investigations, ‘was ——- 
of practicale and reliable methods of forecasting 
streamflow resulting from snowmelt. The emphasis in the analytical work | 
that was carried on by the Corps of Engineers’ Snow Investigations Unit in | 
"pursuit of this objective was in the rational evaluation of the water balance on 7 
basin areas involving snow, and the determination of the storage of water in 


snowpack in relation to the other factors involved in the 
of these investigations are available. 3, 


are useful in deriving practical forecasting equations. General 
; ‘peetoie of the hydrologic cycle and causative meteorologic factors involved 
7 in the deposition of precipitation, as appliedtoa specified basin, are essential — 
for the intelligent selection of suitable parameters fi for forecasting 


-Funoff by index methods; 


__ Due to the fact that the Columbia River Basin is located in the witiee = 
latitudes, in the zone of prevailing westerly air flow, the sources of 
moisture ‘supply are the maritime air masses of Pacific origin, particularly 
, during the winter season. Precipitation occurs with the eastward movement a 
a large-scale frontal systems from the Pacific. These systems, which are © 
usually | occluded by the time. they reach the coast, “produce varying amounts 
of frontal precipitation. Additional precipitation results from orographic 7 = 
Fe Fr the moist air masses as they pass over the Coast, Cascade, and 
Mountain topographical barriers. Inspection of upper-air ‘charts® shows 
= during times of significant winter precipitation over the Columbia River 
Basin, the average air flow at the 500 and 700-mb level is generally from the 
west, but “May vary from southwest to northwest. Further inspection reveals 
= that the basin is generally well centered in the westerly stream of air which, 
on the average, covers a width of about 15° to 20° of latitude. = © 
_ The initial orographic uplift of the air stream over the Coast and cascade 7 
Ranges results in deposition of a portion of the total precipitable water of the 
i air mass before it enters the Columbia River Basin, east of the Cascades, i 
_ The lower elevations and predominant maritime air-mass characteristics of a5 
the drainage areas in the Coast and Cascade Ranges result in higher tempera- . 


2 3 “Summary Report of the Snow Investigations, Snow Hydrology,” U.S. Corps of Engrs. > ic. 
North Pacific Div., Portland, Ore., June 30, 1956, 

4 “Application of Snow Hydrology to the Columbia Basin,” by iineiediiaie and Peter ea 


Boyer, Proceedings, ASCE, Vol. 85, No. HY 1, January, 1959, 
Weather Bur., daily series synoptic weather maps, northern hemisphere sea 
level and 500 millibar charts, Dept. of Commerce, Washington, D, C. 
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are in the higher and more continental areas of the Columbia Basin nearly all of + 
precipitation falls as snow, 
‘The winter runoff from the Coastal streams is considered to be a reliable — 

index of the total winter precipitation falling on the Columbia River Basin, | 
_ because it samples the total incoming moisture supply i in the air masses dur- 
F by times of strong zonal flow. Such an index applied to the Columbia Basin 
a as a whole would be a measure of the total moisture supply to the basin and i 
mi therefore may be used as a forecasting parameter. However, it would = 
i - evaluate the distribution of precipitation on the various tributary streams or 
; ve sub-basins by the application presented herein for the Columbia River Basin a 
. a _ With this concept in mind, a project was carried on in 1954 as part of the ’ 
1 work of the Snow Investigations Unit, whereby forecasting relationships were * 
men for the Columbia River at The Dalles, using graphical correlation _ 
techniques . The period of record used in those studies was the 22- yr period a 
. of 1929 through 1950. _A report on the study has been issued. 6 Subsequent to = 
é a <= that time, nine additional years of record have become available for study i 
3 (1951 through 1959). Furthermore, the North Pacific Division office acquired» 
i. use of an electronic digital computer in 1956, which could be utilized for 
—-miultiple- regression analysis. Accordingly, the analysis presented herein has | 
_ taken advantage of these conditions and is independent of that data previously - 


= 
set of correlations of runoff with a coastal winter-flow index and 
significant runoff ‘parameters, by which ‘runoff o of the > Columbia | River at The 
Dalles for the p period “April through “September may be determined will be 
_ presented. Both winter and springtime effects on runoff are evaluated, and 
_ results derived by arithmetic statistical methods are presented in eraphical 
form, for the evaluation of effectiveness of the parameters in estimating _ 
runoff, and the range of values normally experienced. The accuracy of fore- a 
casts is briefly compared with other methods of forecasting in common use — 
the Columbia Basin. Sy 
Considerable use is made of statistical techniques in deriving the relation-_ 
ships presented herein. It is assumed that the reader is acquainted with basic 
| fo gener terminology and analytical procedures, and for the sake of brevity, 
no general summary of statistical methods is included. For those who oe 


_ information of sian init textbooks on statistical analysis are available. 


Columbia River Near The Dalles, Oreg.— — The United Sts States Geological ‘Sur 
vey, of the has continuous records of stream-— 


— 


oat - 6A Coastal Flow Index Method of Forecasting Seasonal Runoff for Columbia River 
near The Dalles, Oregon,” by David M. Rockwood, Snow Investigations Research Note 
Note 2 No, 23, U. S. Army Engr. Div., North Pacific Div., Portland, Oreg., September 30, © 
7 «Statistical Methods,” by Herbert Arkin and Raymond R. Colton, Outline 
Series), Fourth Ed., Barnes and Noble, Inc.,, N. Y. x 
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flow for the River near The Dalles, June, 1878. The 
> record for that station is the longest in the Pacific Northwest and is rated 

“excellent” ” by the Geological Survey. It is the furthest downstream of any 

rated gaging station on the main stem of the Columbia River, and the drainage 
- area is 237,000 sq miles. Annual and April-through-September runoff volumes © a 
vere determined from this record, Appropriate corrections for major reser- 7 _ 
storages, including Kootenay, Hungry Horse, Flathead, ‘Oreille, 

: Franklin D. Roosevelt, Coeur d’Alene, Chelan, and Brownlee, were applied to | 

Observed flows. There are many irrigation diversions and small reservoir 

‘ storages” which affect the flows at The Dalles, but their net effect on annual 
and seasonal runoff is small in comparison with the total volume. This a 

particularly true with regard to year- -to-year differences. . Their effects were 
i Coastal Winter- Flow Index.—It has been pointed out that the 1 use of a . coastal 
winter-flow index is predicated on evaluation of runoff from areas which are 
at low. elevations relative to the Columbia Basin as a whole, and also which 
- - are well located with respect to the air stream during average winter air-flow — 
‘ conditions. In addition, records of sufficient length for analysis of data are , 
"required. The 31-yr period 1929 through 1959 was chosen for study, principally — 
5 the basis of continuity of streamflow records for coastal streams. The 


combination of streams selected for the coastal winter-flow index is as fol- eg 


Ga ing Station Miles 


illamette River at Albany, Oreg. 
Wilson River near Tillamook, Oreg. 
“Lewis River at Ariel, Wash, 
Chehalis River near Grand Mound, ‘Wash. 
Fig. 1 is a map of the Columbia River Basin and the coastal sill ‘showing fy 


End Willamette miver at Albany was selected because of its long period of = 


as the coast range. Due to. its moderately high < average winter temperatures, 


: ‘oo 0% of the annual runoff occurs during the period October P through 


period through 1950, ‘no correction incorporated for reservoir storage 
because of the relatively small amount of reservoir capacity above the station 


_ March. Records for the station are termed ed “good” by the U USGS. For the 
— 140, 000 a ac- ft of storage out of 10, 000, 000 : ac- ft annual il runoff). For the 


‘Dorena, and Lookout Point Reservoirs were applied to observed 
Wilson near Tillamook samples the air flow north of the Willamette 
‘River Albany, and it represents runoff from the west side of the coast 
* range. Only 6% of its area is above 3,000 ft and about 80% of its annual runoff b, 
occurs from | October through March. ‘records a: are termed “good, ” and 
- there are no upstream diversions or ‘regulation, Lewis River at Ariel was 7 
used for sampling a west slope Cascade drainage area in Washington. Its — 
elevation is about 2, 800 ft, somewhat higher than the other ‘index 
_ drainages, and about 62% of its annual runoff occurs from October through 
_ March. The records are termed ‘ “good, ” and appropriate corrections were 
eh for changes in reservoir ee A ened elevation area in Washington 


ih 
4 
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— 
= 


is represented by Grand Mound. Here only of 
area lies above 3,000 ft, oh aie of the annual runoff occurs from October | 


a _ Table 1 ‘shows average monthly runoff for the period October through | wa 
4 


March, during the period 1929 through 1959, for each of the four coastal 
winter- -flow index stations. Table 2 shows the total October-through-March — - 
flow for each station, expressed in percentage of the 22-yr average for the 4 
"period through 1950. The four station average is indicated in column 6, 
- Table 2. A minor recession correction, for the purpose of correcting the ob- , 
served runoff for water temporarily in storage in the coastal flow index y 3 
streams at the end of March, is shown in column 7. This correction w was de- 7 
termined by r recession analysis for Willamette River at Albany, in a manner 
s similar to that shown elsewhere. 8 The coastal winter-flow index (CFI) values - 
shown in column | 8, Table 2, are e those used in the multiple linear regression a 
Columbia Winter Runoff Index. —Nearly all of the winter pre- 


Cumbia River ‘Basin falls as snow. Sicensation during the winter is largely © 
base flow recession, augmented by occasional low-level rain and snowmelt — 
runoff. _ Only 28% of the annual runoff occurs in the period October a 
_ March, Nevertheless, significant variation in winter runoff does occur as -— 
result of unusual weather conditions. Inasmuch as the basin moisture input 
indexes were correlated directly with April through September runoff for the 
_ Columbia River at The Dalles, it was necessary to account for the ee j 
= of winter runoff. The “October through March runoff for the Columbia River ; 
The Dalles. corrected for storage in Kootenay, Hungry Horse, Flathead, 
_ Pend Oreille, Coeur d’Alene, Franklin : & Roosevelt, Chelan, and Brownlee © 
 eneienes is tabulated in Table 3. The winter runoff correction index is the 
eee of the October through March runoff and the costal winter-flow index 
(CFI); the values are tabulated in Table 3. 
Winter Temperature Correction Index.--In order to account for differences | 2 
‘gg evapotranspiration losses, a winter temperature index is includedinthe | 
ake The index used | in this procedure isa 5- station mean for the months 


of October , November, ‘February, and March. Climatological records for 
Revelstoke, British Columbia; Missoula, , Mont.: Spokane, _Wash.; Yakima, 
Wash.; Boise, Idaho, were used for this purpose, and values of mean 
monthly temperature are tabulated in Table 4. The temperature correction 
index is the product of the 5- .station, 4-month mean temperature, times the 


are shown on Fig. 1. 
ow Coastal Flow Correction Index. —April 1 of eachyear is a convenient - 
point of separation between winter and spring effects of runoff in the Columbia _ 
River Basin. Springtime effects are considered in this study to be confined — 
_ to the months of April, May, and June. Precipitation during April is normally 
A associated with the zonal atmospheric flow pattern characteristic of the winter - 
ss season. Accordingly, the index of moisture supply during April was the coastal- ry 
4 flow runoff. In this case it is only for the Chehalis River near Grand Mound ps 
aa Only the one station was used because the index streams of the Willa- 7 a 
, = and Lewis Rivers may have significant amounts of snowmelt runoff 4 


“Hydrology ASCE, Manuals of of Engrg. ‘No. 28, New York, , 1949, 
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TABLE 1,—MONTHLY RUNOFF IN 1,000 ACRE-FT@ 


Water Year> October through | Januar February 


548 
2 
1,582 

394 


Li 


2,093 1 


Accumulated 


E _ Average for 
Period End- 


1959 3300 


— 


124° 


— 
| 
— 
2,548 | 2195 
rw — 
— 
_s 


‘(RIVER RUNOFF 

TABLE 1.-CONTINUED 


WILSON RIVER NEAR TILLAMOOK, OREG, D. A. = 159Sq Miles* 


| 


| 


— 


5 


Average for 
= 


Period Ending 


— 
WASH, D. A, = 731 Sq. Mi 
(c) LEWIS RIVER AT ARIEL, 
Pr 4 926 | 314 194 
4 1948 — 363 . 410 
4 > 1,569 (3360 337 2 
1945 251 = 


377 17 
Accumulated 


Average for 
Period Ending 


214 a 


162 


077 
735 


1,006 


572 
242 
290 
711 
405 
405 
066 


"259 
4 


1g 


186 


Average for 
Period Ending 


9 1 


Data from U. S, Geol. Survey Water Supply Papers, _1929- 1956 unpublished records 
U.S, Geol. Survey, 1957-1959 
b Stream flow corrected for reservoir storage 
a © Estimated by correlation with Silete River at Silete, Oreg, _ ; 


ABLE 
—_— 
ll 
7 1,382] 1,658 2,089 
(d) CHEHALIS RIVER NEAR GRAND MOUND, WASH, D. A. = Miles 
19 422 | 368 
882 
ag 296 
19 1730 234 j =. 
— — ‘| 433 | 486 261 
1983 1,045 508 426 
™ 


RIVER RUNOFF 
uring the month 1 of April, while tt the Chehalis. River is is virtually | unaffected by 
snowmelt runoff during April. Table 5 lists the April runoff for Chehalis _ 
River near Grand Mound, together with the recession correction applied at 
> beginning and end of the month to account for changes in transitory storage am 
These corrections were derived from the recession analysis. The “generated — a 
runoff” (Table 5) represents the actual streamflow generated from precipi- 
TABLE -COASTAL WINTER-FLOW INDEX, OCTOBER THROUGH MARCH, 
IN PERCENT OF BASE PERIOD AVERAGE, 1929- 


Willamette} Wilson is is | 4- 
River at |Rivernear ‘Station 

Tillamook 


@ 
on 


ad 
for) 
= 


(+ +e + 1 


lend 


> 


oo 
oy 


.0 
0 
0 
5.0 
0.0 
9,0 
3.0 
9.0 
35.0 
4,2 
3.8 
4.4 
0.7 
5.5 
3.5 
2.1 
1.7 
0.7 
5.3 
1.3 
3.7 
3.3 
39.4 
4,2 
38.9 
7.0 
6.1 
9.4 
6.6 
3.0 
A 


KON 
OF 


-1959 Mean Index - 


tation during the month of April. The April correction index is equal to the 


-Precipitation | Index. —During the spring and summer 
normal: winter atmospheric flow pattern may be altered, so that the runoff 
‘Cola coastal streams not adequately reflect the precipitation over 


; Columbia River Ba Basin in asa whole. The effect of spring- precipitation must 


* Hs» 
— 
a 
: 
4 
&g 
— 
1957 123.0 19.0 
90 
1953 | 123.0 | 3 
1952 | 138.0 | 1 13.2 
1950 | 125, 18.0 
1949] 121.9 | 1 
q 1944 4 | 8. 
1942 | 102.3 6.6 
1940 | 79,3 35 
149.8 6 
1937 63.4 5.6 
| 81.9 | 
1933 | 108.1 
61,5 
1929, | 61. 
— 
& 
= 


fore, a epring- for the ‘Columbia Basin asa whole was 
computed, using USWB and Canadian climatological stations having records 


for the period of study and presently reporting on the Columbia River | Co- | 


October ‘March 
2 (3 


i 


KH WO 


0 


wo war 


he 2 In 1,000,000 Acre 
b Corrected for storage in Kootenay, Hungry Horse, Flathead, Pend © | "Frank- 
_ lin D, Roosevelt, Coeur d’Alene, Chelan, and Brownlee Reservoirs ee Ee 
© product, Coastal Winter-Flow Index (Table 2) times October through March Run-— 


operative Hydrometeorological Reporting Network. There are sixteen ‘such 
a stations, the locations of which are shown on Fig. 1. Inspection of the geo- 
: _ graphical locations of the stations reveals that they are fairly well distributed | 
Bord respect to basin coverage. . The arithmetic sum ‘of the precipitation in | 


inches ; at all stations for the months of May and June for each year was com- 3 


= — i 
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puted, and the average for each month we was » determined. ‘These values are ° the 
-_spring-precipitation indexes for each of the months and are listed in Table 6 : 
4.—WINTER TEMPERATURE INDEX® 


November 


Temp. Cor- 
Average, 
Degrees F | Index 


1950 


(1948 


@ 5-Station Average: Revelstoke, ‘Columbia; Missoula, Mon 


oe 


oa 
o 


eee 


9 
9 
9 
1 
2 
3 
.0 
6 
6 
6 
2 
>) 
6 
8 
0 
8 
9 
8.2 
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Ft 


ge 


BSF 


Cr 
a 


OLR DWWOM 


WY 


— a 


Yakima, Wash,; Boise, Idaho, 
7 4 b Mean Monthly Temperature in Saionie Fahrenheit, 
© Product, Coastal Winter- -Flow Index (Table 2) times aver: anni temperature (divided 
w 100), 


_ The total : index is the sum of the May and June precipitation for each year, and — 


the May-June correction. index is” the product of this total times the CFI. 


ANALYSIS AND RESULTS» 
«At is indeed fortunate that this index is feasible for one - the principal — 
rivers of the United watch water supply is. a prime 


iad 


&g 
‘Water Year | October | 
826 7 36.5 | 37.0 
Bae 
1944 | 5 36.1 | 38.9 
| 398 | 39.5 o 
a t.; Spokane, Wash.: 
+ 
. 
i 


March, 
requisite optimum use of multiple- reservoir s storage e. Table 
_ presents: comparative area-elevation data for the Columbia River east of the 
Cascade and the coastal winter-flow index drainages. On the basis of the 


- average for the four coastal winter-flow index drainages, a smaller ‘percentage 


Monthly Recession Generated April Correction 
Water Runoff, 1,000 | Correction, Runoff Index? 000 
Acre Ft 1,000 Acre Ft | 1,000Acre Ft | Acre Ft 


Product, Coastal W inter- Flow Index (Table 2) 


times generated runoff 


area lies above 3,000 ft on the drainage areas 
«6, 000 ft in the Columbia River Basin. In addition to the elevation differences, 
there are also differences in air-mass_ characteristics which affect the. 


seasonal istribution | of runoff between the basins. Also, when considering 


TABLE 5,~APRIL COASTAL FLOW INDEX CHEHALIS RIVER — 
1953 100 10 90 Ad — 
1935 | 120 -54 86.2 — 
— 
i 
4 


virtually no runoff. _ As a result nearly all of the runoff in the Columbia an 
is from areas above 4, 000 


es comprise the arid Columbia and Snake River plains which produce ; 


Water Year | Spring Precipita- May-June 


008 ‘ad 
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.93 
28 
00 
32 
13 

.80 

.48 
69 

.86 
89 

61 
67 
36 
55 

+ 


2.003 


| 


16-Station Average, i in inches 
© product, Coastal Winter-Flow Index (Table 2) times May-June Precipitation Index — 


On coastal. index streams about 90% of the runoff is generated from areas 
oe 4,000 ft in elevation. These differences are reflected in the average 
flow by the fact that 28% of the 
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— 
TABLE 6,—COLUMBIA RIVER BASIN SPRING-PRECIPITATION INDEX (MAY - JUNE)? 
| 
4 
g 
tm 


Mareh, while about 73 the annual runoff is produced on the coastal 
: tributaries for the ‘same period. The difference between the two represents, 
principally, the differences in snow-water storage in the respective basins. 
_ An approximation of annual runoff amounts reveals that on the average, the 
a total runoff west of the Cascade Range in Washington and Oregon (about 63 i 
3 miles) is about 150,000,000 acre-ft, as ; compared with 140,000,000 acre- ft 
of of runoff for the Columbia Basin east of the Cascades (237, 000 Sq miles). | 
Derivation of Forecasting Relationships.—The initial work on the coastal 7 
flow index method for forecasting Columbia River runoff® was based ona 
_ Sraphical statistical solution through use of a coaxial diagram.9 9 One of the 
_ ‘principal reasons for the choice of this method was that it could ‘account for 
inter- -relationships between the “primary” and “secondary” independent v vari- 
ables to be expressed in the forecasting relationship. For example, the effect 7 
a _ on runoff of spring precipitation is a functionof the winter r snow accumulation — 
as well as the magnitude of the e spring precipitation itself. - 


Areas __Elevation__ 
q 


‘Columbia River near The Dalles 
Willamette River at Albany, Oreg. 
is River at Ariel, Wash, 
Wilson River near Tillamook, Oreg. _ 5 
_ Chehalis River near Grand Mound, Wash. 


- Average ge for four coastal winterflow 


= 


_ With the advent of the electronic digital computer, the solution of a 
“linear ‘regression equations became fast and relatively simple for statistical _— 
analysis. By combining certain independent variables as product functions, it 
- is possible to account for variation in the dependent variables resulting from — 
inter-actions between independent variables. This technique was utilized 
7 by the Snow Investigation Unit in the multiple linear regression analysis of 
snowmelt rates, where there was known to” ond a physical law relating two. » 


independent variables asaproduct. . 
a _ For the multiple linear regression analysis presented herein, the primary : 


_ independent variable is the CFI, which indexes the total moisture input to the 
Inasmuch as is primarily in the form of the CFI 


variables of winter runoff, winter temperature, April runoff, and May- a. 
= precipitation, are all dependent on the magnitude and areal extent of the re 
snowpack in their effectiveness in producing runoff. Accordingly, the secondary _ 


independent variables are e included as product functions the (CFI a 


McGraw-Hill Book Co., "New York, 1949, A. 
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rational combination as to consider the variations in the 
we properly. After having arrived at this concept, it was found that cor- >. 
relations involving the “secondary” independent variables as the product — 
functions with the CFI actually provided stronger ‘correlations than when they 
were inserted as simple independent variables in the multiple correlations, 
Multiple correlations based on this principle were computed for several — 

_ combinations of independent variables for the period of record. The parameters. 
"were finally selected for use in the forecast equation on the basis of indexes 
== rationally account for the various phases of the water balance for the the : 

basin, They do not necessarily represent the “best fit” of historic data of all — 

‘parameters: investigated, which might be selected without due regard to 

_ physical explanation in terms of known hydrologic principles. 

_ The forecast equation was computed by use of the “stepwise” multiple 

ae regression program10 for the computer. All variables affecting runoff, — 

including spring precipitation, | were inserted in in the regression analysis, in 

_ order to obtain the most realistic weighting © of each of the variables. The 

equation so derived is equivalent to the forecast equation for conditions , 


known > as of July 1 of each 1 year. Forecasts as of April 1, , May 1, and June 1 - 
were computed by using mean values of each of the variables affecting runoff — - 


after the date of the forecast. The basic forecast equation <A ——* 


0.2836 - 0.678 + 6. 50 Xs + 29. 


When is thr September runoff for the Columbia River at The 


‘Dalles, in millions of acre-ft; Xy is the coastal winter-flow index (CFI) in 
percent (Table 2); XQ is the winter ' runoff correction index, in millions of 
-acre-ft (Table 3); X3 is the winter temperature correction index in degrees, 
F (Table 4); X4 is the April generated Grand Mound runoff correction index, — 
in thousands of acre-ft (Table 5); and X5 is the May- - June precipitation cor- 
rection index, in inches (Table 6). Table 8a lists the values of each of the 
parameters, - together with the Sevenent and observed runoff, for each of the 
years of study, based on the solution of Eq. Results of this 
present a measure of the effectiveness of all parameters used in estimating 
seasonal runoff, ‘and th ereby represent forecasts that would have been made 
as of July 1 of each year. Tables 8b, 8c, and | 8d s' show similar data for fore- 
‘casts which would be made from data ‘available as of June 1, May 1, and April 
-lofeach year, respectively, May, and 
‘The mathematically derived forecast equation can be shown in graphical © 
form to illustrate the relative magnitude of individual variables in accounting Ds | 
for runoff variance, the ranges of each of the variables for the period of :_ 
record, and the effect of the product function in weighting the “secondary” _ 
_ independent variables. Fig. —-2(a) is the graphical solution of Eq. 1, for the § 
winter components of the forecast equation, for which the CFI is plotted as = ; 
the abscissa and April through September runoff components are plotted as 
ordinate values. . The parameters of f October through March runoff and me 
winter temperature index are shown above and below the central base line. 
a values shove & nail central base line represent the evaluation of the © 


_ 10 “Operating Instruction for ‘Use of the Stepwise Regression Procedure on the | IBM 


650, File 4260,” by M. A. Efroymson and O, Schriker, Esso Resenewe and Engrg. Co., 
Process Research Div., Linden, J. 
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runoff. Ordinate below the base line for specified conditions of 
and winter temperature | index provide the solution to the third term in Eq. 1. 
shows similar ‘graphs for the effects of the indexes of ‘April 1 runoff 
for Chehalis River near Grand Mound and May-June precipitation, combined — 
as product functions of CFL, 
The means’ values and extremes of each of the | index values for the 3l-yr “ 
period, 1929 through 1959, , are indicated on the graphs. The forecast runoff 
for a specific set of conditions can be obtained by the algebraic sum of the > : 
ordinates of each of the four components shown in Figs. 2(a) and 2(b), plus 
the » regression constant, 29.6 (all in unit of milliolns of acre- -ft). In application, :. 
of course, it is simpler to solve Eq. 1 arithmetically, but the primary purpose 
of the graphical solution is is to illustrate the effect of | the product functions. a 


analysis as simple variables, rather than as product functions, the parameter 

lines in Figs. 2(a) and 2(b) would be parallel. It is reasonable to assume, for 
these particular indexes, that their effect on runoff is increased with increasing» 

~ amounts of winter snow accumulation, mainly asthe result of the greater snow 7 

‘covered area that would generally prevail. Accordingly, the diverging shape 


of, the parameters as shown in Figs. 2(a) and 2(b), resulting from the product 


Date of Forecast Average Error,| Standard Error Correlation — Coefficient — 

106 Acre Feet Estimate Coefficient, of Determina- 
06 Acre Feet | R | tion D (= R2) 


0,821 


functions, is a more reasonable solution than using — unrelated inde- | 
_ Summary of Results.—The comparisons of forecast versus actual runoff 
listed on Table 8, are plotted on Figs. 3. Table 9 lists statistical ‘comparisons 1 
of. April through September runoff forecasts which would have been made, as 
of April 1 through July 1, based on solution of the forecast equation for the 
_ $1-yr period, 1929 through 1959. Early season forecasts are based on assumed © 
‘mean values of runoff index parameters for those variables which would _ 
the equation subsequent to the date of forecast. The mean April through Sep- 
= tember runoff for the period is 96,420,000 acre-ft, and the standard deviation — : 
of the runoff is 19,030,000 acre-ft. The standard errors of estimate and the 
_ correlation coefficients shown above were computed on the basis of following 


‘numbers of degrees offreedom: 
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RIL THROUGH SEI 


CAST RUNOFF, AP’ 


FORECAST RUNOFF, APRIL THROUGH SEPTEMBER IN 1,000,0 


> 

“March, 1961 


9 t 
Standard error = 5.55 x 108 
_-—- Correlation coefficient, R = 0.960 a 


> OBSERVED RUNOFF, APRIL THROUGH SEPTEMBER, IN 1,000,000 AC. FT. 
COLUMBIA RIVER AT THE DALLES, CORRECTED FOR MAJOR RESERVOIR STORAGE 


= 


Average error = 4.69 x 106 acre-feet 
Standard error = 6.47 x 106 acre-feet 


_ Correlation coefficient, R= 0.947 

(See Table 11) 


3 

| 

coLUMaIA AT THE DALLES, CORRECTED FOR MAJOR RESERVOIR STORAGE . 
3(a),(b). FORECAST VS, ACTUAL RUNOFF 


1,000,000 AC, FT. 


= 
Average 6.48 x acre-feet 
error = 8.48 x 108 acre-feet 
Correlation coefficient, R= 0.906 
(See Table 12) 


FORECAST RUNOFF, APRIL THROUGH SEPTEMBER IN 


MAY 1 FORECAST 


OBSERVED RUNOFF, APRIL THROUGH SEPTEMBER, IN 1,000,000 AC. FT. 


COLUMBIA RIVER AT THE DALLES, CORRECTED FOR seca RESERVOIR STORAGE 


'HROUGH SEPTEMBER IN !,000,000 AC. FT. 
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RI 


Average error = 6.72 x 106 acre-feet 
Standard error = x acre-feet — 
9.34 x 108 “feet 


Comeletion R= 0.680 — 


FORECAST RUNOFF, AP! 
FORECAST 


5 OBSERVED RUNOFF, APRIL THROUGH SEPTEMBER, IN 1,000,000 AC, FT. 


ST Vs vs. RUNOFF 


AT THE DALLES, CORRECTED FOR MAJOR RESERVOIR 


March, 1961 
‘The en errors for each of the Nee efficients for the ne variables i 
in Eq. 1 are asfollows: 
Variable 


of Results With Other Index Procedures. —The prime forecast 
of seasonal runoff in the Columbia River Basin has been April 1, partly be- 
cause snow survey data as of that date has the most general coverage, i 
also because that date represents the time at which a fairly reliable forecast | 7 
7 of streamflow may be determined and yet is s early enough to take into account 


4 operational problems that may be incurred onthe basis of the forecast runoff. 7 a) 
In 1954 the technical staff of the Water Hpaemy* Subcommittee, Columbia 


of four procedures in use at that time for 
the Columbia Basin, based on April 1 forecasts are shown. lla 
ee The procedures were developed by four United States government : ponpenomndl 
: including the Corps of Engineers, USWB, Soil Conservation Service, and USGS. 
_ They were based primarily on snow survey or winter precipitation data, but 
with minor adjustments for other runoff factors. Inasmuch as they were de- 7 
veloped for different forecast periods, a consistent measure of aque 


was the standard error of estimate, expressed in percent of the average runoff 


for the —r"" period. These values ranged from 11.0% to 13.2% for April 1 
-forecasts.+ lla For the CFI method, the standard error of 9,340,000 acre- ft, 
in terms of ‘the average April through September runoff of 96,420,000 acre-ft, _ 
is 9. 7%. The correlation coefficients of the four procedures ranged from 0.71. 


to 0.88, based on April 1 forecasts. The omar winter flow index method has 


The oi this of forecasting snowmelt runoff from the 
; Columbia River at The Dalles is based onthe geographic location of the basin — 
with respect to the large-scale air-mass and atmospheric flow conditions - 
_ during the winter season, as well as the differences of the topographic char- 
~ acter, principally with respect to elevation, of runoff-producing areas, between | 
»,: the Columbia Basin and the adjacent coastal regions of the Pacific Northwest. 
_ Because tl the method is based on an index procedure, some consideration — 
_ should be given to the various indexes presently used for forecasting runoff. 
_ In general, precipitation indexes, either from rain gage or snow course water- 
J equivalent measurements pies based on samples which measure id an sand 


st Seasonal Runoff of Columbia River neat an 1e 
Dallas, stent ” Columbia Basin Inter- -Agency Committee, Water Mgt, Subcomm 
Portland, Oreg., 1954, +1954, (a) Table 9, 49, 
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RIVER RUNOFF 
. the moun mountains ns of the wes west, where there is great areal ‘variability of pre- 
cipitation, the relation between gaged and actual precipitation amounts may 
Be. from year to year as the result of meteorologic factors which, for the : 
most part, are too complex for a rational analysis by existing techniques and 
_ intensity of the gage network. As an index, then, precipitation measurements 
aa account for the | random variability \ caused by inherent deficiencies of 
measurements | by gages or snow courses, as well as the variability caused 
by the character of individual storms. The coastal winter-flow index, though 
having limitations, possesses certain advantages over gaged precipitation 
_ measurements. It samples a much larger proportion ofthe total water supply, 
illustrated by the fact that in the Columbia Basin, the ratio is about one 
part in 14 in terms of runoff. The accuracy in the measurement of streamflow | 
is an inherent limitation, but such deficiencies also apply to other index pro- 
cedures» for "forecasting runoff. The principal limitation is the variability 
from year to year of the average air-mass characteristics, and the average 
atmospheric circulation pattern during times of significant precipitation. 
The -winter- -temperature index and Columbia River winter runoff account 
for part of the variability caused by the air-mass differences. They occur a 
from year to =. As indicated by the comparison of forecast | errors, the 


_ better than precipitation indexes as measured at presently operated rain gages 


_ or snow courses, , for the Columbia River Basin as a whole. a Paes” pat 
_ ‘The effects of spring precipitation are accounted for in this procedure by 
the April generated runoff index for the Chehalis River near Grand Mound, | 
and the May- June precipitation index, based on a 16-station mean. Use of 
‘statistical techniques for weighting individual indexes of runoff should 
predicated on inclusion of all variables in the multiple regression analysis, 
By correlating runoff directly with just those values of the variables known 
_as of April a. it is possible to obtain an apparently better correlation for r April 
1 forecasts than using the correlation for all variables, and then computing 

forecasts for April 1, using mean spring precipitation. Omitting the spring 

_ precipitation effects in the derivation of April 1 forecasts is unrealistic and 
simply forces a best fit of the historic data with incomplete information. It 
does not necessarily provide the best weighting of individual winter indexes. 
eo. In the particular case of the winter indexes of runoff for the coastal © 
winter- -flow index method, correlating them directly with April through Sep-_ 
tember runoff ‘provides a standard error of estimate of 9,190,000 acre-ft, i 
compared with 9,340,000 acre-ft for April 1 forecasts, obtained from the 

_ general forecast equation and using mean values of spring precipitation in- 
dexes. ‘The: apparent reduction inthe standard error of estimates is fallacious. 3 
Similarly, direct correlation of indexes known on May 1 of each year would 
provide a standard error of estimate of 000 acre- -ft, compared with 

‘The of product functions in multiple ‘linear regression analysis of 
seasonal runoff variance may have application | to forecast procedures based 
on. snow course or precipitation gage indexes. Product functions should be — 
used only where there is known to be a physical process relating two inde- | 

pendent variables asaproduct, 


Application to Other Basins. Since this procedure deals with largescale 


meteorologic phenomena, | its application is intended primarily for large 
basin areas , rather _ smaller ee areas, However, where there are 
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at pack in the middle latitudes on basins smaller than the Columbia River pro- © 
vided that the index area is well centered air stream dominating 


The coastal in conjunction with other runoff parameters, 
"provides a basis for forecasting seasonal (April through September) runoff 
volume for the Columbia River at The Dalles. The historical streamflow and 
. climatological records used in the correlations have a sufficient period of 7 
record: (31 yr), considering me number of variables used (six), to insure a 
reasonably sound statistical analysis. ‘The use the “stepwise” ‘multiple 
linear regression program for the electronic digital computer expedited the © = 
a correlation analysis used in developing the basic forecast equation. The basic | 
. approach of correlating runoff (the dependent variable) with indexes of all of 
significant factors affecting runoff (the independent variables), is 
considered to be hydrologically sound. Apparent improvements in fit of his- 
= data by direct correlations with runoff of variables known as of of the date 
of forecast (that is, April 1) are considered to be fallacious. Some be 
The use of product functions of the primary and secondary 
variables improved the forecast relationship, in accounting for runoff variance. ‘ 
; Product functions are believed to be valid in cases where the effectiveness of 
5 one independent variable is known to increase with increasing values of a 
- second independent variable. Comparisons of forecasting methods for Game 


in the Columbia River Basin, have been made by the technical staff of the 
| Water Management Subcommittee, Columbia Basin Inter-Agency Committee 

ca Use | Of these comparisons shows that the coastal winter flow index method 
m « as accurate with regard to fit of historical data as other procedures based 


River at The Dalles, developed by government agencies for operational te 


owmelt. 
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_ EFFECT OF AQUIFER 1 TURBULENCE ON WELL 


; _ JOE L. MOGG. 1. Mr. Philip has called attention to a concept that was new — 
to the author. He pointed out that inertial forces, rather than the inception of 
turbulence, initially caused the head loss to vary w with an exponent of the ve-_ 
locity that was greater than one. 

Mr. Philip’s excellent ‘and numerous plus some others were 

reviewed by the author. G. Schneebeli2 describes experiments that also show 
that | the head loss begins to vary with an exponent of the velocity greater than 2 
one after a critical Reynolds number is reached. However, dye introduced in 
the flow, for the purpose of identifying streamlines, showed that the flow was 
laminar and not turbulent in this non-linear region. Turbulence was evident 

at higher Reynolds numbers when the dye became readily dispersed. Thus, -_ 
_ these experiments show that there is a broad region of flow between laminar- - 
linear flow and turbulent-non-linear flow that is laminar, but not linear. The io 

cause of the non- linearity could well be due to the inertial forces. 
wi The conclusion is, then, that turbulence seldom occurs in the vicinity of a 


p 
flow. 


- & November, 1959, by Joe L. Mogg (Proc. Paper 2265). 

Field Engr., Edward E. Johnson, Inc., St. Paul, Minn, abe 

_ 2 «Experiments on the range of validity of Darcy’s law and the appearance of turbu- 
‘lence in a flow »” by G.  Schneebeli, La Houille Blanche, ‘March, , April, 1955, 
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NEW TO LocaL. FLOOD PROBLEMS@ 


Closure by Herbert D. Vogel 


HERBERT D. VOGEL, F. ASCE. _—The observation by Bernard Golding, 
M. ASCE, to the effect that communities make their own observations 


and determinations through local engineering forces and consulting engineers, | a 

is sound in principle. The primary need isto obtain adequate data from what- 

_ ever sources may be available and make the best possible use of it. _ a 
7 Legislation has been enacted (1960) by the Congress to permit the United 


States Army, Corps of Engineers to carry out studies of this kind in water- 


sheds and on rivers and tributaries under their jurisdiction. This provides a 
_ tool for the use of communities throughout the entire country and should en- 
- courage the approach to solution of flood problems that has been indicated. 
ig However the studies may be made, it is a certainty that results will be of 
benefit to individual communities during years | of growth and that savings | of 
great magnitude » will be made by the Federal Government in the future. 
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THE FOURTH ROOT nef DIAGRAM 


e by T. Blenc 


: 2 of the yr Nett, and its u: use, , the writer must confess to a third “object in mind” q 4 
_ to add to the two stated by Gerald Lacey, discuss its reason, and then answer - 


include detailed replies automatically. The third object in 
4 
from the of the science of engineering fluid mechanics so 
that they can develop knowledge of pipe flow in terms of their own data and . a 
the logic that they acquire in the exercise of their profession. The context i 


and the references were intended to aid this object, whereas the diagram pro- 
_ vided the means for co-ordinating data. The reason for the object is that the 
_writer is one of those who deplores the long neglect of disciplined instruction _ 
in basic science and mathematics in North American schools, continued into 
college engineering curricula. In his opinion the neglect has been responsible q 
-- the displacement of the applied science of hydraulics by #o-cates fluid 
mechanics “necessarily based on -_evasions, part-truths, dogmatism, and so 
forth. From such a subject no advance in fundamental understanding is pons 
‘ble; in fact it misdirects research, For example, to get away from pipes, , the 
a engineering student is taught a Bernouilli equation that allows for friction by 
- pretending that pressure is the same in : all directions about a point, omitting 7 
7 the velocity head of fluctuating motion, often ignoring the kinetic energy cor- 
et factor, _and lumping all errors and omissions into a “loss of head.” 
_ may imagine the wasted effort ‘if he were asked to conduct fundamental 
- eoearch into the hydraulic jump on a 45° spillway chute at a place at which * 
4 velocity had become terminal. Here the shear stress on the floor equals ll - 
ES stress and one principal stress is zero, the shear stress along the 


floor in the jump is high, the kinetic energy and momentum correction factors | 
” i. and the turbulent velocity head at the control sections all merit attention. i 
he had no fluid mechanics training at all, but good fundamental dynamics and 
7 physics, his chances of success would be reasonable because he would know — 
how to wate his problem in (ran of its nerves factors, translate it into the | 


a gJanuary 1960, by T. _Blench (Proc, Paper 2340), 
Prof., , Dept. of Civ, Engrg., Univ, of Alberta, Edmonton, Alberta, 


7 | 
a 
ig _— should be tackled inductively or deductively with or without expert aid. = 7 
lsc’ if equipped to study hydraulics as presented 


March, 
in 1 texts40, 41, 42 that presume an adequate foundation of mathematics and 
physics, or draw attention to assumptions, idealizations and limitations and 
do not avoid an argument by “it can be proved that” or “Seen to so- and- : 
‘To avoid unintended offence to textbook writers the author states that he 2 a 
places blame squarely on neglected basic science education, has 
errors because of it, has experienced the difficulty | of teaching inadequately — 
prepared students, and cannot a publisher accepting a non-conforming 
_ college book. He knows there is a changing attitude to education and that at 
least some colleges, including his own, are making” desirable innovations. 


7 ~agaed of the preceding opinions an attempt will now be made to ——« 
—— 


‘The practical problem may be stated as “A Newtonian liquid flows el 7 
conditions of uniform dissipation of head (energy per unit weight) per unit | 
4 length in a circular pipe of given material. Derive a formula to express the — - 
ee mean speed of flow, U, interms ofS, p, y, and d and the boundary properties. e. 
An omniscient being, OB, if asked to solve and explain, would probably | 
aed from the laws of intra and inter molecular action and combine them with 7 
eee form of Newton’ s second law. These would show that tt the motion — 
- would not be steady, in general, but the flow would “roll itself up” into eddies = 
_ whose statistical mechanics could be worked out ‘go that the probability of a 
i specified velocity at any point during any prescribed time interval could be | 
deduced. Thence the mean behavior at every ‘point could be determined d and 


‘fi lly t t t t i 


from ‘which is. for dimensional reasons, and in which OB would 
¥ _ have to explain how to measure “roughness. ” His answer would be the = 
namically correct as it would have been deduced correctly from unas-— 
~ gailable premises in terms of undisputable laws. An ordinary hydraulics : 
-vestigator, HI, of the classic school would admit he did not possess OB’ Ss 
omniscience ‘so would proceed like his iis brother physicists who also ha had to oun 
‘ravel the mysteries of “fluids” (for example, electricity, heat) whose internal 
_ ‘mechanism was not observable. That is, he would correlate measureable 
; _ quantities covering the whole available range, and test formulas fitting mad 
results: to see if they suggested simple relations : among dynamical entities - 
force, _ rate of work, and so forth. His knowledge of general physics would © 
Cause him to believe that if he were lucky enough to hit a formula identical —_ 
with» OB’ ’s true answer it would be simple and would certainly express some — 
simple relation among dynamical entities. Such a relation might be reducible 
to minimal type, analogous to “light travels SO as to make the time of transit 
7 a minimum,” “free flow over a weir occurs so | that the energy of flow is a0 
- minimum at a critical section... .” Although he would hope to discover the 
dynamically correct formula by his indirect approach he would expect the _ 
Odds to be against him, particularly if the measureable quantities were not | 
well defined or did not cover an adequate range. As a scientist he would not 
a “convince himself that he had made a discovery until he.had checked it against | 


_ 
40 «Fluid Dynamics, by L. Prandtl, Blackie & Son, 1952, pp. 125-130. 
41 «The Mechanics of Turbulent Flow? by B, A. Bakhmeteff, Princeton Univ, Press. = 
42 Pluid Mechanics and Heat Transfer,” ” by J. M. Kay, Cambridge Univ. Press, New 
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DISCUSSION 
available information and verified th: that t it contained no physical contradic- 
; tions. Then he might allow himself to state that he had found a relation that 
; ‘seemed to have, or suggest, certain dynamical meanings within a specified 
; . range of variable. The relative newcomer on the North American scene, the © 
_ enthusiastic fluid mechanics investigator, EFMI, does not have HI’s inhibitions — 
; but feels free, on occasions, to conjure formulas from the air, consider physi- 
cal absurdities in them to be unimportant so long as the answer looks practi- _ 
_ cally good, produce only enough data to justify his assumptions and and extrapolate 
_ HI has developed two main lines of attack. One is 3 to regard the flow formula 


for various values of an unmeasureable relative ‘This leads 
to the standard friction factor diagram and the discovery of different phases | 
of turbulent f flow. The practical and basic value of this diagram is that it is 
plotted from readily measureable variables, except for “ roughness.” So, f 
data can be separated into groups of approximately constant roughnesses, it 
gives a means of defining roughness in terms of whatever formula HI likes 

to try as his current guess at OB’s and, of course, it shows him the range of 
conditions over which “roughness” is relevant to the flow. Then, if he finds 


that facts justify describing roughness by a “ roughness height” e in the rough 


Eq. 5 if he favors the fourth root type of formula, or by the logarithmic form- é 
ula that EFMI likes, , or by any reasonably fitting formula; whatever formula 
he likes can be used, after algebraic conversion as with Eqs. 2 and 3, to draw ' 
lines of equal relative roughness on the friction factor diagram in the ni 
= and to label them with appropriate e values. For the purely practical 
man, any system of e’s is good enough provided the formula on which it is 
based is” a good practical fit, for fixed e, in the range of data of interest to 
him; but, of course, HI will never be satisfied till he thinks he has hit OB’ s 
dynamically correct formula that would fit the full range of its qyaamical 
applicability. The second line of attack fo follows logically from | the first oa ; 
cause the phases of turbulent flow shown by the f diagram direct al 
immediately to the question whether they are associated with some paoenrny . 
of the velocity distribution, Observations of velocity distribution, although _ 
very difficult | near the boundary at - which conditions are important for ex- 
_ planations, showed the existence of a thin layer of non-turbulent flow for so- © 


called “smooth poe and indicated, | as 3 implied by D’ Arcy long 6 280," 


ordinates all welnalty data plotted from their turning points as origin would 
indicate one curve there (See Fig. 3). A logarithmic curve is a good fit to the 
inner zone and to the wall zone at which one curve is not strictly a fit. 
-Prandti40 has indicated conditions (not those of a circular pipe) under which 
a velocity distribution curve could be deduced to follow a logarithmic law | 
5 
over a range. A parabola is a good fit, according to Stanton,” over most of a 7 
pipe but does not attempt to fit the non-conforming wall zone. A hyperbolic © 


| fanetion, according to the writer, fits a trifle better than a parabola. No doubt 


mathematicians could find comparably good fits. Now, with a formula fitted 
to a velocity curve a different way of looking at a flow formula becomes — 
venient. flow formula may be regarded as the consequence of integrating 
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1961 
the velocity distribution mn (not the deficiency one) over er the cross. section to to ob- 
tain Q, by the cross ‘sectional a area to obtain qf the formula 


_ must be practically useful over the same ranges of variables, However, as the 
velocity formula remains speculative, the chances of deducing OB’s dynami- 
; cally correct flow formula from it seem remote. EFMI wrongly accepts the 
logarithmic fit to the velocity deficiency distribution over the whole pipe as 
a absolute dynamical truth and uses it with various assumptions to deduce a 
logarithmic flow formula. This procedure he wrongly attributes to Prandtl ~ 


3 


With the preceding background we can now , outline attacks along HI ‘HI lines 
_that should automatically answer specific parts of the comments. 
First let us— consider velocity distribution in a circular pipe under the 
given conditions and use dimensional Then we we may write: 
R ey v). 
which y is the distance of of from 
_ that the boundary can be } wepences by a geometric surface), R is pipe radius 
(with the same proviso), u is the time-mean velocity at y,T 9 is the time-_ 
mean shear stress _at the boundary. Given the mathematical ability we could 
always eliminate -e from Eq. 23 by inserting its value from Eq. 24, and the 


a fit then the flow formula deduced fro it by integration 
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DISCUSSION 
max To P/V) (25) 
in(y/R, Umax/Vx» Vx 
Now, suppose the smooth- delineate shown by the f diagram, in which “) 
can be dropped from Eqs. 23 and 24. Then we could eliminate Umax from 
them and Eq. 26a would reduce to: 
It does not seem eenteieieinatiaimiaiis this equation and Fig. 102 of reference 
15(a) quoted by Mr. Rajaratnam because that figure is fitted by a single line ~ 
implying that y/R can be neglected. Also, if that line is really the ultimate > ” 
truth, then we must believe that du/dy has a definite value at the centerline 4 
of a pipe, that is, there is a cusp in the velocity distribution. It is interesting — 
to note that the line has an equation that gives a mixing length Geteintion 


serting measured values in the formula. that ‘defines mixing length. 43 ‘3 The 
writer emphasizes that he is not seriously disputing practical goodness of fit © 
4 of me line over the range of its data, but is disputing its dynamical exactness: 


has ‘been discovered beyond doubt, unlimited extrapolation isp permissible, 


iy Reverting now to Eq. 26a we note that it is consistent with the general ex- 


ng 


perimental finding that, for any boundary, there is a universal velocity de- 
& ficiency curve is mot too small, it does not 


“yy, 


wee 
experience shows that viscosity can cease to affect certain 
problems appreciably _ when -Reynold’s number exceeds a certain ‘amount. 
as we did to derive 26b from Eq. 26a, but eliminating 


would ‘not the possibility of nature 
for this smooth boundary case. ‘ase, However, ‘the function would be ¢ discontinuous 
because the existence of a laminar z zone has been demonstrated experimental- — 
- ‘i This seems to sum up all that is definitely known about velocity distribu- ; 
tions. For practical reasons, or for analysis, various deductions may be made 
"trom various assumptions, but the results cannot be dynamically correct gen-— 
erally except by fortunate cancellation of errors. For example, the meen 


... “The Mechanics of Turbulent Flow,” by B, A. Bakhmeteff, Princeton Univ. ‘Press, 
Princeton, N, J., 1936, p. 78, Pig. 49, 
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"March, , 1961 
0 of the laminar sublayer . . given by the equation . ‘awl 1)” - that Mr. Rajarat- 

nam quotes on the authority of reference 15 is one of many quantitative defi- 
nitions possible in terms one of many possible assumptions. It is derived by 
assuming that the laminar zone (reference 15, Fig. 102) is fitted by a straight 
— To = yu du/dy, instead of a parabola as it should be for dynamical and 
functional correctness, and that the turbulent zone is fitted by the logarithmic © 
_ formula that we have just demonstrated cannot be dynamically correct. The | 
. of these two is found by trial and error solution of V«6/y = N- 
2.5 log N+ 5. 5. Some proofs seem to make a kind of physical argument out 


of the constancy of V, é6/v , and thence claim special validity. Actually — SS 
= curve that fits the closely bunched points of Fig. 102 reference 15 will 
have a single intersection with any single curve that is supposed to fit the 
laminar film so will give a fixed V, é/v for the intersection. It should not be — 
forgotten that the laminar film is of indefinite and fluctuating size that cannot — 
wa measured, so its “thickness” is merely a conventional measure. However, 7 
the convention may be based on a dynamically | correct relation or not, so the 
definition might lead back to something useful. The writer chose to define 6 
in terms of an argument of Prandtl’s adapted from a slightly different case;® 
_ allowing for the | actually fluctuating motion in the laminar film the argument — 
ean be used almost word for word. The different cetations will not agree nor 


= results deduced from them. erfect ded ry. ) 


Another class of - quite ‘perfect detection concerns flow | formulas and 
s 


es 
- and averages the terms over the | cross section wy multiplying by 2 mr x ae 


integrating over the whole cross section, dividing by No matter 
- what one assumes the average value of Uma, is itself, and of u is U. The 


_ average value of the right side depends on the function and, of course, a. 
7 experimentally determined values of the parameters in the function. Ifthe 
7 function is taken as A (r/R) 2. in which A= 6.56 agrees with Stanton over (80% a 


r of R, the average is A/2, from the property of a paraboloid of revolution that _ 

- the mean ordinate is half the maximum, If the function is the standard loga- 

: rithmic, - minus 2.50 loge y/R, that attempts to fit the non-conforming wall 
zone, the average is 3.75 and is practically more accurate than 6.56/2 = 3. 28 
that deliberately ignores the wall zone. The author has found that 13 [cosh(r/ R) 

_ ] fits Nikuradse data very closely over 85% of R and gives an = of 

3.435 (See Fig. 3). From the viewpoint of aiming at finding OB’s dynamically 

.. perfect flow formula it is important to note that any function, fn(r/R), , will” 

which B is a constant, e even the whole is repeated for a 
broad channel the only difference in Eq. 31 will be in the value of B (see Eq. 

7 (34) . However, as was shown via Eq. 26a, the right hand side of Eq. 27 is not 

quite a function of r/R alone, so B should be slightly variable, and some doubt - : 


B= constant. To preceed from Eq. 31 towards deducing flow formulas we 


us 
the definition f = 2 g — and the he dynamical fact that T> = 4 g dS, so a 
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converts to: sto: 
Umax/Y 
Eq. | 32a shows that no further prc progress towards dotucing the form of f is 

_ possible without knowledge of the relation of U max/V, _ to the factors on 
which it depends; this: | knowledge is not contained in the velocity deficiency a 
Eq. 27. The HI method of obtaining the knowledge would start by seni 
data and studying them. The writer would not like to assert, without extensive 

research into literature, that such data as were obtained preceded a “mer 
to believe the logarithmic type of answer of engineering textbooks. For “ cough 
boundary” one can quickly obtain Umax/V, and e/r data of the Nikuradse’s 

_ experiments on artificially roughened pipe from the diagrams in reference a. 

If the reaction between the two variables is to be proved logarithmic the Jia 
are plotted on single log paper and will give an excellent straight line. If — 
are to be be proved d exponential the data are plotted on double log paper and will 
give an 1 equally excellent straight line of slope 1/8. If the equation of the for- 

- mer - line is inserted in Eq. 32a the standard logarithmic expression for fin — 
terms e/r will ‘result. the of the latter line is inserted it will 


latter’ s functional form ‘even were it lucky enough to be the ultimate 
_ dynamically correct one. (An important lesson of these alternative plots is | 
that function can agree a | suitable power one over large 


seem to have been attempted “Instead, the sequence of indirect 
duction seems to have been that explained neatly in reference 12, Sections (25 
and 28. Briefly, Prandtl deduced velocity distribution along a a plane wall in an 
\ infinite fluid with constant shear stress throughout and far enough from the 
wall : for viscosity to be ineffective (reference 40, pp. 125-129). Having deduced 
it as logarithmic he proceeded, by an ingenious but not quite convincing argu- 
; - ment, , to extend it to fit the zone where viscosity mattered and came out ~ 


t= k u/V, = log V, y/v +C, .. 


He then pointed out that it would fit ybmradee’s ’s ‘experiments on smooth pipes" 
fas an approximation” and that “for flowinpipes” there was actually a “devi- 
. = from theory.” The writer does not know of his having recommended or 
supported, as dynamically exact, integration of this flow formula to give a 
_ mean value of U, and thence an expression for f. However, it is this integra- 
tion with-the constants fitting Nikuradse’s smooth pipe f data, that gives the 
logarithmic smooth boundary line on the Moody diagram, and its formula is | 
found associated in textbooks and articles with the names of Prandtl and von | 
Karman in a way that ; seems to have caused engineering readers to believe | 
: that these hydraulicians deduced and recommended it. So, when Mr. Rajarat- 
ham says “it is widely known that for values of the Reynolds manber higher 
than 105 the flow resistance is better approximated by the Prandtl- ‘Karman 
equation in the writer feels that he has been misled by ESMI and the state-— 
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‘ment should read “it is widely mis- stated by textbooks that. . . by the sie 
‘named Prandtl- Karman ee All that the Nikuradse smooth-pipe f value 


enough Reynold’ If his “smooth”, pipes had been replaced by even 
- “smoother” ones his smooth line should have been expected to start turning = 
up at higher Reynold’s numbers still; if they had been a little less “smooth” 3 
then presumably he would have used them to obtain the constants in the loga- 2 
- rithmic formula although they now pass off as extensions of the Colebrook- 
_ White ones and not “smooth.” Problem 5 was framed deliberately for the in- v 
mentee reader so that he could decide for himself whether there was some- 
thing dubious about the logarithmic smooth boundary formula at Reynolds" 
Numbers higher than those used in the lab. Problem 4 was to convince him © 
that, if the Moody Diagram > is "wrong, it may be alarmingly ‘wrong for very 
A definite quantitative anawer 10 | the queries about canal distributions is 
prevented by inadequate data. As the writer’s interest in flow formulas arose 
from the possible - dynamical implications of Mr. Lacey’s regime formulas, 
he has never been impelled to arrange observations of his own to determine 
the constant in pipe velocity deficiency formulas that would suit channels. His > 
inquiries from laboratory workers on channels and his long field experience 
with canals support, very strongly, the view that the functional form of the | 
: _ Velocity distribution down the centerline of a broad canal is that along the a 
radius of a pipe. Narrowness probably | does not change the functional form _ 
‘ — by introducing an extra non-dimensional ratio such as b/d and altering ie 
| the value of a parameter. Exact determination is difficult in laboratory flumes 
_ because they are narrow and the absolute depth is small. In the field current © 
meters are “useless — near the free surface and the bed, and many beds are 
‘mobile. Moreover, even inlined canals the water is often turbid and laboratory 
_ Observations show that turbidity can alter, radically, the coefficient 2.5 in the ' 
———- deficiency formula—one result given privately by a well-known | 
=_— showed a change to 4.8 with a suspended load concentration of 15. _ 
grams per litre, but a very good agreement with 2.5 over 85% of the depth for 
clean water. Finally wind is a major disturbing element in the field. A mathe-_ 
- matical approach to the problem seems to show a contradiction between the 
two common practices of assuming that the velocity deficiency expressions 
are the same, complete with constants, over the radius of a pipe and the depth 
of a broad canal and that, inuse of the Moody diagram or the Manning formula, 
a channel can be replaced by a circle of the same hydraulic radius. Ifthe | 
velocity expressions are the same then the right hand - — $1 becomes: 
that will not give the same value of B for the canal and the pipe. (Here r is 
depth from | surface and R is overall depth). But if this is so then Eq. 32a 
_ shows that the velocity distributions are different for the two cases, in which | 
- it is supposed that V. x and f are the same. This raises the question whether © 
_ the parameter in the veloute formula might differ between pipe and channel 7 
7 in such a way as to make B+ fn (r/R) the same for both. None of the functions | 
of this article seems adequate for the purpose, but a mathematician might | 
_ discover one. _ However it should not be forgotten that B is not ———— 
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7 Tact. Physicany the writer has not been troubled by the difference between 
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surface and center flow of canal and pipe. He has visualized mean surface 


other. But another way he has visualized the mean motion at the surface of a 
channel as if the mirror image were superimposed to turn it into a pipe. 

Perhaps the whole preceding discussion will explain why the writer jae 


the of Speculative velocity distributions to obtain a flow 


to obtain the simple laws of electromagnetism. The cautions procedure seems > 
a to be to use Eq. 32c to fix the constants of a a practical velocity distribution in 
ss: The commentators’ interest in the justification for using channels to aid 
x in testing for a universal flow formula obviously calls for a little poten 
in amplification of the references. The writer’s interest in dynamics and 
_ physics, his years of practical canal engineering, andhis study of the success | 
_ with which Mr. Lacey developed his regime theory along the logical lines of 
- classical physics44 have all helped him to think of the flow formula in terms se 
of either the equating of resistance to driving force or the selection, by na- > 
ture, of one particular flow for any one problem from among the infinite num-_ 


with an interest in ‘physics, to believe that a formula giving the relationship by 
change in ‘the dynamical set- -up. Complications may be added by geometric | 
“either kind of boundary, or in channels of mobile boundary, all carrying New- 
_ rolling-up into an eddying form that does not seem to be altered essentially — 
channels with rigid boundaries? The answer seems to be that there should be 
tried to prove the point from somewhat advanced dynamical principles but 
‘is, a a search had to be made fora thread running through the various competing o 


ber of flows that one can imagine as ‘kinematically possible. ‘It is difficult, — 
dynamical entities in simple observed circumstances will be com- 
: 6c or that its functional form will be altered seriously without a radical © 
conditions, by extraneous phenomena being measured with the ones of. 
_ interest, but behind them is the simple relation. Fundamentally, in pipes of 
-tonian fluids, the boundary “puts on the brake” and the flow responds by 7 
by alteration in the nature of the brake. Why, then, should the functional form © 
7 of the flow equation be different for the three cases or, for that matter, for | 
no difference, but that the term that expresses the “braking action” will have - 
different “detailed expressions for the different cases. At one time the author | 
- lacked the necessary skill. However, the normal HI procedure that he likes 
to think of as “plotting formulas” instead of points, had to be followed. That 
3 formulas for smooth and rough rigid pipes and for mobile boundary canals 
(the last made available in physically suggestive form by Mr. Lacey’s monu- 


mental work, 44 The thread discovered is the equation given in the paper or un- ; 
a der the heading “Universal flow formula,” with its values of x. The “theory 
of flow,” such as it is, amounts to believing, till something more definite | fad 
a turns up, that the odds against this functional form having no physical signifi- 
“cance are very high. Accepting this belief the form merits a try in practical 
a use and, far more important physically, ought to initiate a search for a i 
namical theory that would prove or disprove it. ‘lt is not believed that the ex- 
pressions for ‘x are beyond improvement; in fact, for the regime case, it Jf 
seems very likely that small b/d ratios will require an adjustment. Obviously : 
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“Flow in Alluvial Channels with Sandy Mobile Beds, Gerald Lacey, 
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Many more data are required from the field t to answer Problem 5 finally and 
thereby test the theory decisively over the practical range for smooth bound- 
ary. The data will have to be numerous, and from different observers, be- 


Cause accuracy of measurement will not be too good and the of acci- 


dental bias is high even in laboratory work. 
(1) The Prandtl proof is not valid for a canal. — outstanding reason is. 
a that it rests on the assumption of constant shear stress, but canal shear is 
oe proportional to depth from surface. Flow of wind along the earth’s surface — 
b probably pourttes nearly ideal conditions not too close to the ground and e 
aa (2) Referring to the Webster and Metcalf data, the variation of n ane not 
: _ seem sufficient, in terms of reasonable errors of observation and the uncer-_ 
tainty of where to measure d in a corrugated pipe, to pass judgement on the 
im matter. _ Note that, for _a known discharge, my varies as as the 2. 67 power of di- 


(3) ‘The reference to Morris’ work is appreciated. The author should have 
diagram applied for rigid boundary roughnesses that can be expressed 
_ in terms of a single roughness height.” The oversight is regretted, => La 
‘The recent references by Mr. Van Sickle merit study in original. ‘The 
_ writer believes they emphasize his own view that deviation between the loga- 7 
rithmic and the exponential formulas is too small over most practical ranges — 
to permit a decision to be made between them. For the rough boundary loga- 
rithmic formula the best present test, in terms of real data, seems to be as 
in _“Inconsistence of Manning “Formula and Moody Diagram,” and for the 
smooth formula the best test is to study the few very large conduits, with 
_ Reynold’s Number greater than about 5 x 106, in which conditions are a bit 
= off smooth and there are enough data to show that the transition curve is con- 
sistent with Blasius (in that it could bend into the Blasius line) and inconsistent 
ie the logarithmic curve in that it heads right across it. Thus Fig. 1 of 


reference 29 is for Reynold’s Number well below 106 and indicates only tt that — : 

5 the aluminum pipe appears to be just starting to go into transition on its way _ 

: _ towards roughness; some of the points are below the log line a little and the - 
Figure has" the logarithmic line incorrectly drawn in its upper portion. In| | 
reference (31, Fig. 3 shows four sets of observations with marked inconsis- 

_ tencies that puzzled the commentator, Mr. Campbell, and caused the mer, 

_ Mr. Burke, to admit that he believed conditions were not of uniform flow. Mr. - 
Campbell noted that the one set of the four that gave the approximate trend - 

: of the logarithmic smooth boundary line fell below it. A very interesting fea- _ 
ture of the reference is the plot of observed velocity distribution in its Fig. 7, 7 


This plot deviates from the so- -called Prandtl one (that was not for pipes) in. . 


(herein), buttoagreaterextent. 
‘Finally, re references 35 and 36 | are very informative but do not seem to add 
eniiieenie,. for p purposes of the | present discussion, to the discussion of the 
paper of Hickox, Peturka, and Elder. They include Raam’s method of esti- 
: mating f for rock tunnels—a method that the author likes because it defines 
a practical method of measuring a parameter obviously related to roughness 
and correlates it with f. This procedure ought to make the f-n diagram, or oa 
any other, useable outside the range where roughness height becomes large 
with diameter and is. On a par, for practical with ex- 


3 the 1e inner 90% o of the pipe in the manner of the hyperbolic cosine curve of Fig. . mS 
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“treat canals. . The scope of the material covered is the widest yet pub- - 
lished, and their review of the data collected by others and by themselves is 
‘ most valuable contribution to the study of the behavior of man-made alluvial — 


"correlations ‘studied relatively well. as was pointed o out, he larger 


TABLE 8.—-COMPUTED MANNING 


Discharge, cfi 


canals, conveying flows in excess of 1000 cfs, seem to depart markedly from - 
_ ~~ general trends indicated by some of the correlations. In particular, Fig. 8 
shows a limitation of the mean velocity to approximately 3 fps in the Punjab 
and Sind canals when R2s becomes larger than 0.005. 


computed by application of ‘the Manning formula to the Punjab canals, one 

_ obtains values substantially higher for the larger canals than for the smaller — 

canals. For comparison purposes, an average value of n for the smaller 
_ The Blench-King formula, Eqs. 15 and 29, is also inaccurate. Table 9 shows 7 
_ the results of the application of the Blench-King slope formula and the Lacy 


1960, by D. Simons M. L. Alherteon (Proc. Paper 2484), 
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WATER CONVEYANCE CHANNELS IN ALLUVIAL MATERIAL4 
i MARCEL BITOUN,!¥ M. ASCE.—The writer was interested in this paper _ 
| 
data reveal that 
— 
4 
4 
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It — ‘that in of th these formulas slopes s sub- 
Although the authors indicate that a reasonable estimate of the an 
tg slope can be obtained by using the correlation of V and R2S or by using the 
< Blench- King ‘regime slope formula, this accuracy does not appear to be suf- 
ficient for design purposes, as the determination of the slope of irrigation | : 
canals is one of the most criticalitems. = j= 
A reliable formula for flow on movab 
_ sideration to the grain size of the bed material and to the bed configuration. > 
_ Such a formula was proposed by H. K. Liu, M. ASCE and S. Y. Hwang, A. M. 

- ASCE (69). The validity of its form was verified on the Punjab canals, and 
_computed slopes for the large canals are shown in Table 10. 
= The coefficients Cp listed in Table 10 are slightly aitenent: from the co-— 
efficients recommended by Liu and Hwang. They were derived by the 


‘TABLE 9. -COMPARISON. OF REGIME SLOPE FORMULAE 
= 


_ Measured Slope 

— 


. 
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0.2 
0.2 
0.2 
0.2 
0.1 
0.1 
0.1 
0.1 


0600 


adjustments to Se for apy application of the formula to different areas is attributed 

by the writer to the fact that the mean grain size is not sufficient to com- 

pletely define a bed material 


veg, 


definitely better than that of the ene by 
¥ the: authors, It is within acceptable limits for design purposes and it is “e- 
that its use be considered. 
_ Another anomaly is the fact that the Punjab canals, reputed to be stable, 
if a are characterized by widely different tractive forces on the bed and — 
different width-depths ratios. According to Lindley, Lacey, and Blench, for 
canals constructed in similar soils, , the -width- 
ratio and the velocity is 
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DISCUSSION 


correlation of W/D and Vv ‘is shown in Fig. 170). The width used i in 
correlation is as defined by Blench: W= A/D. 
It is apparent from Fig. 17(a) that the various canals are characterized by 
widely different values of K. For example, i in ‘comparing large canals No. : 
3, 7, 8, and 9, the bed material mean diameter varies only from 0.41 mm to — 
—60. 43 mm, and V varies from 2,54 to 3.25 fps, whereas K varies from 4.75 to 
4. 75. For the ensemble of the Punjab canals, K is comprised within an even F 


_ However a comparison of groups of canals operating under similar aa 


eter 

- oe ‘the same mean diameter of bed material, the tractive force on the 
bed “inereases | up to a certain value as the width- depth ratio decreases, 
2. The maximum value of the tractive force exhibited by the canals and 

the mean diameter of the bed material both increase as the width- depth ratio. zo 


> ‘TABLE 10.—SLOPES BY THE LIU-HWANG 
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=O 
opes computed by the -King formulas as on Table 9, 


opes computed by the Lacey formula as shown on Table 9, baer 


Table 11 shows the values of 450 w/D for the group 
of ‘canals operating at velocities of approximately 1.4fps. 
Tt is suggested that some of the canals could have been constructed with 7 
4 lower value of W/D without danger of bed scouring. For example, canal No. + 
3, in which the tractive force is 0.026 psf could have been operated with a 7 
. force as high as that of canal No. 14 (0.044 psf) without being eroded, 
_ Thus, , instead of searching for relationships common to all the canals a 7 
—tavestigaton , envelope curves corresponding to the limit of stability should 
be defined. Such curves have been drawn on Fig. 17(). They represent the 
“minimum width- depth ratios that characterize the stable Punjab canals and 
should correspond to tractive forces on the bed close to 7“ limit of — 


Analysic of the Diniah 
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Canal No|Material | Ch | x | |Slope(per|Slope (per | King? 
| o41 | 18:3 | 0.49 | 0.325] 022 
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FIG, 17.—PUNJAB CANALS 
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‘compatible a. mean grain size of bed “material. These curves are 
_ gented here only in order to illustrate the concept of limiting conditions. 
| Further investigations are necessary to determine accurate numerical values. _ 
The writer grouped the forty-two Punjab canals according to their pee y. 
_ median grain sizes, then selected the highest tractive force value for each 7 


group and t tabulated these” values Opposite their respective grain sizes 


0.088 


0.052 


0.087 


(2) 
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values of T canbe compared to the ones recommended by E. w. Lane 
stable canals. Fair agreement is shown by Fig. 18. 
or was previously mentioned that, for similar mean velocities, the wate ma- 


4 Group 1 corresponds to d 
Group 2 corresponds to d 


Group 3 corresponds to d5 


Fig. ‘17(b). Considerations of fluvial dynamics indicate that when balanced re-_ 
gime is established , the material that forms the bed of an alluvial channel is: 
_ characterized by ¢ a coarser mean diameter than the original soil in which the 


2 decreases. This” is put in . evidence by the family o of ‘envelope « curves ; shown on 
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" = suspension by turbulence and carried away on sail load. If the material in 
7 a place contains a sufficient proportion of grains witha size ‘such that they can- 


be removed, an armor of a certain thickness develops on the bed and pre- 
ae . vents aneed segregation. The median size of the armor increases with the 7 


‘Lane- Recommended values for canals 

high content of fine 


a) With clear 
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forces acting on the bed. ‘However, although this bed is not being eroded in a - 
‘ stable alluvial canal, it is moving as bed load. Stability requires that the bed a 
ae movement of the canal bed be replaced by a — of s solids into the — 
canal at its headworks, in size and in quantity. 


iod that precedes establish- 
channel was excavated. During the transition per tedtoaprocessof 
of regime conditions, the material in place is subjecte 
t= 
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DISCUSSION 
4 
lia, the case of the Punjab canals, it could be that if some of the wider canals _ 
(with relation to depth) had been constructed somewhat narrower, their bed 
material would have been segregated to a larger size without causing further = 
instability. They | did not silt ey because the solid charge was not ex-— 


_ It would thus be possible to predict the extent to which segregation will 


‘occur in ina new canal on the basis of the tractive 1 forces imposed, the gradation »? 
- the m material in place, and the solid charge p permitted to enter the canal at 

the headworks (or imposed by the parent river when sand control devices are 

not present). The design c can be conducted so as to allow segregation affecting 


7 


° 


sufficient amount in the material originally in place. This would lead to lower | 
_ Width- depth ratios. However, the width is limited in practice by the fact that > 
erosion of the banks could then occur. Stability of the banks could be deter- y 
_ mined by relating the tractive force on the banks to the friction angle of the = 
= in place and the slope of the cut (orfill), 
mene stability cannot be insured unless the bed load movement is de- 4 
‘ee for the chosen size of canal and segregated bed material size, and oi 
aan the supply of sediments into the canal at the headworks is maintained | 


- equal in size andin quantity to the material migrating downstream | as bed load. . 


— os bed material up to a large grain size, if such grains are available in . 


variation of sediment concentrations should be correlated to the variation of — 
aoe the canal imposed by irrigation requirements. An integration 4 
of the aggradation and degradation a provide a yearly cy cycle resulting in 
In that respect, it would have been diaeiitine to find among the data pre-_ 
sented detailed information on seasonal variations of flow in the canals in- 
vestigated and on the sediment allowed to enter the canals. 


Com _An approach to the regime theory should not be based on a statistical 


= the year, stability of the canal must take this into account. The 


—— 


— ~ 


SS results” ‘equivalent to those ‘obtained from 3 an .n approach based on the tractive 7 
force theory. These approaches, however, should be supplemented by due 7 
consideration given to solid transportation as one of the determining factors 
of stability. _ Addition of a bed load formula to the basic equations of stability _ 

is necessary and should enable a complete determination of the regime all 


equilibrium. stability implies a an integration of allowable ag- 


water requirements. 


69, “Discharge Formula for Straight Alluvial Channels, H. Liu and 


Ss. Y. Hwang, Proceedings, ASCE, November, 1959. 


Mi MONIR M. KANSOH. 11_This study of existing stable canals in the United 
States is actually a continuation of the regime studies: that: were initiated 


— 


| 

— 
— 

— 

a 
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less than 0.20 mm with a concentration varying from about 3500 ppm 
on to about 200 ppm with an average of about 1600 ppm, nearly similar in 
, conditions in the Sind in which the “mean size of bed material is within - 
_— the limits 0.0346 mm to 0. gues mm” and “the suspended sediment load ranges 
from 3590 ppm to 156 ppm.” ” This isto be compared with the conditions in the i 
Pale in which “the average diameter of the bed material is approximately — 
4 -43 mm” and the mean silt intensity | is, like that in the canals investigated _ Zz 
by Simons and Bender, about 328 ppm. - According to the authors, the tractive 
force theory may not be applied in Egypt where water is not clear and where : 
es the bed and sides of the canals are, , in most cases , composed of fine cohesive | 
materials. Design of stable canals should, then, depend on regime formulas. 


= Egyptian canals, because the silt is finer than that in the Punjab, the © 

velocity is sometimes _ taken as of the value given by Kennedy 


using units, 


cially for canals in Lower that the formula 
be used in v which m = uke whereas | p varies iiditin 0.26 (for fine stilt) to se ol 
(for coarse silt), using metric units. 
‘The form of the canal cross section is neniaed to be inane: 
it is, as is usually the case, a trapezoidal cross section with a bed width B, 
‘D and wi with at an 6 with the horizontal, then 


= Q/(p DO cosee 
p-D™// 


Va To fix ‘the: dimensions of a stable canal in order to convey a discharge Q, 
the critical velocity | as given by the empirical regime formula is substituted 


ow. 
the Punjab by Kennedy and that yielded the empirical 
— 
] 
: _ Observations on thirteen canals in different parts of Egypt by K. Ghaleb be 
45 
— 


into the Chezy equation with the aid of one of ‘the fc formulas given by Bazin, 
Manning, Kutter, and others. The Manning formula is often 
because, in 1 addition to its simplicity, it gives satisfactory results. 


(foot units) or 


The depth D can be ‘computed after p putting = in terms ¢ i D as aa 
4 The bed width can be directly 
and even avoided by using diagrams, 
_ relation between BD and D, or between W W and D might a: as well be ~ : 
tained from « empirical r regime "equations like Eq. . 4, or Eqs. 13 and 14. 
After a_ careful examination of a large number of recognized canals in a 
Egypt, Molseworth Yennidunia recommended the the relations 
) 
_ for D _— 1.62 m and below 1.62 m respectively, § Ss ae the water slope in 
canal in centimeters per kilometer. 
_ They also recommended the ratios | 


7 


"when BD is greater than 2. 00 m, in which C1= = 1. -00 and Cg = 1.75 for ‘deep | 


drains while C; = 0.90 and C2 = =1, 45 for shallow drains. 
These relations are vepresemes in Fig. 19, that gives the desirable ‘pro-— 
portions of of ‘stable trapezoidal canals” ane drains with L :1 side slopes. The 
‘dotted ‘curves that correspond to the value 
2. The ratio-curves corresponding to canals of different nt slopes, 
and shallow drains. 


‘The point of intersection of the appropriate K-curve and the proper ration 
the ‘desirable e depth ai and t bed width of 


& 
— 

— 

a 
when Bis less than2.00 mand 
— | 
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FIG, 19. OF STABLE TRAPEZOIDAL CANALS AN AND 
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DISCUSSION" 


Interpolation is sometimes necessary when the computed value a 

lies between two K-curves, or when the slope lies between two of the slopes 


q 
indicated by the ratio-curves. 
The same relations of B and D, for are the basis of a standard 


of stable ‘trapezoidal canals with on- 1-1 side slopes when the 
solids in suspension do not exceed 2000 ppm. This diagram, Fig. 20, is suit- 
able for one class of canals is whose Manning roughness = 0. 025, 
_ To use the standard diagram, Fig. 20, a vertical straight line indicating 
the is extended to intersect the right hand inclined straight line" 


oa te left hand ‘inclined straight line corresponding to the same selected 
a point that directly indicates the desirable depth and bed a 
width ofthecama, 
_ The procedure may sometimes be altered by fixing the depth and trying to 
“tind slope and bed width corresponding to the given or required discharge. It 
may be the bed width that is first known while the depth and the a are to 


determined, 
_ Hesitation usually accompanies the selection of the ‘suitable pe indi- 
cated by one » pair out of the two groups of inclined straight lines. To reduce 
this hesitation, HL El Defrawy suggested a particular selected slope for each — 
range of discharge, depth and bed width. These e slopes are shown in thick lines” 
_ A. M. El Banna proposed a diagram, Fig. 21, that is a combination of the 
standard | diagram, Fig. 20, and an old nomogram known as Wallace diagram. 
-- For using Fig. 21, it is sufficient to extend the straight line in both direc- 
tions joining the value of the discharge Q with the equal number on the scale 
line. The straight line will directly give the desirable depth, bed width, ond 
- slope. In choosing the number on the scale- line, a little adjustment is is some-— 
times necessary so that the bed width may bea round number. bau 


_ Adopting the dimensions given by this diagram, permits a velocity of flow 


ad in which 
m= = 1,00 and D is the depth of flow when the canal conveys the maximum de- 
_ The authors state that Q could be > correlated with | D but, for design pur- 
_ ‘poses, R is a more meaningful measure. This may explain why Q is given in. 
terms of P and R in Figs. 3 and 5. In fact, neither P, R nor the top width is 
of a practical value in fixing the dimensions or proportions of a canal cross 
as is the ratio B/D or w/D, and so on. 
The many diagrams given in the paper are, in the writer’s opinion, not 
more useful than any one of Figs. 19, 20 or 21. The former diagrams could — 
_ preferably have been amalgamated into a single diagram, nomogram or abac, © 
_ to give a direct solution for the problem of fixing the slope, dimensions and/or _ 
_ proportions of stable canals under different conditions. Sl 
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all Nevertheless, the information given a ‘t relation or proportion 
of the average width and the depth of stable canals. The width depth ratio can 


From Fig. 2 the nearly equals 1.12 W. Noting Fig. 3 an and Eqs. 


668 00. 50 
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FIG, 21.—PROPOSED DIAGRAM FOR THE DESIGN OF TRAPEZOIDAL NON 
CANALS IN EGYPT WITH 1 ON 1 SIDE SLOPES 


az Similarly, taking the values 0 of the hydraulic radius R from Fig. 1 (roughly 
0.81 D) and from Eq. 12, from Fig. 5, and Eq. 21, for sand bed ane natural — 
berm, and from Eq. 22, for coarse that 


Hence, required ratio may be taken | as 


for sand bed and nd natural berm, or “bi 


- 
— | 

| 
— 4) 

— 

(85) 
— 
/D=7.8Q°” ... 


— 
_ The last kind of material was mentioned, with this description, many times — 


_ in the paper. The writer wonders if there exists a coarse cohesive material. 


One of the important points is the limitation of the validity and applicability - 
of regime equations that, as was stated, are valid only within the range of the gal 
Se data or are only valid for the limited range of conditions on which ~ i 


a they are based. Two of these conditions were mentioned, in the same page, q +4 


whereas one, a very important one with regard to the range of applicable 
slopes, was ignored. It is true that the range of slopes was mentioned as 
varying from 0.000058 to 0.000388, but this information is most. confusing and 


misleading as explained by the following example. 
____ Considering a uniform flow in a stable trapezoidal canal of a bed width | B, % 7 a r) 
depth D and sides’ sloping ; at anangle @ to the horizontal, and combining Chézy 


Meaning formulas with the Kennedy in order to facilitate 
computation, may be put inthe form 


=p 
instead of the 0 ong form of Eqs. .72. Then 


(B+D¢ ot 6) D 
| 
90) 


(91) 


(cot 6) - - - 2 (cosec 
n is” given an average value 0.025 and if 03 is considered to be 45° and, ; 
supposing that the canal conveys water fine such that Ke 
equation can be less than 


1,486 
0.396 x 0.028 


15008°/* - 2.828 (3b) 
= 
1500 s°/ 4 0.000665 - 


- (0. 


we 


5, 
— 
| 
| &§ 
> 
or. 
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“that int case met water fine silt 
uniformly — in a canal of n = 0,025 and @ = 45°, the slope S should vary from 
0. 000058 to 0.000232 only, whereas the slope is stated, in page 66, to vary 
_ from 0.000058 to 0.000389 without limitation or reference to any conditions, — 
Even in the of slope 0. and 


0, 000232, the ratio B ia drops quickly f from infinity to to zero. ro, Reasonable v values 


a of pmyt then be expected only w within the the range of slopes f from about 0.00010 — 


to about 0. 0. 00020, f for example, under the described circumstances, a 
~— In addition to the care necessary to design stable canals, measures must 
- be is taken to reduce silt | troubles, for example, by feeding the canals and 
os branches from the outer side of curved reaches of rivers and main canals, 
well as providing the with sand 


e 


om em wg 


| Real values of B and D necessitate that the ratio = should be ositive that 

| 
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DRAG FORCES IN VELOCITY 


— 


al authors at the planning stage of the investigation and designed and tested the 
_- force-measuring device that is described briefly in the paper under Mr. - 
_ Moore’s direction, for a previous investigation. ” One of the primary aims in | 
the development of this device was to make it adaptable for use in subsequent — 
_ investigations, and it is gratifying to the writer that apparently this aim was» 
‘The a authors very carefully point out that this was only an exploratory 
oe. presumably because they felt that the limited amount of data did 
not permit more than general conclusions to be reached. There are, , however 
number of points that could perhaps be given additional attention. 
Eq. 4, in which, presumably, p is pressure intensity in pounds per square 
inch and r the cylinder radius in inches, gives the local force, f, in pounds 


& per inch of height. For use in Eq. ‘3 however, -f must be in pounds to give a 


_ dimensionless Cp. It would seem likely that the authors actually multiplied the © 
sade from Eq. 4 by some unit of height, possibly the orifice spacing, to 


determine the local force to be used in ‘Eq. _ 3. Could the a authors clarify this 


on the writer’s knowledge, the variable shutter gate used for achieving the 
- desired velocity gradients has not been reported previously in the literature. | 
Could the authors provide more details on the construction and operation - 
the gate and the _ range of effects achieved? How does it affect transverse “- 
aa profiles? Is it possible to achieve a negative velocity gradient, that — 
is with the maximum velocity at the channel bottom? Where were the velocity 
profiles: of Fig. 4 taken with respect to the gate and the measuring cylinder? | 
| For how great a distance downstream did the gradient persist before boundary | 
layer fo formation produced a more normal velocity distribution? = —— ry 
a The statement is made that “the measured total drag included, of course, 

_ the wave drag at the free surface. An approximate evaluation of the wave drag 
was made, based on the profile of the wave surface around the cylinder.” 
Could the authors give more details concerning the procedure for evaluating © 
_ wave | drag and, in addition, some quantitative indication of how “satisfactory” — 
the agreement was between measured and computed total drag? The total ca 


= for the various runs would be of interest. 


a a july, 1960, Frank D. Masch and Walter L, Moore (Proc, Paper 2546), 
6 Hydr, Engr., Turner and Collie Cons. Engrs., Inc., Houston, Tex, 
> Tan Experimental Investigation of the Drag Characteristics of a Model H- Section — 
Bridge Pier,” by Donald VanSickle, Thesis presented to the Univ, of Texas at Austin, 
in 1958, in partial fulfillment for degree of Master of Science, 
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cylinder Fic shows curves for only 10 of these with numbers 
4 q 


1 4, 7, and 14 through 24 being omitted. Similarly, on Fig. 5, ee local ar: 
‘coefficient is shown for eight points in one graph, ten points in another, ‘and 
thirteen points in the others. Because the depth was the same for all cases, 
” meetin the same number of points would have been available for all. 
Could the authors explain the missing piezometers? er 
ae In Fig. 4, the upper two curves indicate that velocity measurements were > 
taken at the water surface, because the depth was Stated to be 0.8 ft. Were 
readings actually taken at the surface, or was the depth for these curves” 
‘The label for ordinates in Figs. 4 and 5 is more correctly § given as “Dis- 
- from the Bottom, Feet” as described in the text. _ aoe 
cual Based on the stated flow depth of 0.8 ft and the average velocity of approxi- 
7 ry d 2 fps, the Frouge Number for the authors’ tests would be on the order 


the dash pot the authors find that these oscillations of 
7 the cylinder caused similar pressure oscillations on the manometer bank? If 
so, ‘was any difficulty encountered in taking readings on 13 (or 24) manometers- 

a In Fig. 1, the arrows indicate ‘that on the ur upstream side of the pier, at the _ 
_ channel bottom, the induced flow is inthe upstream direction. Was any attempt — 


made e to inject dye from the bottom piezometer to check this? = 
_ Had the authors considered the possibility that the surface wave could be 


"upstream stagnation point, and a “standing: depression” at the downstream 
ts stagnation point, the differential head between the stagnation points and the 
surrounding fluid could induce flows similar to those observed. 
a In Fig. 3, those portions of the curves lying between 90° and 270° are shown 
7 as being positive, Because the cosine inthis region is negative, the pressures © 
= = must also be negative to produce a positive product. Are these pressures 
x negative with respect to surrounding fluid pressure or some other datum and 
how was the “datum pressure” measured? Because the cosine between - - 
_ and 90°, and between 270° and 325° ' is positive, negative pressures must ob- 


: - tain between these pened to give the negative product shown. Do the authors a 


the downstream | stagnation points for cur ves 10 through 13 is at least twice 


- amples of engineering problems in which velocity gradients are of significance. 
. It should be pointed out that in both types of problems pressure fluctuations 
and flow accelerations are also present. While the authors are evidently aware 
Of this, they seem to give the impression that the velocity gradient is more 
important than | acceleration of the flow. Do they feel that their results could, 
- in general, _ be said to apply to flow with acceleration even though the experi- 
_ments were performed in a steady, non-uniform flow? | 
_ The writer certainly does not wish to detract in any way from the value 
v of the authors’ work. While he feels that knowledge of the variation in cg _ with 
velocity gradient is important, he also feels that knowledge of the <adieiion ; 
fe: in cg with acceleration of flow is of at least equal importance. In all probability 7 


both and flow acceleration are of significance, and must be 
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c 
= ‘cerning the effect of accelerations on the sattatanee of cylinders, this is the ‘ 4 


. first paper, to the writer’s knowledge, that is : concerned with the effect of ve- | 
locity gradients on resistance. For this reason be is a most welcome contri- 


Although ‘neither the aforementioned the authors’ paper give 
information that can be used directly by designers, they do illustrate the im- 7 

portance of velocity gradients and flow acceleration on drag forces. The writer 


wishes to congratulate | the authors for making t the 1 results of t their 


lated efforts to determine reliable design criteria. 
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tion,’ by Ss. Russell I 


ceedings, ASCE, Vol, 82, December, 1956, p. 1113. 


9 “Ocean Wave Forces on Circular Cylindrical Piles,” by as. Wiegel, K. E, Beebe, 
and James Moon, Proceedings, ASCE, Vol. 83, April, 1957, p. 1199, 

_ 10 “Water Forces on Accelerated Cylinders, ” by A. D. K. Laird, C. A, Johnson, and 
Ww. Walker, Proceedings, ASCE, Vol. 85, March, 1959, p. 99,0 
“Water Eddy Forces on Oscillating Cylinders,” by Alan D, K. Laird, Charles A. 
‘Johnson, - Robert W. Walker, Proceedings ASCE, Vol. 86, November, 1960 p, 43, — 
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PREDICTING RUNOFF | ON SMA SMALL EXPERIMENTAL WATERSHED@ 


ais _ JAIME AMOROCHO, 4 F. ASCE.—The author has presented ir 
that show the inherent limitations of the unit hydrograph approach a asa a method 
of predicting the outflow response of a watershed to individual inflow episodes. _ 
However, some of the interpretations given in the paper to explain the appar- 


ently irregular behavior of small watersheds in this regard deserve further 
In accordance with the unit hydrograph principles, the ordinates of any 
storm hydrograph are directly proportional to the ordinates of a unit hydro- 
graph of the same basic duration. The factor of proportionality is the effective — 
intensity of the storm expressed as a multiple of the unit rainfall excess. 
This may be said more formally by stating that the unit response to any net 
inflow episode is an invariant function of the physical characteristics of the 
catchment. The term “net inflow” is 1 meant to signify the portion of the gross : 
inflow that becomes surface runoff. The principle (of superposition implicit 
in this general ‘assumption is extended to any time sequence of inflow rates — 
> and precludes the possibility of using different unit response functions (unit i= 
hydrographs) as the inflow rate changes. = = © 
yy, In very few if any of the field studies performed to date can it be claimed 
that the unit hydrograph theory has been tested strictly. This is because in : 
practically no instance has it been possible in nature to ascertain with any 
degree of assurance that a storm pattern that caused a flood shydrograph — 
actually measured, was ‘proportional in all its elements to the storm used to 
_ derive the unit hydrograph under analysis. Other factors, such as the uncer- 
_ tainty i in the determination of the infiltration r rates 3 used fo for the the construction 
‘of the ‘unit hydrograph < and for the predicted flood hydrograph ‘contribute v very, 
importantly to the weakness of such tests. Hence, the fact that the approximate = 
outlines of recorded flood 1 hydrographs |! have been reproduced in a number of 8, 
instances on the basis of a unit graph previously | derived may appear to be 
largely fortuitous. For every instance of approximate agreement there have | 
= very many instances of gross discrepancy. It is a well-known fact that 2 


the unit graphs derived from large floods usually differ from those derived — 


‘The watersheds investigated by the author are, by virtue of their small 


areal extent, particularly well suited for a critical check of the application 
7 a the unit hydrograph concepts. As a direct result of the reduced areas that ; 


make the occurrence of very complex geometric patterns of rainfall with 4 
steep, non-stationary intensity gradients less likely within the catchment — 


_ boundaries, the conditions required for close overall proportionality of in- 
flow can be expected to become more possible. ase tel watersheds normally | 


4 August, 1960 by Neal E, Minshall (Proc, Paper 2577). 


Asst, Prof.,, , Univ. of California, Berkeley, Calif, 
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: = ferent locations, preter no major discontinuities exist in the storm ila 
b>. ™. over the catchment. In these cases the outflow hydrograph becomes rather 
= insensitive to inflow intensity changes recorded at individual precipitation 
| stations , and its shape tends to reflect the catchment characteristics only. 
" However, because major rainfall discontinuities are more likely in extensive - 
watersheds » the conditions conducive to of inflow are, 
The author has derived ‘unit on the of floods 
w short storms of different, “reasonably uniform” intensities. The hydro- — 
4 eras shown in Fig. 6 were presumably reduced to the same unit es 
fe : excess duration in order to render them consistent with each other. This is 
4 Za not indicated in the paper, but because the observed rainfall excess durations . 
an as reported in Table 3 are approximately the same, no large error is to be — 


M 


expected if the conversion was omitted. _ The striking fact is that the the hydro- > 


represent « distinctly different outflow responses. 

_ A basic premise of the unit hydrograph procedure is that the ‘relationships | 

between storage and outflow are linear.5 This results from the requirements 
for a linear solution with respect tothe net inflow rate of the storage equation — 
fora hydrologic unit. To illustrate, consider abasin receiving a constant rate , 


‘The equation c can be stated as 


rate of change in This equation has an immediate 
ind provided S and Q are related linearly. For if we can write 


that 


this equation, the a= in the exponent is a characteristic 
parameter of the basin. The constant K has the dimension of time and repre-_ 
sents the lag due to a “linear reservoir” as defined by Eq. 3. The Ps 


a(t- t 


The subindex refers to the v enete at the moment of cessation | of the | inflow. 
Now, because by Eq. 5 the discharge at peak is — 


7 «al is seen Eqs. 5 and 8 are linear of I, that is the basic 
of the unit hydrograph method, J. C. I. Dooge,5 M. ASCE, has 


Journal of Ge of Geopk h nysical Research, J. C. I. Dooge, February, 1959, 
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DISCUSSION 
an exhaustive analysis of the unit ‘theory in which variability 
In natural watersheds, the assumptions of linear reservoirs and linear 
. ae processes can not be taken for granted. The terms S and Q are 
normally related by non-linear functions which contain additional variables 
an as friction and velocity gradients. Introducing also the variability 
f the inflow rate, the storage equation has the 
‘Tt can ‘be easily surmised that th this equation will, in general, not have a 
c = Practically all cases in nature are of this complexity; therefore, it may y be 
a ; stated quite generally that the classical unit hydrograph concept can almost | 
never be applied strictly. It should be emphasized that the preceding conclusion — 
governs not only the proportionality between the ordinates of the runoff hydro- 
graph and the unit hydrograph but the additive superposition of the effects of — 
Some rather interesting conslusions may be derived from Figs. 5 and . 
4 that illustrate the previous statements. The derivation of unit hydrographs ~ 
from storms of varying intensities give rise to inaccuracies that depend on yD 
7 ‘Subjective factors that inevitably creep into the process of separation of the © i 
different component runoff episodes. Inaccuracies also result, as mentioned 
previously, from the determination of the infiltration rates, that, in the pres-— : 
% ent case, was made on the basis of a completely empirical criterion, Never-— S 
~~ “4 theless, some measure of verification of the latter | was afforded by the ob- 
. ‘served actual runoff, and the clear trends apparent in the alignment of points : 
plotted in Fig. 5 appear reliable. From these trends it may be concluded (1) | 
4 that the Shape of the unit graph is a non-linear function of the inflow rate; and 
me) that - depending on whether the unit graph corresponds to an initial period ~~ 
«of inflow or to inflows that occur after some runoff created by previous rain _ ’ 
is taking place, the response is different although the intensity may be the 2 
same. This may be interpreted as a “field proof” of the pepe yi of the — 
principle of superposition in this case BA ag 
It is obvious that the exact relationships suggested by the the Fig 5 
apply only t to the particular watershed described inthe paper. Other 
_ ships will exist in different catchments, but as long as their empirical deter-_ 
: . if mination is possible, the type of analysis that follows should have general ap- 
‘ plication. The only requirement is that the overall geometry of the inflow _— 
@ tern should not vary greatly within natin basin boundaries for different rain — 


cross-plotted to define a direct between the peak and 
the time to peak. Scaling values from the figure, it is found that the following 4 


q 


i It: 25, 


For unit response early in storm: = 15.25t, 


For unit response late in storm: pinta 
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4 . time to peak. These equations represent the loci of the peaks of hydrographs — ai 7 


of curves that must satisty the following conditions: 


. To m xima passing he’ lines by equations 


To encompass an area of unity between each curve and the time axis 
the definition of the unit hydrograph). 


— * oe To eels value at time equal zero and to be asymptotic to the time 


— be continuous and have a positive second near the 


anda derivative equal to sero atthe origin, 
__ These conditions are met for example by some asymmetric frequency dis- 

ns It may be noted by reference to Fig. 6 that the curves tend to be-— 

come symmetrical as the time topeak increases. We might, therefore, choose 
a position in which this symmetry is closely approximated and fit therein a — 
normal frequency distribution to represent the response, The 
at time zero would have a finite value, but the error thus introduced in meet-_ 

= conditions 2 and 4 would probably be small. This basic function aie 


then be e skewed a and peaked simultaneously by introducing appropriate ool 
_ meters: in such a way that Eqs. 10 or 11 would be satisfied for each value _ 
a tp. A similar approach could conceivably be used with other functions, in 
_ order to arrive at some empirical expression fo for the unit response as a func- 7 
__-It is desirable, however, to use a function that is ‘meaningful in terms of . 
theoretical ‘analysis of watersheds. J. E. Nash® and Dooge® derived inde- 
the general expression of the outflow function for a a catchment wit 
_ pression has the form of a Poisson distribution and represents the unit hy- r. 
drograph of a linear catchment of this 
n whi == m 


4 
= ‘The valies of a and n can be computed from observed outflow hydrographs | 
& by the method of moments as described by Nash, A simple procedure, that 


_ assures that the maximum value of q satisfies Eqs 10 or 11 ‘exactly, although 7 
the fit may be less precise elsewhere, consists of finding the solutions for a i 


and n n from the formed 12 and 


‘obtained from the data of Fig. 5. The effect of the imperfect fit of the din 
to the unit response derived from the data is of little practical consequence 
a as willbe discussed later, = 
For q maximum, the Eq. must zero: 
ily ong 


«6 Proceedings, , J. E, Nash, hae of Toronto, Internatl. 
Hydrolo Se tember, 1957, 
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and combining Eqs. 16, 12, 13 on 14 
1) 


An explicit expression for n in terms of | experimental parameters 5, 


For n n< 9 it is that 


(n- 


Therefore, Eq. 20 becomes 


Likewise, from Eqs. 13 and 14 we get 


From the e data presented in the paper we find that 6 = = ‘* 170; * = = 0. wacal 


2.0090 1°: 326 | 1. 


= = 4. 1569 4.15691" +0 3518 eee (25) 

Fig. an shows a oumaorton between the hydrographs presented by the 

- author in Fig. 6 and the curves computed by means of Eqs. 12, 24 and 25. 
Bee It is noted that the fit becomes increasingly better for higher values of the 7 
; inflow intensity. Similar fits computed by the method of moments were found _ 
to be slightly better for the curves corresponding to low intensities but be-— . 

_ came increasingly difficult at higher intensities. The precision of the fit is 

within the order of magnitude of the error to be expected from the field data, 

4 as may be surmised from the scatter of the points shown i in Fig. 7. Ih general, 
4 slightly higher runoff volumes would correspond to the areas below | the era 


ay curves during the first hour than to the areas below the curves shown in Fig. 
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dilaatuoes by functions that define the number of equal theoretical linear ee 
_ reservoirs and the theoretical lag time of one linear reservoir in terms of 
the inflow rate. The substitution of these expressions in Eq. 12 gives an equa- a 
- tion that can be considered as the approximate solution of Eq. 9 when I does | 
ya not vary in time. An interesting implication of these equations is that a water- 
shed may be regarded as equivalent to a linear catchment only when the in- 
a flow is constant. As soon asthe inflow rate changes, the number of theoretical — 
reservoirs becomes larger or ‘smaller and the response varies accordingly. 
‘Therefore, neither n nor y may considered as invariant — 


| 
= 2.65 in 


he | FROM FIELD DATA (Figs) 
—— COMPUTED BY Eq. 


FROM BEGIN OF EXCESS RAINFALL, (minutes) 
FIG, 11,—COMPARISON OF UNIT RESPONSE CURVES FOR DIFFERENT 
INFLOW INTENSITIES 


parameters of a catchment. This conclusion explains the inapplicability of the 
principle of superposition in this type of hydrologic problem, and may not _ 

only be extended to the critical evaluation of the unit hydrograph method but % _ 
other related flood procedures, such as method, 
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preceding analysis suggests a “practical procedure e for obtaining the 


function for small watersheds in which uniform 


— 


4 to periods of reasonably uniform rainfall excess. The values of peak discharge _ 
and time to peak should then be plotted against rainfall intensity and from | 
these - relationships, by an n analysis similar to that described in this discussion. = 
the expressions for n and a in terms of I could be derived. It should be 5 ll 
in mind that whatever form these expressions may acquire, their nature is 7a 

a statistical and empirical. ‘Hence, from a practical standpoint there is _ 


no object in scrutinizing the physical meaning c or the dimensions of the para- | 

‘ 


Berkeley, Calif., these ci concepts are » being studied in model c: 
be - ments. A full theoretical and experimental treatment of this subject, initiated — 
ini 59, is being carried out under a of the the University’s s Re- 


recorded. This procedure would consist of first selecting a number of 
and 4: > Lan 


_ PREDICTING STORM RUNOFF ON SMALL EXPERIMENTAL 4 


ate D. Dougal 
D, DOUGAL,” A. A. M. ASC CE.- Research and have 
indicated the many variables that must be considered in developing a usable i 
. rainfall-runoff relation. | _ Many that have been developed present unique prob- © 
- lems in their use in estimating runoff on ungaged watersheds. The author’s | 
- work extends known techniques both in methods of approach and in providing © 
information on a particular area of the nation. akeoacnicemenaar 
_ ‘The writer has studied® the problem of a modified infiltration index and the ? 
rate of infiltration rather than either total retention or runoff estimates, but — 
3 is apparent that either of the latter could be obtained indirectly. The ante-— 


cedent precipitation index approach was used, and season (or week-of-year) - 
_ duration of excess rainfall and average intensity of rainfall were the other _ 
independent variables given primary consideration. Certainly the Ralston 
‘Crem watershed near Iowa City, Iowa is classical in hydrologic research. 
_ The use of the antecedent precipitation index, computed by Eq. 1, appears an 
.: be widely applicable. As recommended by Kohler and Linsley,? the index ag 
should be used in conjunction with season of the year or temperature. The > 
a - author states, for the various methods and lengths of period tested in com-_ 
ies the index, that the best correlation with retention was a 30-day index 
plus 20% of the departure from normal of the 30 to 60 day precipitation. It is 
later stated that an antecedent retention index was used in an attempt for rre- 
= finement, but the results showed no noticeable improvement over the 30- day 
_ index finally selected for use. The recession factor “k” in Eq. _1 was evaluated 3 
as 0.95 for the Edwardsville, Illinois watersheds 
‘The selection of 0.95 for the recession factor, a a relatively high value, “—— 
_ dicates that previous precipitation affects soil moisture in the watershed to a _ 
"substantial and somewhat lasting degree, because the API values recess quite 7 
slowly. This means that if appreciable precipitation has previously occurred 
its effects on soil moisture linger on. Perhaps the author might comment on > 
_ the effect of the impervious clay- ‘pan layer on the relatively low amount of soil 
_ The writer initially investigated three values of the recession yn factor “k” 4 
for the Ralston Creek watershed: namely, 0.84, 0.90, and 0.94. In the result- ~ 
ing graphical correlation, the values of 0.84 and 0.90 were about equally ca- 4 
_ pable in enabling ‘first approximation month- of-year curves to be | drawn, but 


7 Staff Engr., Iowa Natural Resources Council, Des Moines, Iowa, 
8 “A Study of Factors Which Affect Infiltration Rates on the Ralston Creek W. ater-_ 
- shed,” by Merwin D, Dougal. Thesis presented to the Graduate | College of the State 
Beene J of Iowa, at Iowa City, Iowa in February, 1958, in partial fulfillment of the re- “a 
ments for the degree of Master of Science, 


-edicting t the Runoff from ‘Storm Rainfall,” by M, A. Kohler and R. K, Linsley, 


Paper No, 34, Dept. of Commesos, Weather Bureau, D. C., 
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March, 1961 
little correlation was , evident for the. value of 0. 94. The mean of the first two oS 
0.87 was subsequently selected as most applicable for the Ralston 
_ Creek watershed. This watershed is about 3.0 square miles in area, relatively 
ie rolling to hilly. The soil is of loessial origin over glacial till with a Clinton 
silt loam soil surface . The Ralston Creek watershed experiences high infil- a 
In regard to the methods and lengths of period available in computing the 
wae precipitation index, the writer would like to add several others. 7 


fore a . storm: ‘in which to recess the initial index. The author does not state 
initial value was used in his study. The — | four proced- 


a Beginning w with a zero value of API 60 days prior to the storm. 
a. (2. - Beginning with the normal 10-day precipitation for the period prior t to 
| _ $0 days before the storm (first storm on May 1, use as the initial value the 
10-day normal prior to April 1, and begin API computations 
Aprill). 
Beginning with the observed average 10-day precipitation for the 30 
period prior to 30 days before the storm (first storm on May 1, use as initial 
_ index the 10-day average precipitation of March and begin API computations — 7 
: Same procedure as (3), except the 10-day average precipitation is — 


to a 3-week period before the storm. 
c= The four procedures, when applied in the Ralston Creek study to periods © 


prior to storms indicated below, yielded the values indicated in Table 5. The 
results indicate that the initial value need only be reasonably assumed. To 

illustrate, for the storm of August 4, 1924, the initial value for procedure (2) 
Was 58 in., the: day normal for June. _For procedure (3) the observed 
-day average for June was 2. 47in. As indicated previously, the 

respective API values after a 30-day recession period were almost identical. 
E It seems possible to the writer that the use by the author of a different “k” 


| factor ina 30- ~day index might have given identical results to -™ 30- 0-day index 
plus the 20% departure used inthe study. 


_ A comment by the author on his method used in comune g API values for — 
‘successive storms within a day would be appreciated. For instance, in Table 
1, the same API value is listedfor three storms occurring on August 14, 1946. 
Presumably the API value is higher for successive storms on the same cay >, 
by direct addition, or by recessing on an hourly basis. 
_ The writer concluded, in the analysis of API in the Ralston Creek diate, 
= that it is simpler tocarry the API computations forward continuously from S 
first storm in a particular year to the last storm. Less computations are 
a involved than if separate API computations are made for each storm. In addi- 
a ™, , this continuous method would be used in flood forecasting techniques. 
The writer agrees that a large volume of data is needed to develop graphi- _— 
ler ‘- ‘correlations in a hydrologic analysis. Over 258 storms in a 1927 to 1956 
record were used by the writer and all were utilized to advantage. However, 
; _ about 100 of the storms were selected as having “ideal” storm characteristics 


7 and were used in initial correlation efforts of API and season, month- -of-year, _ 
_ The author indicates the variation in season to season that can be expected. — 
the final grouping some periods to periods, it is 
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"DISCUSSION 


_ ation in plotting the data. The writer attemptedto use temperature asa para- 
_ meter (departure from the mean weekly temperature for the week prior to the 
storms under consideration). It was concluded that if all other variables were 
- about constant, temperature effects were of importance, but normally were 
too obscured by | other variables to be of value. 
‘The writer believes that the abscissa in Fig. 3 could be more aptly termed — 
_ “Duration of excess rainfall, in hours”. because the plotted points indicate Ps 
_ single total storm event and not the variationin retention during a storm. Th 
author also assumed that infiltration capacity is independent of rainfall in- _ 
+ This is a perplexing problem to analyze because in cultivated areas 
the intensity may effect surface cover conditions, especially in the spring — 
months. The comparison of estimated and observed runoff indicated in Fig. 40 
" cialis results for runoffs above one inch, and it is above this that = ; 


August 4,1924 
March 14, 1934 
May 3, 


28,1947 


‘oa API (k = 0,90) 


engineer is normally concerned in estimating storm runoff. For use in esti- 

mating past yields, individual estimates might be in error, but an a 1 
average should be obtained over a period of days or weeks. - aT Pe. 
_ The writer, in the Ralston Creek study, was unable to correlate certain 
_ gtorms that appeared on further analysis to possess a common characteristic, [ff 
that is, a fair amount of intense rainfall had occurred within a few hours a 


_ hr period maximum) before the beginning of the storm being scrutinized. This _ 4 
was finally sees to an additional variable, namely, the variation or dif- . 


easily reasoned that there is considerable. difference in soil surface conditions 
_ whether the API is falling or rising to the index existing immediately before 
occurrence of the storm rainfall. Toillustrate, assume an API value of 2.0 in. 
exists immediately prior to each of two hypothetical storms. For one storm : 
the API of 2.0 in. has resulted from the log arithmic recession of a 5 in. to 
‘ 6 in. rainfall a week or more before the storm in question. For the second 
hypothetical storm the API of 2.0 in. resulted primarily from a 2 in. intense 
rainfall a few hours before the storm in question. For these two hypothetical 
storms, all other conditions being the same, the obvious difference in soil 
moisture at the surface and in soil surface conditions will cause infiltration 
values for the two storms to be much different. The writer introduced shift- _ 
curves to account for this phenomena that enable, in effect, a shift to be made 
4 to a higher API value when a rising API situation was encountered, thus giving | 


interest to inquire of the author’s method for allowing for the seasonal vari-_ 
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The author has accomplished an admirable evaluation of hydrographs on 
small watersheds, The variation of the peak rate of runoff with both rainfall — 
—— and pattern of rainfall intensity is of interest. In regard to the vari- a 
ation with intensity, is it possible that the high intensity causes surface depths 4 
‘ = to reach a point at which not only is the flow velocity increased appreciably, E 
. but surface vegetation is altered; that is, flattened, disrupted, and so on? The © 
-author’s results, as indicated in Table 4, reveal a good correlation, with most st 
values for large amounts of storm rainfall being within 10%. 
_ The writer would like to inquire if the author compared his synthetic hydro- a 
_ graphs with those that would be computed by present techniques used by the a) 
S. Department of Agriculture. An examination of Fig. 9 and Fig. 5 indicates 
that the former is based on the occurrence of high intensity rainfall late in a 
storm, Perhaps the author ean state which | of the early or late, 
more commonly experienced in nature. 
_ The results of the author’s study are stimulating, Additional research ef- 
forts in other parts of the nation should be / encouraged to aid the hydrologist 
and the engineer in his design capacity. . Arecompilation and analysis of Ralston 
; Creek data has been completed at the University of Iowa, Iowa City, Iowa, and 
is to be published (1961) in cooperation with the lowa Highway Research 
a and other agencies. This report should also add to the field of knowledge 
L. McFALL!® and BENA. JONES, JR. 11 validity of the pro- 
_ posed method was tested with data from a central Illinois watershed main- 
by the Agricultural Experiment Station of the University of Illinois. 
: i is located near Monticello, approximately 120 miles northeast of 
Edwardsville, Ilinois. From the limited data available, two storms over the 
82-acre watershed were selected for comparison. 


Total Runoff, in inches | Excess © 


June 27, 951 0.7 — 0.51 | 0,58 
July 9, 1951] 1. 60,82 


The watershed is predominantly - silt loam soil and has slopes that are 
_ generally less than 4%, typical of the dark-colored, moderately permeable 
gently sloping prairie soils of east central Illinois. | 
, _ The time retention curves proposed by Minshall were used in computing 
runoff rates. The 30-day antecedent precipitation index (API) was computed 
. 3 the same formula used by Mr. Minshall except the K value was changed | 
‘to 0.90. (There is some question in the minds of the writers as to whether a 


80-day index should be used, but data are lacking to test this parameter.) 


Instr., Univ. of Dlinois, Urbana, Ml, 
Assoc, Prof., Univ. of Dlinois, Urbana, nl, 
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DABLE 6,—STUDY OF CENTRAL ILLINOIS WATERSHED 
_ Date of 30-day Peak Runoff Rates, in 
Storm | API | __ inches per hour _ _ 
observed | computed | observed [| computed tion, in 
: 
i 
| 


both storms was ‘less than the computed peak runoff. But the observed az. 
_ runoff was higher gone the eon value for one storm and the reverse was 


=a comparison tool than the peak discharge rate and the time-rate ‘change 
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. families of curves for each variation in soil retention characteristics. a 
data and retention curves presented were studied, and it appears feasible to # 
~ construct a coaxial correlation that would contain, in addition to the seasonal 
_ parameters, information about the soil characteristics. This proposed method _ 
will greatly facilitate the determination of precipitation retention when working _ 

a with widely varied soil conditions. More data and exploration are needed to © 
confirmthisopinton, 
ani The tests conducted indicate that the method proposed by Minshall for sy 
_ thesizing hydrographs will lead to peak runoff values that are higher than ob- © : 

ii served flows. It is impossible with the limited data available | to predict whether _ 

= total runoff will be higher or lower than observed values, In these cases the 

- difference in values does not exceed sass These are good results Fl 


{4 . ‘The results of the study are shown in Table 6 and one of the hydrographs 
a 
ae j c: _ the shape of the portion of the hydrograph that falls below the 15% runoff 4 | 
— 
| 
218 ‘aie from data from clay-pan prairie soils, the use of a soil retention factor,in 
y-pan p 
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the limited knowledge : about some of the factors, The w: ettens believe Minshall’s 


method shows promise and its applicability to a number of soi soil, cover, , and 
conditions should be investigated by other workers. 
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